Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,156 C→G 6.9% intergenic (‑90/+505) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1560CG 6.9% 738.7 / 10.1 250intergenic (‑90/+505)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (55/175);  new base G (13/7);  total (68/182)
Fisher's exact test for biased strand distribution p-value = 2.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360008‑360302
                                                                                                                                                     |                                                                                                                                                  
 gCACTTCACCGCTTGCCAGCGGCTTACCATCCCGCGCCACCATCCAGTGCTGGAGCTCCTATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCGGGCAGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                  <  1:3230948/150‑1 (MQ=38)
  ctactcaggagagcgttcacccacaaacaacccataaaacgaacgccccagtctttcgactgagcctttctttttatttgatgcctggggaattAATACCTGGAATTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                                 >  2:541034‑M1/95‑150 (MQ=255)
   tactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                <  2:2221220‑M1/57‑1 (MQ=255)
   tactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                <  2:46034‑M1/57‑1 (MQ=255)
      tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCaa                                                                                                                                             >  1:868165/1‑150 (MQ=255)
      tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCaa                                                                                                                                             >  1:4272798/1‑150 (MQ=255)
       caggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                            <  1:2497959‑M1/61‑1 (MQ=255)
       cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                            >  2:4057812/1‑150 (MQ=255)
          cggttGCCAGCGGCTTACCACCCCGCGCCACCATCCAGGGCAGGAGCTCCTATTATAAGCTGTAAAACATGATAATTCCGACCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         <  1:3862973/147‑1 (MQ=25)
           gCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                        <  2:4242769/150‑1 (MQ=255)
           gCTTGCCAGCGGCTTACCACCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGAAAAACATGAGAATTCCGATCCCCAATCCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                        <  2:2545200/150‑1 (MQ=37)
            cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCCACCTTcc                                                                                                                                       >  1:3275158/1‑149 (MQ=255)
             ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  2:2771984/150‑1 (MQ=255)
               gCCAGGGGCTTACCATCCCGCGCCACCATCCCGGGCTGGAGCTCCGTGGATACGCTGTCAAACATGAGAATGCCGATCCCCAAGTCCTGGATTTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                    <  1:2370998/150‑1 (MQ=21)
                 cAGCGGATTACCATCCCGCGCCACCCTCCTGTGCTGGAGCTCCGATTATAAGGTGACAAACATCAGAATCCCGACCCCCACTGCCGGGCTGTTTATGGGGGGCGACCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                  <  1:945466/150‑1 (MQ=18)
                  aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                 <  2:1695711/150‑1 (MQ=255)
                     ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGGTAAAATACCCTc                                                                                                                              >  1:3483453/1‑150 (MQ=37)
                      cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATCCGCTcc                                                                                                                             <  2:1453350‑M1/76‑1 (MQ=255)
                         tACCCTCCAGCGCCACCATCCTGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                          <  1:3629873/150‑1 (MQ=255)
                         caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCAACGTTCAaatacgcagcccg                                                                                                                          <  1:1181588‑M1/79‑7 (MQ=255)
                         caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTGGCCACCTTCAAAACCGCTCCcgg                                                                                                                          <  2:67226‑M1/79‑1 (MQ=255)
                           aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGGAACGTTCAAAGCCGGTCgcggcg                                                                                                                        <  2:4265703‑M1/81‑1 (MQ=255)
                             caacagataaaacgaaaggccccgtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCGTCGCACCGTTCAACTCCGCTCccggcggc                                                                                                                      <  2:3032383‑M1/83‑2 (MQ=255)
                             aTCCCGCGCCCCCAGCCTGTGCTGTAGCTCCGATGATAAGCTTTCAAACATGACATTTCCGACCCCCATTGCGTGGATGTTTTGGGGGGGGGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  1:1868303/150‑1 (MQ=12)
                              tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                     <  2:4211995/150‑1 (MQ=255)
                              aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                     <  2:1766945‑M1/84‑1 (MQ=255)
                                    ccccccaTCCCGGGCCGGAGCTCCGATTATAAGCTGACAAACATTATAATTCCGACCCCAATTTCCTGGCTGTTTCTGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                               <  1:3424279/146‑1 (MQ=17)
                                     ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                              >  1:3585621/1‑150 (MQ=255)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCAACCGTCAAATAAGCACCCCGCGGCtttttgcaa                                                                                                             >  1:2679555/1‑144 (MQ=25)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACCGTCAAATACGCGCCCCGCGGATTTGtacaa                                                                                                             >  1:2796548/1‑146 (MQ=25)
                                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                            <  1:4297579/150‑1 (MQ=255)
                                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCCGTTCCGGTCCCGGCGGAGTTGTCCTAc                                                                                                            <  2:4178909‑M1/93‑1 (MQ=255)
                                          tccccgTGCTGGTGCTCCGTTGATAAGCTGACAAACTTGATATTTCCGAACCCAAATTCCGGGCTTTTTTTGGGGGGGGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                         <  1:3360574/145‑1 (MQ=17)
                                          gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattATTTCCTGTCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATCCGCTCCCGGCGGCTTTGTCCTCCTCa                                                                                                         <  2:798844‑M1/96‑1 (MQ=255)
                                          aTCCTGTGCTGGAGCTCCGATGATATGCTGACAAACATGAGAATTCCGACCCCCAATTCCTGGCAGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                         <  1:2318901/150‑1 (MQ=18)
                                            aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                       <  2:3978897‑M1/98‑1 (MQ=255)
                                            aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTACTACTCAgg                                                                                                       <  2:3337770‑M1/98‑1 (MQ=255)
                                              ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                 >  2:2499876‑M1/51‑122 (MQ=255)
                                              ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGAGTTGTCCTACTcaggcg                                                                                                     <  2:4021669‑M1/100‑3 (MQ=255)
                                              cgTGCCGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                     <  2:3611716/149‑1 (MQ=255)
                                                tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCCACCTTCAAATACCCAACCCGCGGATTTGTCCGACACACgagat                                                                                                   >  1:4178292/1‑149 (MQ=21)
                                                 ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCATAAAAAACGCTCCCCGCCGCTTTGTGCAACACAGGCGAGc                                                                                                  <  1:4155751‑M1/103‑1 (MQ=255)
                                                  cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                 >  2:3012410/1‑150 (MQ=255)
                                                   agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                <  2:4260217‑M1/105‑1 (MQ=255)
                                                   aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                <  2:3544064/150‑1 (MQ=255)
                                                    ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                               <  2:1972622/150‑1 (MQ=255)
                                                      ctttcgactgagcctttggttttatttgatgcgtggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             >  1:2134028‑M1/43‑150 (MQ=255)
                                                      ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGGTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:2306637‑M1/108‑1 (MQ=255)
                                                      aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:3813323/150‑1 (MQ=255)
                                                       gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCATTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                            <  1:1322095/150‑1 (MQ=255)
                                                              tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACCTTAAAATACGGACCCGGCGGCTTTTTGGTAAACAGGTGGGCGTTAACACacaaa                                                                                     >  2:370003/1‑150 (MQ=18)
                                                               tatAAGCTGTAAAACAGGACAATTCCTATCCCAAATACGTGGCATTTTATGGCGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                    <  1:456309/149‑1 (MQ=14)
                                                                agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                   <  2:2504106‑M1/118‑1 (MQ=255)
                                                                  cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCAAaacaac                                                                                 <  2:1570209‑M1/120‑1 (MQ=255)
                                                                    tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                               <  2:4285960‑M1/122‑1 (MQ=255)
                                                                    gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                               <  2:92621/150‑1 (MQ=255)
                                                                    gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                               <  1:1674890/150‑1 (MQ=255)
                                                                     ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTCCTACACAGGAGTGCGTTCACCGACAAACCACAGa                                                                              <  1:3995563‑M1/123‑1 (MQ=255)
                                                                        tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                           <  2:3915038/150‑1 (MQ=255)
                                                                         ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACCAACCCCAGATaaa                                                                          <  2:3705473‑M1/127‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                        <  2:3823060/150‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                        <  2:2671572/150‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                        <  2:242535/150‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAAACTTAAAAAACGCGACCCGCGGAGTTTTGCAACTAAGGAGAGCGTTAAACGACCAACAACACAAAAAAc                                                                        >  1:2253553/1‑150 (MQ=12)
                                                                            tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                       >  1:4260217‑M1/21‑150 (MQ=255)
                                                                            aCATGAGAATACCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                       <  1:3846884/150‑1 (MQ=255)
                                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACACATAAAACGa                                                                      >  2:1523453/1‑150 (MQ=255)
                                                                                tgatgcctgggtaattAATTCGTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                   >  2:3657642‑M1/17‑150 (MQ=255)
                                                                                 gatgcctggggatttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGAAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                  >  2:2953981‑M1/16‑150 (MQ=255)
                                                                                  gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCCCAACGTTCAAATCCCCTCCCGGCGGATTTGTCCTACGCAGGCGCGCGGTCTCCGCCAAACCCCAGATAAGAAGAAAGGc                                                                 >  1:2860867/1‑150 (MQ=14)
                                                                                  atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGc                                                                 <  2:3372534‑M1/136‑1 (MQ=255)
                                                                                      ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                             <  2:311023‑M1/140‑1 (MQ=255)
                                                                                       tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                       <  2:3506708‑M1/66‑1 (MQ=255)
                                                                                       tggggaattAATCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTAATCAGGAGAGGGTTAACCGACAAACAACAGATAAAAAGAAAGGCCCAGt                                                            <  1:2463217‑M1/141‑1 (MQ=255)
                                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                            >  1:1051097‑M1/9‑133 (MQ=255)
                                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                            <  2:1051097‑M1/125‑1 (MQ=255)
                                                                                         gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          >  1:1766945‑M1/8‑150 (MQ=255)
                                                                                         gATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          <  2:1105129/150‑1 (MQ=255)
                                                                                          aTCCCCATTTCGTGGCTGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                         <  1:2666175/150‑1 (MQ=11)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                         <  2:4136492/150‑1 (MQ=38)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                         <  2:2155885/150‑1 (MQ=38)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTCAAAACCCCTCCCGGCGGCGTTGTCCCACGCCGGCGAGCGGTCCCCGACAAACCACAGATAAAACCAAAGGCCCCGTCtt                                                         >  2:3679593/1‑150 (MQ=11)
                                                                                           accccAATTTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                         <  1:490327/149‑1 (MQ=11)
                                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGGTAACCGACAAACAACAGATAAAAACAAAAGCCCACTcttta                                                       <  2:2563305‑M1/146‑2 (MQ=255)
                                                                                             attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACACCAGCTAAAACGCAAGGCCCAGTCTTTCg                                                      <  2:4199050‑M1/147‑1 (MQ=255)
                                                                                               cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                    <  1:2882237/150‑1 (MQ=33)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                   <  2:1819434/150‑1 (MQ=32)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAAAGGTCAAAACCGCGCCCGGCGGGTTTGTCCTACACAGGAGGGCGTTCACCGCCCAACCACAGATAAAACGCAAGGGCCAGTCTTTCGTCt                                                   >  1:1999055/1‑150 (MQ=11)
                                                                                                 atttCCTGGCTTTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                   <  1:843478/148‑1 (MQ=9)
                                                                                                 aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                  >  2:2412019/1‑150 (MQ=32)
                                                                                                               tagggcgggcgCCCGCCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACCCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGtttt                                     <  1:2958077/147‑1 (MQ=1)
                                                                                                                             tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                       <  1:2418664/150‑1 (MQ=11)
                                                                                                                              gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTGGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGATGCCTgggg                      <  2:2931970/150‑1 (MQ=11)
                                                                                                                              cccccccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                      <  1:35423/149‑1 (MQ=14)
                                                                                                                                cccccACCCTCCGGGCCGGTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  1:2718110/147‑1 (MQ=12)
                                                                                                                                 cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   <  1:3012410/150‑1 (MQ=33)
                                                                                                                                 cGCCACCCGCCGGGCCGTTGCATCGCAACGTACAACTCCGCGCCCGGCGGATTGGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   <  2:1903363/150‑1 (MQ=11)
                                                                                                                                  gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGACTGGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                  <  2:2094991/150‑1 (MQ=11)
                                                                                                                                                   tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                 >  1:3932847/1‑134 (MQ=35)
                                                                                                                                                   tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                 >  1:3339285/1‑134 (MQ=35)
                                                                                                                                                     |                                                                                                                                                  
GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360008‑360302

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