Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,156 C→G 7.3% 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 7.3% 465.1 / 3.6 159intergenic (‑90/+505)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (49/96);  new base G (8/6);  total (57/102)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360009‑360302
                                                                                                                                                   |                                                                                                                                                  
gCACTTCACCGCTTGCCAGCGGCTTACCATACAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACAAGAGAATTCCGATCCCCAATTACTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                    >  1:2051355/1‑148 (MQ=255)
cctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                    <  1:3135229‑M1/53‑1 (MQ=255)
   actcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                                 <  1:2211333‑M1/56‑1 (MQ=255)
      caggagagcgttcaccgacaaacaacagagcaaccgaaccgcccaggctttcgactgagcctttcgatttatttgctgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                              >  2:2387012‑M1/90‑148 (MQ=255)
         cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           <  1:1189416/148‑1 (MQ=255)
              gCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCGGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  2:2798570/148‑1 (MQ=255)
              cgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACGTTCaa                                                                                                                                      <  2:1138522‑M1/67‑1 (MQ=255)
               ccAGCGGCTTACCAGCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                     <  2:3199522/148‑1 (MQ=255)
                cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                    <  1:2883837/148‑1 (MQ=255)
                 ggCGGCTTACCATCCCGCGCCCCCATCCTGTGCTGGAGCTCGGTTGATAAGGTGTCAAACATGCTAATTCCGAGCCCCATTTCGTGGCTGTTTTGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                   <  1:2408258/147‑1 (MQ=21)
                     cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                               <  2:1755815‑M1/74‑1 (MQ=255)
                      ctaaccaccCCGCGCCACCATCCCGTGCAGGAGCACCGTTGATAAGCTGAAAAACATGACAATCCCCACCCCCAATTCCTGGCAGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                              <  1:2137112/145‑1 (MQ=18)
                       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                             <  2:789045/148‑1 (MQ=255)
                        caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                            <  1:3250154‑M1/77‑1 (MQ=255)
                          ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGGCGCAACCTTGAAATCCGCTccccg                                                                                                                          >  1:2172956/1‑146 (MQ=37)
                          ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCCAGGTTCTAATGCGCTCCcgg                                                                                                                          >  1:2857485/1‑148 (MQ=37)
                          aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                          <  2:283330‑M1/79‑1 (MQ=255)
                            caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                        <  2:2548410‑M1/81‑1 (MQ=255)
                             tCCCGCGCCACCATCCCGTGCAGGAGCTCCTATGATACGCTGAAAAACATGATAATTCCGATCCCCATTTCCTGGATTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                       <  1:142202/148‑1 (MQ=18)
                             tCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGTCGCAACGATGAAATACCCCACcggcgg                                                                                                                       >  1:910621/1‑148 (MQ=25)
                              cccgataaaaccaacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      >  2:3054591‑M1/66‑148 (MQ=255)
                              ccagataaaacgaacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCTGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      >  2:3054601‑M1/66‑148 (MQ=255)
                              ccAGCGCCACCATCCAGGGCAGGAGCTCCGATTATACGCTGTAAAACATGAGAATTCCGATCCCCATTGCCTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  1:2823113/148‑1 (MQ=21)
                               cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                     <  2:1050899/148‑1 (MQ=255)
                                aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACCTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                    <  2:1442962/148‑1 (MQ=255)
                                 gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAACTTCAAATCCGCTCCCGGCGGAttt                                                                                                                   >  1:1466987/1‑148 (MQ=255)
                                 gataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                   >  2:205156‑M1/63‑148 (MQ=255)
                                 gataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGGTCAAATGCGCTCCCGGCGGAttt                                                                                                                   <  1:3141109‑M1/86‑1 (MQ=255)
                                   taaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                 <  2:2389568‑M1/88‑1 (MQ=255)
                                     caccaTCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCATTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                               <  1:1131884/148‑1 (MQ=255)
                                     caccaTCCAGTGCTGGAGCTCCGATGATACGCTGTCAAACATGCGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                               <  1:584299/148‑1 (MQ=255)
                                     caccaTCCAGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGAGATTTCCGATCCCCAATTCCGGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                               <  1:1131864/148‑1 (MQ=21)
                                      accaTCCAGTGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                              <  2:2643273/148‑1 (MQ=255)
                                        cgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTGGCGTCGCAACGTTAAAATCCGCTCCCGGCGGCTTTGTCCTAc                                                                                                            <  2:2423777‑M1/93‑1 (MQ=255)
                                           aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                         <  2:490288‑M1/96‑1 (MQ=255)
                                            aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAg                                                                                                        <  2:1723407‑M1/97‑1 (MQ=255)
                                             ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCAGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                       >  2:2919620‑M1/51‑148 (MQ=255)
                                             ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACCCACCCGGCGGATTTGTCCAACTCAgg                                                                                                       <  2:1272477‑M1/98‑1 (MQ=255)
                                             ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCTACTCAgg                                                                                                       <  2:1468343‑M1/98‑1 (MQ=255)
                                               cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATACGCTCCCGGCGGATTTGTCCTACTcaggtg                                                                                                     <  1:2973614‑M1/100‑3 (MQ=255)
                                                gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                    <  2:2637076/148‑1 (MQ=255)
                                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                    <  2:2577104‑M1/101‑1 (MQ=255)
                                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGGTGTCCTTCtcaggcga                                                                                                    <  2:2577119‑M1/101‑4 (MQ=255)
                                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGTCGCAAAGTTCCAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                    <  2:498297‑M1/101‑1 (MQ=255)
                                                 cagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAACCGCTCCCGGCCGATTTGTCCGACTCaggaggg                                                                                                   <  1:1964621‑M1/102‑3 (MQ=255)
                                                  agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                  <  2:471704‑M1/103‑1 (MQ=255)
                                                    tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                <  2:1780524‑M1/105‑1 (MQ=255)
                                                      gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATCCCCTCCCGGCGGATTTGTCCTACGAAGGAGAGCGGTc                                                                                              >  1:1951138/1‑148 (MQ=18)
                                                       ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCTACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGCGCGTTCa                                                                                             <  2:1923160‑M1/108‑1 (MQ=255)
                                                       ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:736003‑M1/108‑1 (MQ=255)
                                                       ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:1146776‑M1/108‑1 (MQ=255)
                                                       cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:2332783/148‑1 (MQ=255)
                                                       cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:2331659/148‑1 (MQ=255)
                                                         cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                          <  1:1470766‑M1/63‑1 (MQ=255)
                                                            ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                        <  2:1200560‑M1/113‑1 (MQ=255)
                                                                gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                  <  1:1936197‑M1/55‑1 (MQ=255)
                                                                gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                  >  2:1936197‑M1/32‑86 (MQ=255)
                                                                  ctttcgatttatttgatgcgtggggaattAATTCCTGGCTGTTTATGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                  >  1:751106‑M1/30‑148 (MQ=255)
                                                                    cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                <  2:1612993/148‑1 (MQ=255)
                                                                    cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                >  2:2231954/1‑148 (MQ=255)
                                                                      gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGa                                                                              >  2:3415242/1‑148 (MQ=255)
                                                                       tttttattttatggatggtgaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                             >  2:1903044‑M1/25‑148 (MQ=255)
                                                                       taaaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCTCAACGTTCAAATCCGCTCCCGGCGGTTTTGTCCTACTCTGGAGAGCGTTCACCTACAAACATCCGAt                                                                             <  2:2488477/146‑1 (MQ=17)
                                                                        cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTGGCAACGTTCAAAAACGCTCCCGGCGGATTTGTCCGACGAAGGAGAGCGGTCACCCACAAACAAAAGATa                                                                            >  1:3505901/1‑148 (MQ=17)
                                                                              tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                      <  1:2272671/148‑1 (MQ=255)
                                                                               gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACCTTAAAAAACGCTCCCCGCGGCTTTTTGCTAAACCGGAGAGCGTTGACCCCCAAAAAACAGATAAAAAGaa                                                                     >  1:2616663/1‑148 (MQ=14)
                                                                                   gcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                     <  2:2357270‑M1/100‑1 (MQ=255)
                                                                                   gcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                     >  1:2357270‑M1/13‑112 (MQ=255)
                                                                                   gcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                 >  2:1699351‑M1/13‑148 (MQ=255)
                                                                                   aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                 <  2:3493408/148‑1 (MQ=255)
                                                                                   aTTCCGATCACCATTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCTCAACGTTCAAATACGCTCCCGTCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACTGATAAAACGAAAGGc                                                                 <  2:2604923/148‑1 (MQ=11)
                                                                                    ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                <  2:2915811/148‑1 (MQ=255)
                                                                                    ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                >  2:2924845/1‑148 (MQ=255)
                                                                                    cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                <  2:2959541‑M1/57‑1 (MQ=255)
                                                                                    cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                >  1:2959541‑M1/12‑68 (MQ=255)
                                                                                      ccGTTCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                              <  2:1512569/148‑1 (MQ=255)
                                                                                      ccGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                              <  2:1966427/148‑1 (MQ=38)
                                                                                         ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                         >  1:3241364‑M1/7‑134 (MQ=255)
                                                                                         ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                         <  2:3241364‑M1/128‑1 (MQ=255)
                                                                                         ggaattAATTCCGGGCTGTTTCGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                           >  1:2577119‑M1/7‑148 (MQ=255)
                                                                                          tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAAAACGCTCCCCGCGGATTTGTGCTACTCAGGAGAGCGTTCACCGACAAAAAAAACAAAAAAAGAAAGGccccctat                                                          >  1:1147638/1‑143 (MQ=11)
                                                                                          gaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          >  1:877452‑M1/6‑148 (MQ=255)
                                                                                           ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTGGCAACCGTCAAAAACGCGCCCCGCGGCTTTGTGCTTCTCACGGGGGCGTTCACCGACAAAAAAAACATAAAAACAAAGGCCCACTcag                                                         >  1:1952131/1‑146 (MQ=2)
                                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCCGGCGAGCGTTCACCGCCAAACCACAGATAAAACGAAAGGCCCAGGCttt                                                        >  2:2844361/1‑148 (MQ=11)
                                                                                              taaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACGTCAAAAACCCTACCGGCGGGTTTGTGCTAATCCGGAGTTCGTTCAACGACAAACAAAAGATAAAAAGAAAGGCCCACTCattag                                                      >  1:3541459/2‑146 (MQ=2)
                                                                                              tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                      >  1:2291774‑M1/2‑148 (MQ=255)
                                                                                              aaaTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                      <  1:2291771/147‑1 (MQ=17)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                     >  2:695124/1‑148 (MQ=32)
                                                                                                  tCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                  <  1:1881703/148‑1 (MQ=11)
                                                                                                                             cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTATTTGTCCTACTCAGTAGAGCGTTCACCGACAAACAACAGATAAAACGAAATGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                       <  1:2456505/148‑1 (MQ=2)
                                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                     <  1:736003/148‑1 (MQ=32)
                                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATAAAAAGTAAGGCCCAGTCTTTCGACTGAGGCTTGCGTGTTCTTTGATGCCTgggcg                     >  2:2810562/1‑146 (MQ=11)
                                                                                                                                     ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCACCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAACGCCCAGTCTTTCGACTGAGCCATTCGTTTTATTTGGTGCCTGgggcataggg               >  2:1392102/1‑148 (MQ=11)
                                                                                                                                      cTCCGGGCAGTTTCTACGCAACGTTCAAATCCTCTCCCGGCGGATTTGTCCTACTCAGGAGATAGTTAACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGATTTATTTGATGCCTGgggaattaatt              <  1:3501759/148‑5 (MQ=11)
                                                                                                                                       tCCGGGCCGTTGCTTCGCAACCTTCAAAAACCCTCCCGGCGGATTTTTTCTACTCAGGAGAGCGGTCACCGCCAAACAACACATAAAAAGAAAGGCCCAGAAGTTCGACAGAGCCTTTCGGGTTATTTTTTTGCTGgggaattgggca             >  2:1527712/1‑146 (MQ=2)
                                                                                                                                           ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaa         <  2:3250154/148‑1 (MQ=255)
                                                                                                                                            gCCGTTGCTGCCCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat        <  1:1923160/148‑1 (MQ=255)
                                                                                                                                                ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACTTGGGAATTAATTCtt    >  1:1741292/1‑148 (MQ=255)
                                                                                                                                                ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt    >  1:1742373/1‑148 (MQ=255)
                                                                                                                                                  gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa  <  1:479403/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360009‑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: