Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,222,6920GC25.4% 492.3 / 90.0 213intergenic (+148/+225)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base G (52/107);  new base C (54/0);  total (106/107)
Fisher's exact test for biased strand distribution p-value = 4.15e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

AAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAT  >  CP009273/2222543‑2222841
                                                                                                                                                     |                                                                                                                                                     
aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                                                       <  2:3859666/150‑1 (MQ=255)
 aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                                                                      <  2:2928616/150‑1 (MQ=255)
   attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGc                                                                                                                                                    <  2:658448/150‑1 (MQ=255)
   attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCAAAGCCGCCCCCCCCCGCACCCCCCGGc                                                                                                                                                    >  1:4839136/1‑150 (MQ=38)
    ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                                                   <  1:4526147/150‑1 (MQ=255)
    ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCTGCCCCCCCGGGct                                                                                                                                                   >  2:2981716/1‑150 (MQ=255)
    ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGct                                                                                                                                                   >  2:2915215/1‑150 (MQ=255)
      aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                                                                 <  1:2718898/150‑1 (MQ=255)
       attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                                                                >  2:3083586/1‑150 (MQ=255)
       attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                                                                <  2:3976666/150‑1 (MQ=255)
         taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                                                              <  2:468987/150‑1 (MQ=255)
           atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCACTGGCGa                                                                                                                                            >  1:2731444/1‑150 (MQ=255)
            tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGaa                                                                                                                                           >  2:2758730/1‑150 (MQ=255)
            tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAATCACCCCCCCCCACCACCCCCCCGCACAGGCGaa                                                                                                                                           >  1:1817825/1‑150 (MQ=37)
             gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGATGGGGCGGCCCCACGCCGCCGCCCCCCGCCCCCCCCGGCACTGGCGAAc                                                                                                                                          >  2:2477690/1‑150 (MQ=37)
             gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAAc                                                                                                                                          >  2:3230074/1‑150 (MQ=255)
                 tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACGCCCCCGCCCCCAGCACCCCCCGGCGCGGGCGGACaaaa                                                                                                                                      >  1:1469306/1‑150 (MQ=25)
                   aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCGGGCTCATGCGGACaaaaaa                                                                                                                                    >  1:3063363/1‑148 (MQ=255)
                    ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAAAAAAATcc                                                                                                                                   >  2:3036379/1‑149 (MQ=255)
                     gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAAAcgc                                                                                                                                  >  1:4269072/1‑150 (MQ=25)
                      tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                                                                 <  1:4679672/150‑1 (MQ=255)
                      tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                                                                 <  2:1795322/150‑1 (MQ=255)
                      tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                                                                 <  2:2751038/150‑1 (MQ=255)
                      tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                                                                 <  2:300750/150‑1 (MQ=255)
                       ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccc                                                                                                                                >  2:547784/1‑149 (MQ=255)
                       ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccg                                                                                                                                <  1:972475/150‑1 (MQ=255)
                       ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCCCCCCCAGCACCCCCCCGCGCAGGCGAACaaaaaagccc                                                                                                                                >  1:2292194/1‑144 (MQ=25)
                         gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCt                                                                                                                              <  2:372251/150‑1 (MQ=255)
                          aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                                                             <  1:3934678/150‑1 (MQ=255)
                            aaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcg                                                                                                                           <  2:3541462/150‑1 (MQ=255)
                             aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGCCAATAGCCGCCGCCCCCAGCACCCCCCGGCTCAGGCGCACAAAAACCCCCCGgcgc                                                                                                                          >  1:4206824/1‑150 (MQ=21)
                                aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                                                       >  2:4468209/1‑150 (MQ=255)
                                   agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAAACCCCCCCCCCCCAGCACCCCCCGGCTCAGGCGAACAAAAAAGCCCCAGCGCGGTGcc                                                                                                                    >  2:2130514/1‑150 (MQ=255)
                                   agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCCGCCAATCGCCGCCGCCCCCAGCAACCCCGGGCACAGGAGAACAAAATCACCGCGGCGCCGTGcc                                                                                                                    >  2:1113554/1‑150 (MQ=21)
                                    gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGccc                                                                                                                   <  1:3689758/150‑1 (MQ=255)
                                       tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCCCGCTCAGGCGAACAAAAAACCCCCAGCGCGGTGCCCACg                                                                                                                >  1:4807915/1‑150 (MQ=18)
                                        atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                                                               <  2:4089139/150‑1 (MQ=255)
                                         tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCGCCCCCAGCACCCCCCGGCACAGGCGAACAAAATCCCCGCAGCGCGGTGCCCACGGc                                                                                                              >  1:390063/1‑150 (MQ=18)
                                           tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGGACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                                                            <  1:910989/150‑1 (MQ=255)
                                            cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                                                           <  2:697336/150‑1 (MQ=255)
                                             aGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                                                          <  2:89902/150‑1 (MQ=255)
                                                  ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                                                     <  2:2881718/150‑1 (MQ=255)
                                                   tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATAGCCGCCCCCCCCAGCACCCCCCGGCACAGGCGAACAAAATCGCCCCGGCGCGGTGCCCACGGCGTATACCATa                                                                                                    >  1:348699/1‑150 (MQ=18)
                                                    cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATACTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                                                                   <  2:4366170/150‑1 (MQ=255)
                                                    cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAAGCCGCCCCCCCCAGCAACCCCCGGCTCAGGCGAAAAAAAAAGCCCCGGCGCGGTGCCCACGGCATATACCATaa                                                                                                   >  1:1987845/1‑149 (MQ=18)
                                                     aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                                                                  >  2:2409593/1‑150 (MQ=255)
                                                         tggtGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTATGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                                                              <  1:354738/147‑1 (MQ=255)
                                                         tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                                                              <  2:447514/150‑1 (MQ=255)
                                                         tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                                                              <  2:2436219/150‑1 (MQ=255)
                                                           ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                                                            <  1:3986244/150‑1 (MQ=255)
                                                              cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAAAGCCGCCGCCCCCAGCACCCCCCGGCTCTTGCGAACAAAATAGCCGCAGCGCGGGGGCCACGGCGTATCACATAAGGCTACCcgg                                                                                         >  1:4030952/1‑150 (MQ=17)
                                                                cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                                                       <  2:782579/150‑1 (MQ=255)
                                                                 tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAAGCCCCCCCCCCCAGCACCCCCCGGCTAATGGGAACAAAAAACCCGCGGCGCGGTGGCCAACGCGTATACCATAAGGCAACCCGGTCg                                                                                      >  1:1481636/1‑150 (MQ=12)
                                                                 tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCg                                                                                      <  2:426100/150‑1 (MQ=255)
                                                                   cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                                                    <  1:4313430/150‑1 (MQ=255)
                                                                    tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgc                                                                                   <  1:3000325/150‑1 (MQ=255)
                                                                      tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCCCCCCCCGCACCCCCGGGCTCTGGCGAACAAAAACCCCGCTGCGCGGTGCCCTCGGCTTATCCCATACGGCTACCGGGTCGGgcgc                                                                                 >  2:4254525/1‑150 (MQ=17)
                                                                          gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                                                             <  1:4122137/150‑1 (MQ=255)
                                                                             cATTGCTTCTTATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCATTACGGCTAACGGGTCGGGCGCGAGGTaa                                                                          <  1:1258845/150‑1 (MQ=255)
                                                                             cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                                                          <  2:2153149/150‑1 (MQ=255)
                                                                              aTTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGCCCAACGCCGCCCCCCCCAGCACCCCCCGGCGCAGGCGAACAAAAAACCCCCAGCGCGGTGCCCACGGCTTATACCAAACCGCAAACCCGTCCGGCGCGAGGGAAc                                                                         >  1:2835158/1‑150 (MQ=255)
                                                                                tGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                                       <  1:2794672/150‑1 (MQ=255)
                                                                                tGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGCTCATGCGAACAAAAACCCCCCTGCCCGGGGCCCACGGCTTATACCATACGGCTACCCGGTCGGGCGCGCGGTAACaa                                                                       >  2:3780318/1‑149 (MQ=14)
                                                                                 gCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCTGCACCCCCCGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGCGTATACCTTACGGCTACCGGGTCGGGCGCGAGGTAACaac                                                                      >  2:2676122/1‑148 (MQ=14)
                                                                                    tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                                                                   <  1:350276/150‑1 (MQ=39)
                                                                                           aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                                                            <  2:875491/150‑1 (MQ=33)
                                                                                           aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                                                            <  1:1505615/150‑1 (MQ=33)
                                                                                                tGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                                                       <  1:2193155/150‑1 (MQ=32)
                                                                                                  gcgcAAGGAGCGCGGAGGGGGGGGCCAAAAGCCCCCGCCCCCAGCCCCCCCCGGCTCTGGCGAACAAAAAAGCCCCAGGGCGGTGCCCTCCGCGTATCCCATAAGGCGACCCGGTCGGGCGCGAGGTAACAACCCAGTAAAAACCGccca                                                     >  2:4588017/1‑149 (MQ=2)
                                                                                                      aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc                                                 <  1:81249/150‑1 (MQ=32)
                                                                                                             gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc                                          <  1:423616/150‑1 (MQ=32)
                                                                                                                 aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg                                      >  1:3453005/1‑150 (MQ=32)
                                                                                                                 aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg                                      <  2:756849/150‑1 (MQ=32)
                                                                                                                  gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                                     >  1:1347932/1‑150 (MQ=32)
                                                                                                                   ggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                                    >  1:1447204/1‑150 (MQ=32)
                                                                                                                    gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg                                   >  1:3815020/1‑150 (MQ=32)
                                                                                                                      gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGc                                 <  2:1962750/150‑1 (MQ=32)
                                                                                                                      gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGc                                 <  2:4390839/150‑1 (MQ=32)
                                                                                                                       cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                                <  1:489742/150‑1 (MQ=32)
                                                                                                                       cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCt                                <  1:2398609/150‑1 (MQ=32)
                                                                                                                        ggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg                               >  1:4233102/1‑150 (MQ=32)
                                                                                                                          ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcc                             <  1:275356/150‑1 (MQ=32)
                                                                                                                           cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                            <  1:1049884/150‑1 (MQ=32)
                                                                                                                            aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                           >  1:2437325/1‑150 (MQ=32)
                                                                                                                            aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                           <  1:2142185/150‑1 (MQ=32)
                                                                                                                            aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                           >  2:2218704/1‑150 (MQ=32)
                                                                                                                            aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                           <  2:3497684/150‑1 (MQ=32)
                                                                                                                             aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                          >  1:873055/1‑150 (MQ=32)
                                                                                                                                gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGcc                       <  2:1495383/150‑1 (MQ=32)
                                                                                                                                        ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa               <  2:465128/150‑1 (MQ=32)
                                                                                                                                          ccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc             <  1:1487509/150‑1 (MQ=32)
                                                                                                                                             gCACCCCCGGGCTCTGGCGAACACAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc          <  2:2855136/150‑1 (MQ=17)
                                                                                                                                             gCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc          >  1:59393/1‑150 (MQ=32)
                                                                                                                                              cACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCt         <  2:1381375/150‑1 (MQ=32)
                                                                                                                                              cACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCt         <  2:574320/150‑1 (MQ=32)
                                                                                                                                              cACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCt         <  2:1469306/150‑1 (MQ=32)
                                                                                                                                                 ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCgg      >  1:1679879/1‑150 (MQ=32)
                                                                                                                                                 ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCgg      >  2:3293747/1‑150 (MQ=32)
                                                                                                                                                    cGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCaa   <  2:507860/150‑1 (MQ=33)
                                                                                                                                                     gggCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTACCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTGCCGGGAACGCCTCGGCAAt  >  2:3212242/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAT  >  CP009273/2222543‑2222841

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