Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,406 T→C 19.0% Y275Y (TAT→TAC)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,4060T→C19.0% 150.2 ‑5.3 ‑0.0268Y275Y (TAT→TAC) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (15/36):  ref base (111/106):  total (126/142)
Fisher's exact test strand distribution p-value = 5.20e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.29e-01

Only 100 of 289 total aligned reads displayed.

ACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAG  >  NC_007779/3441306‑3441506
                                                                                                    |                                                                                                    
aCTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtat                                                                                                      >  1:255727/1‑101 (MQ=255)
 cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                                                                     >  1:3376171/1‑101 (MQ=255)
 cTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCtata                                                                                                     >  1:9326167/1‑101 (MQ=255)
  tCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTatat                                                                                                    >  1:11546967/1‑101 (MQ=255)
   cggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATAtt                                                                                                   <  1:14330942/101‑1 (MQ=255)
   cggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATAtt                                                                                                   >  1:6474953/1‑101 (MQ=255)
    ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                                                  >  1:14840825/1‑101 (MQ=255)
    ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                                                  >  1:1124094/1‑101 (MQ=255)
    ggcggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTc                                                                                                  <  1:13447622/101‑1 (MQ=255)
       ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                                                                               >  1:4733432/1‑101 (MQ=255)
       ggTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACg                                                                                               <  1:12274079/101‑1 (MQ=255)
          cACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGaaa                                                                                            >  1:5340291/1‑101 (MQ=255)
           aCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAc                                                                                           >  1:3045128/1‑101 (MQ=255)
           aCCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAAc                                                                                           >  1:5502321/1‑101 (MQ=255)
              tctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcg                                                                                        <  1:8857959/101‑1 (MQ=255)
              tctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcg                                                                                        >  1:13871831/1‑101 (MQ=255)
              tctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcg                                                                                        <  1:4576164/101‑1 (MQ=255)
               ctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgc                                                                                       >  1:10430107/1‑101 (MQ=255)
               ctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgc                                                                                       <  1:7646994/101‑1 (MQ=255)
               ctTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgc                                                                                       >  1:9655142/1‑101 (MQ=255)
                ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                                                                      <  1:7206490/101‑1 (MQ=255)
                ttCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCcgcg                                                                                      <  1:11845316/101‑1 (MQ=255)
                 tCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGa                                                                                     >  1:4021795/1‑101 (MQ=255)
                  ccATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAg                                                                                    >  1:13655026/1‑101 (MQ=255)
                     tctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTgg                                                                                 <  1:4985198/96‑1 (MQ=255)
                      ctccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGa                                                                                >  1:1308234/5‑101 (MQ=255)
                      cGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGa                                                                                >  1:8314000/1‑101 (MQ=255)
                       tccgaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAt                                                                               >  1:5551260/4‑101 (MQ=255)
                         agaagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGAttt                                                                             >  1:541507/1‑101 (MQ=255)
                          gaagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGGCACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGAttct                                                                            >  1:15111344/1‑99 (MQ=255)
                           aagaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGACCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTg                                                                           >  1:10880419/1‑101 (MQ=255)
                           aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                                                                           >  1:1205864/1‑101 (MQ=255)
                           aagaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTg                                                                           >  1:10717848/1‑101 (MQ=255)
                             gaggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCg                                                                         <  1:13056218/98‑1 (MQ=255)
                             gaagTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCg                                                                         >  1:14720665/1‑101 (MQ=255)
                              aggTGGAGAAGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGt                                                                        <  1:7345039/99‑1 (MQ=255)
                               agTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTa                                                                       <  1:9833529/101‑1 (MQ=255)
                                 tGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAAc                                                                     <  1:6709014/101‑1 (MQ=255)
                                  ggAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACa                                                                    <  1:7145823/101‑1 (MQ=255)
                                   gagaaGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAActc                                                                   <  1:5798036/98‑3 (MQ=255)
                                   gAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAg                                                                   <  1:1687633/101‑1 (MQ=255)
                                   gAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAg                                                                   >  1:13334781/1‑101 (MQ=255)
                                     gaaGATGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAActccc                                                                 <  1:10324190/100‑5 (MQ=255)
                                     aaaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGcc                                                                 <  1:6995050/101‑1 (MQ=255)
                                      aaGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGccc                                                                <  1:2385759/101‑1 (MQ=255)
                                       agatGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCg                                                               <  1:613476/98‑1 (MQ=255)
                                       aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg                                                               <  1:7842378/101‑1 (MQ=255)
                                       aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCg                                                               >  1:14544600/1‑101 (MQ=255)
                                       aGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCg                                                               >  1:9442750/1‑101 (MQ=255)
                                         cTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGt                                                             <  1:12428216/101‑1 (MQ=255)
                                         cTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGt                                                             >  1:3558040/1‑101 (MQ=255)
                                         atGGTGTGGGCGACCCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGt                                                             <  1:11117541/100‑1 (MQ=255)
                                          tGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTg                                                            >  1:2093468/1‑101 (MQ=255)
                                           gttATGGTCCCCGCCCTGGGGGCCGGTTCCGTGGTTGGTCCCCCCCCCCGGTCCCTTTA‑TCCCCAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGcc                                                          <  1:4377175/99‑1 (MQ=255)
                                           ggTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGc                                                           <  1:14121116/101‑1 (MQ=255)
                                             tATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCg                                                         >  1:11368602/1‑101 (MQ=255)
                                              gtgggcgaCGCGTTGGGGCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCgg                                                        >  1:10119694/8‑100 (MQ=255)
                                               tGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGAt                                                       <  1:12783668/101‑1 (MQ=255)
                                                ggTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATc                                                      >  1:3895185/1‑101 (MQ=255)
                                                  tCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCgg                                                    >  1:12933103/1‑101 (MQ=255)
                                                    cACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTa                                                  >  1:7818437/1‑101 (MQ=255)
                                                     aCGTGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTAc                                                 <  1:955583/101‑1 (MQ=255)
                                                      cgcgCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACa                                                >  1:2841222/1‑101 (MQ=255)
                                                       gcgcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAg                                               >  1:14676753/1‑101 (MQ=255)
                                                       gcgcTGGGGAGCGGATACGGTGACGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAg                                               <  1:1891356/101‑1 (MQ=255)
                                                        cgcTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGt                                              <  1:2582497/101‑1 (MQ=255)
                                                          cTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGc                                            <  1:10935594/101‑1 (MQ=255)
                                                           tGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGcc                                           <  1:5955996/101‑1 (MQ=255)
                                                            ggggCGCGGACACCATCATGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTGCCg                                          <  1:3204542/101‑1 (MQ=255)
                                                                aGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCt                                      <  1:4412078/101‑1 (MQ=255)
                                                                  cGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAc                                    <  1:12241109/101‑1 (MQ=255)
                                                                   ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAcc                                   >  1:11021705/1‑101 (MQ=255)
                                                                   ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAcc                                   >  1:927541/1‑101 (MQ=255)
                                                                   ggATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAcc                                   >  1:6333413/1‑101 (MQ=255)
                                                                    gATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCa                                  >  1:8245959/1‑101 (MQ=255)
                                                                    gATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCa                                  >  1:1209936/1‑101 (MQ=255)
                                                                    gATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCa                                  >  1:9216994/1‑101 (MQ=255)
                                                                     aTACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAg                                 <  1:2686515/101‑1 (MQ=255)
                                                                      tACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAgg                                >  1:3090999/1‑101 (MQ=255)
                                                                       accatcaTGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTGCCGATTTACCAGgc                               <  1:1320867/95‑1 (MQ=255)
                                                                        cGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGgcg                              >  1:4920291/1‑101 (MQ=255)
                                                                           tcaTGGACCTGTCCACCGGTAAGGACATTCACGAAACCCGCGAGTGGATTCTGCGTAACTCCCCGGTGCCGGTTGGTACGGTGCCGATTTACCAGGCGCTg                           <  1:174100/99‑1 (MQ=255)
                                                                           tGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTg                           <  1:1403768/101‑1 (MQ=255)
                                                                                  tctctcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAgg                    >  1:3003495/1‑101 (MQ=255)
                                                                                      tcCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTaa                >  1:15049794/1‑101 (MQ=255)
                                                                                       ccACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAAc               <  1:12398925/101‑1 (MQ=255)
                                                                                          ccGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACggg            >  1:657643/1‑101 (MQ=255)
                                                                                           cGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGa           >  1:8973061/1‑101 (MQ=255)
                                                                                           cGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGa           <  1:3748014/101‑1 (MQ=255)
                                                                                             gTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATc         <  1:877773/101‑1 (MQ=255)
                                                                                             gTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATc         >  1:5651178/1‑101 (MQ=255)
                                                                                             gTCGCTATATTCACGAAACCCGCGAGTGGATTGTGCGTAACAGCCCGGTGGCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATc         <  1:6734103/101‑1 (MQ=255)
                                                                                              tCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCg        >  1:3512726/1‑101 (MQ=255)
                                                                                               cGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGc       <  1:5247993/101‑1 (MQ=255)
                                                                                                gCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTAGCAGGCGCTGGAGAAGGTTAACGGGATCGcc      >  1:9071981/1‑101 (MQ=255)
                                                                                                gCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGACCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGcc      >  1:6859807/1‑101 (MQ=255)
                                                                                                 cTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCg     >  1:2317146/1‑101 (MQ=255)
                                                                                                  tataTTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGa    <  1:11953785/101‑1 (MQ=255)
                                                                                                  tataTTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGa    <  1:1560099/101‑1 (MQ=255)
                                                                                                  tataTTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGa    >  1:10051288/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
ACTCGGCGGTCACCTCTTCCATCGAAGAAGAAGTGGAAAAGCTGGTATGGTCCACGCGCTGGGGAGCGGATACGGTGATGGATCTCTCCACCGGTCGCTATATTCACGAAACCCGCGAGTGGATTTTGCGTAACAGCCCGGTGCCGATCGGTACAGTGCCGATCTACCAGGCGCTGGAGAAGGTTAACGGGATCGCCGAAG  >  NC_007779/3441306‑3441506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.