Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,181,800 G→T 5.8% S255* (TCG→TAG)  thiG ← thiamine biosynthesis ThiGH complex subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,181,8000GT 5.8% 1503.7 / 44.8 389S255* (TCG→TAG) thiGthiamine biosynthesis ThiGH complex subunit
Reads supporting (aligned to +/- strand):  ref base G (169/197);  new base T (9/14);  total (178/211)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 391 total aligned reads displayed.

AGGCGATAACAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGC  >  CP009273/4181653‑4181948
                                                                                                                                                   |                                                                                                                                                    
   cGATAACAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGc                                                                                                                                                  <  1:5010058/149‑1 (MQ=255)
    gATAACAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGcc                                                                                                                                                 <  1:2206923/149‑1 (MQ=255)
     aTAACAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCt                                                                                                                                                <  2:351153/149‑1 (MQ=255)
      tAACAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTc                                                                                                                                               >  2:1321912/1‑149 (MQ=255)
        aCAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCa                                                                                                                                             >  1:4471767/1‑149 (MQ=255)
        aCAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCa                                                                                                                                             <  2:3716297/149‑1 (MQ=255)
             gcCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAAt                                                                                                                                        >  1:5072919/1‑149 (MQ=255)
              ccATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATc                                                                                                                                       <  1:173074/149‑1 (MQ=255)
              ccATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATc                                                                                                                                       >  2:3799324/1‑149 (MQ=255)
              ccATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATc                                                                                                                                       <  2:3902249/149‑1 (MQ=255)
                 tcatGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCgg                                                                                                                                    <  1:2855110/149‑1 (MQ=255)
                     gtcgtcGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAg                                                                                                                                <  2:2484560/149‑1 (MQ=255)
                      tcgtcgCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGc                                                                                                                               <  1:211288/149‑1 (MQ=255)
                        gtcgCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCgg                                                                                                                             <  2:4329490/149‑1 (MQ=255)
                          cgcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGc                                                                                                                           <  1:2471139/149‑1 (MQ=255)
                           gcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCt                                                                                                                          >  1:2279451/1‑149 (MQ=255)
                           gcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCt                                                                                                                          <  1:1059491/149‑1 (MQ=255)
                           gcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCt                                                                                                                          <  1:2526312/149‑1 (MQ=255)
                           gcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCt                                                                                                                          >  1:2634429/1‑149 (MQ=255)
                              ggTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtggt                                                                                                                       <  2:3583246/149‑1 (MQ=255)
                              ggTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtggt                                                                                                                       >  1:1297946/1‑149 (MQ=255)
                                tGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtggtgg                                                                                                                     <  1:3440992/149‑1 (MQ=255)
                                  aGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCTATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCa                                                                                                                   <  1:4194053/149‑1 (MQ=255)
                                  aGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCa                                                                                                                   <  1:1791937/149‑1 (MQ=255)
                                      gCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAg                                                                                                               >  2:5242590/1‑149 (MQ=255)
                                        gAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCa                                                                                                             >  1:2773875/1‑149 (MQ=255)
                                          ggCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCaaa                                                                                                           >  1:3604744/1‑149 (MQ=255)
                                            cATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAAt                                                                                                         <  1:1647962/149‑1 (MQ=255)
                                             aTTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGCTCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATa                                                                                                        >  1:586293/1‑149 (MQ=255)
                                             aTTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATa                                                                                                        <  1:3218267/149‑1 (MQ=255)
                                              tttAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATaa                                                                                                       >  1:3054339/1‑149 (MQ=255)
                                              tttAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATaa                                                                                                       <  1:2755076/149‑1 (MQ=255)
                                              tttAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATaa                                                                                                       <  1:3915311/149‑1 (MQ=255)
                                                tAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACt                                                                                                     >  2:1775457/1‑149 (MQ=255)
                                                 aGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTg                                                                                                    <  2:1710505/149‑1 (MQ=255)
                                                  gcgcCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCTATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGc                                                                                                   <  2:260582/149‑1 (MQ=255)
                                                  gcgcCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGc                                                                                                   <  2:606148/149‑1 (MQ=255)
                                                  gcgcCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGc                                                                                                   <  2:3250565/149‑1 (MQ=255)
                                                   cgcCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCg                                                                                                  <  1:1035113/149‑1 (MQ=255)
                                                      ccGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCt                                                                                               >  2:3185653/1‑149 (MQ=255)
                                                          tctACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGccc                                                                                           >  1:2021433/1‑149 (MQ=255)
                                                               gTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTcc                                                                                      <  2:4727197/149‑1 (MQ=255)
                                                                 cagcagCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCTATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCgg                                                                                    <  1:481921/149‑1 (MQ=255)
                                                                  agcagcCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGa                                                                                   >  2:2949890/1‑149 (MQ=255)
                                                                  agcagcCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGa                                                                                   >  2:2701983/1‑149 (MQ=255)
                                                                       ccGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGAc                                                                              >  2:1618344/1‑149 (MQ=255)
                                                                         gTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGt                                                                            >  1:412556/1‑149 (MQ=255)
                                                                          ttttGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTg                                                                           >  2:2783388/1‑149 (MQ=255)
                                                                            ttGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGcc                                                                         <  1:3099653/149‑1 (MQ=255)
                                                                             tGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCa                                                                        <  2:4481637/149‑1 (MQ=255)
                                                                               ccGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGt                                                                      >  1:4509711/1‑149 (MQ=255)
                                                                               ccGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGt                                                                      <  1:1583880/149‑1 (MQ=255)
                                                                                cGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCCCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTa                                                                     >  2:4051621/1‑149 (MQ=255)
                                                                                 gTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAg                                                                    >  1:3188204/1‑149 (MQ=255)
                                                                                    gATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGcc                                                                 <  2:3365960/149‑1 (MQ=255)
                                                                                      tatgCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTg                                                               <  2:807172/146‑1 (MQ=255)
                                                                                       aCGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGc                                                              <  2:1466083/149‑1 (MQ=255)
                                                                                       aCGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGc                                                              <  1:1908424/149‑1 (MQ=255)
                                                                                          cAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTc                                                           >  2:2440459/1‑149 (MQ=255)
                                                                                            ggCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTa                                                         >  2:3045849/1‑149 (MQ=255)
                                                                                                 aTGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGcc                                                    <  2:4380088/149‑1 (MQ=255)
                                                                                                     cgtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGAc                                                <  1:686001/149‑1 (MQ=255)
                                                                                                      gtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACg                                               >  2:1008261/1‑149 (MQ=255)
                                                                                                      gtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACg                                               >  2:4900592/1‑149 (MQ=255)
                                                                                                        cccAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCTATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGaa                                             <  1:2327001/149‑1 (MQ=255)
                                                                                                        cccAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGTGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGaa                                             >  2:2118109/1‑149 (MQ=255)
                                                                                                         ccagtccagtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGaaa                                            <  2:5092555/149‑1 (MQ=255)
                                                                                                         ccagtccagtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGaaa                                            <  1:856478/149‑1 (MQ=255)
                                                                                                            gtccagtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGc                                         >  2:1998525/1‑149 (MQ=255)
                                                                                                                agtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTc                                     <  2:2373537/149‑1 (MQ=255)
                                                                                                                  ttGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGc                                   <  1:4745711/149‑1 (MQ=255)
                                                                                                                   tGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGcc                                  <  2:4873728/149‑1 (MQ=255)
                                                                                                                   tGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGcc                                  <  2:2248647/149‑1 (MQ=255)
                                                                                                                    gTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCa                                 >  2:1144302/1‑149 (MQ=255)
                                                                                                                     tCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCAt                                <  2:2764193/149‑1 (MQ=255)
                                                                                                                     tCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCAt                                <  2:3776851/149‑1 (MQ=255)
                                                                                                                          aGCGATCGCTGAAGGTTTTCATGCCTATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGa                           <  1:2968398/149‑1 (MQ=255)
                                                                                                                          aGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGa                           >  1:1388007/1‑149 (MQ=255)
                                                                                                                          aGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGa                           >  2:4255334/1‑149 (MQ=255)
                                                                                                                           gCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGAc                          >  1:249735/1‑149 (MQ=255)
                                                                                                                              aTCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACggg                       >  1:1225802/1‑149 (MQ=255)
                                                                                                                              aTCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACggg                       <  1:5001426/149‑1 (MQ=255)
                                                                                                                               tCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGa                      <  1:4315264/149‑1 (MQ=255)
                                                                                                                                cGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGAt                     <  1:1842648/149‑1 (MQ=255)
                                                                                                                                cGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGAt                     >  2:1219182/1‑149 (MQ=255)
                                                                                                                                   tGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGt                  >  2:1392344/1‑149 (MQ=255)
                                                                                                                                    gAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTc                 <  2:417198/149‑1 (MQ=255)
                                                                                                                                    gAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTc                 <  1:4408648/149‑1 (MQ=255)
                                                                                                                                     aaGGTTTTCATGCCTATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTcc                <  2:5002425/149‑1 (MQ=255)
                                                                                                                                      aGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCcg               >  2:2638609/1‑149 (MQ=255)
                                                                                                                                         ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGa            >  2:4394314/1‑149 (MQ=255)
                                                                                                                                         ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGa            <  1:1749602/149‑1 (MQ=255)
                                                                                                                                         ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGa            <  2:3080736/149‑1 (MQ=255)
                                                                                                                                             cATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGc        >  2:3828583/1‑149 (MQ=255)
                                                                                                                                              aTGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCa       >  2:1091604/1‑149 (MQ=255)
                                                                                                                                              aTGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCa       <  2:2055470/149‑1 (MQ=255)
                                                                                                                                              aTGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCa       >  1:1024966/1‑149 (MQ=255)
                                                                                                                                               tGCCGCTGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCaa      <  1:5101140/149‑1 (MQ=255)
                                                                                                                                                gCCGATGCCTCCAGAAATCCGGTCAGCTGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAAt     >  1:3723592/1‑149 (MQ=255)
                                                                                                                                                   gATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGc  >  2:725936/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGCGATAACAGCGCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGC  >  CP009273/4181653‑4181948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: