Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,156 C→G 8.4% 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 8.4% 627.8 / 4.6 218intergenic (‑90/+505)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (54/142);  new base G (15/7);  total (69/149)
Fisher's exact test for biased strand distribution p-value = 3.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTT  >  CP009273/360008‑360300
                                                                                                                                                    |                                                                                                                                                
ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                  >  1:2967702/1‑149 (MQ=255)
   aCTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGACAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                               <  1:649383/149‑1 (MQ=255)
    cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTCg                                                                                                                                              >  1:2819568/1‑149 (MQ=255)
     ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                             <  2:2393237/149‑1 (MQ=255)
       cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCaa                                                                                                                                           >  1:845463/1‑149 (MQ=255)
         ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                         <  2:4050823‑M1/62‑1 (MQ=255)
          gagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                        <  2:2649733‑M1/63‑1 (MQ=255)
          cGCTCGCCCGCGGCAGACCACCCCGCCCCACCACCCAGTGCAGTAGAACCTATGATACGCTGACAAACAGGAGACTCCCCAACCCCATTGCCGGGCAGTTTATGGCGGGGGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                        <  1:1910360/149‑1 (MQ=12)
            gagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                <  1:3665874‑M1/55‑1 (MQ=255)
             ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCAACGTTCa                                                                                                                                     >  1:2597450/1‑149 (MQ=255)
              gcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCgccacggtcca                                                                                                                                    <  1:2989207‑M1/67‑6 (MQ=255)
               gCCCGCGGCTAACCACCCCGCCCCACCACCCAGTGCCGGAGCTCCGATTATACGGTGACAAACATGAGAATCCCGATCCCCAATCCCTGGATGTTTAGGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                   <  1:4013918/149‑1 (MQ=21)
                gttcaccgacaaacaacccctaacacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                  >  2:2508813‑M1/81‑149 (MQ=255)
                gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCGCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                  <  2:3650629‑M1/69‑1 (MQ=255)
                ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAAt                                                                                                                                  <  1:995733/149‑1 (MQ=255)
                ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGTTCAAAt                                                                                                                                  >  1:3962421/1‑149 (MQ=255)
                 cAGCGGCTTACCATCCAGAGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGTCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                 <  2:1246456/149‑1 (MQ=255)
                     ggCTTCCCATCCCGCGCCACCATCCAGTGCCGGAGCCCCTATTATAAGCTGTCAACCATTACAATTCCCACCCCAAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                             <  1:1640351/149‑1 (MQ=18)
                      gCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                            >  2:4021195/1‑149 (MQ=255)
                      cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                  <  1:2508807‑M1/69‑1 (MQ=255)
                       cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGTCGCAACGTTCAAATCCGCGcc                                                                                                                           >  1:3214827/1‑149 (MQ=255)
                           aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                       <  2:3348583‑M1/80‑1 (MQ=255)
                           aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                       <  2:3092572‑M1/80‑1 (MQ=255)
                              tCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                    <  2:715387/149‑1 (MQ=255)
                              tCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                    <  1:4053503/149‑1 (MQ=255)
                              tCCAGCGCCACCATCCAGGGCAGGAGCTCCGATGATAAGCTGTAAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                    <  2:715956/149‑1 (MQ=255)
                              accagatacacccaacggccccggctttcgacttcgcctttcgttttatttgaggcctggggaattAATGCCGGGCAGTTTATGGCGGGGGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                    >  1:792459‑M1/67‑149 (MQ=255)
                              aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                    >  1:2649733‑M1/67‑149 (MQ=255)
                              aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCccggcggc                                                                                                                    <  2:409014‑M1/83‑2 (MQ=255)
                               ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                   <  2:868052/149‑1 (MQ=255)
                                cagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCGACGTGCTAATCcgctcccgggggcgg                                                                                                                  <  2:2384999‑M1/85‑7 (MQ=255)
                                cagataaaaccacaggcccagacgttcgcctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                  >  2:1981076‑M1/65‑149 (MQ=255)
                                    taaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                              >  2:2961107‑M1/61‑149 (MQ=255)
                                     ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                             <  1:1715119/149‑1 (MQ=255)
                                      aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCTCCCGGCGGCTTTGTCCt                                                                                                            <  2:2842798‑M1/91‑1 (MQ=255)
                                      aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                            <  2:818982‑M1/91‑1 (MQ=255)
                                       nccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGTCGCAACCGTCAAATCCGCACCCCGCGGGTTTGtgcga                                                                                                           >  1:880983/2‑145 (MQ=25)
                                       accaTCCAGTGCAGGAGCTCCGATGATACCCTGACAAACATGAGAATTCCGATCCCCAATTCCGGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                           <  1:2524631/149‑1 (MQ=21)
                                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGGTGGCAACCTTCAAATACGCTCCCGGCGGGGTTGtgcga                                                                                                           >  1:1760565/1‑145 (MQ=25)
                                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                             <  1:2720868‑M1/74‑1 (MQ=255)
                                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                             >  2:2720868‑M1/58‑131 (MQ=255)
                                         caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGTTCAAATAACCTACCGGCGGAGTTGTCCTACt                                                                                                         >  1:3722955/1‑149 (MQ=25)
                                            ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                      <  2:2051220/149‑1 (MQ=255)
                                             aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTTtcctacgccgg                                                                                                     <  1:2313252‑M1/98‑6 (MQ=255)
                                              ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                    <  2:1871511‑M1/99‑1 (MQ=255)
                                               gTGCAGTAGCTCCGATTATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCATGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                   <  2:1431020/149‑1 (MQ=255)
                                                 gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                 <  2:504750/149‑1 (MQ=255)
                                                  cagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                >  2:155022‑M1/47‑149 (MQ=255)
                                                  cagtatttcgactgagactgtctaaatatttgttgcctcgcgaaatAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                >  2:845463‑M1/47‑149 (MQ=255)
                                                    gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAACTTCAAATACGCTCCCCGCGGATTTGTGCTACTACGGAGAGCGt                                                                                              <  1:3574884‑M1/105‑1 (MQ=255)
                                                    gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGCCTCAGGAGAGCGt                                                                                              <  2:2106142‑M1/105‑1 (MQ=255)
                                                        cTCCTATGATAAGCTGACAAACATGAGAATTCCGAGCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                          <  1:3895618/149‑1 (MQ=255)
                                                         tCCGATGATACGCTGACAAACATGAGAATTCCGATCCCCAATCCCTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                         <  1:3503291/149‑1 (MQ=17)
                                                         tCCGAGGAGAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                         <  2:3579846/149‑1 (MQ=255)
                                                          ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                        <  2:2819568/149‑1 (MQ=255)
                                                          ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                        <  1:1374216/149‑1 (MQ=255)
                                                            gatgatAAGGTGTCAAACATTTGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                      <  2:76857/149‑1 (MQ=255)
                                                              tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAAACCGCACCCGGCGGCTTTTTCCAAAAAAGGAGAGCGTTGACCGacaa                                                                                    >  1:3247241/1‑149 (MQ=18)
                                                               gagcctttctatttatttgatgccggggcaattACTTCCGGGCATTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                   >  1:1538096‑M1/34‑149 (MQ=255)
                                                                atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                  >  2:1613998/1‑149 (MQ=255)
                                                                atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                  >  2:1181648/1‑149 (MQ=255)
                                                                 gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                 <  2:2864715‑M1/118‑1 (MQ=255)
                                                                  aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGAATCAGGAGAGCGTTCACCGacaaacaa                                                                                >  2:1699209/1‑149 (MQ=255)
                                                                     ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                             <  2:3633572‑M1/122‑1 (MQ=255)
                                                                     ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCACCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAAAAACAg                                                                             <  2:3470465‑M1/122‑1 (MQ=255)
                                                                     cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                             <  2:3495251/149‑1 (MQ=255)
                                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                            >  2:1798434/1‑149 (MQ=255)
                                                                       gTCAAACATGAGACTACCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                           <  1:549096/149‑1 (MQ=255)
                                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACGCCGGAGAGCGTTCACCCCCAAACAACAGAt                                                                           >  2:2477431/1‑149 (MQ=17)
                                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGTCGCAAACTTCAAAAACGCCCCCGGCGGATTTGTCCAACTCAGGAGAGCGTTGACCCACCAAAAACAGAt                                                                           >  1:1720198/1‑149 (MQ=17)
                                                                       cgttttatttgatgcctggggaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCATATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                           >  2:3869496‑M1/26‑149 (MQ=255)
                                                                        gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                              <  2:2754511‑M1/105‑1 (MQ=255)
                                                                        gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                              >  1:2755104‑M1/25‑129 (MQ=255)
                                                                        gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                          <  2:2017875‑M1/125‑1 (MQ=255)
                                                                         cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAAAGGTAAAATACGCGCCCCGCGGCTTTTTCGTACTAAGGGGAGCGGTCACCGACAAAAAACACATaa                                                                         >  1:1981076/1‑149 (MQ=17)
                                                                         cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACCCAACCCGCGGCTTTTTCCGACGAAGGAGAGCGTTCACCGACCAacaacacaaaa                                                                         >  1:1979909/1‑144 (MQ=17)
                                                                          tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAACCGTCAAATACGCTCCCCGCGGATTTGTACTACTCAGGAGAGCGTTCACcgacaaacaacaaaaaaa                                                                        <  1:1297100‑M1/127‑7 (MQ=255)
                                                                          tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                        <  2:629261‑M1/127‑1 (MQ=255)
                                                                            tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                      <  2:3900192‑M1/129‑1 (MQ=255)
                                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCCACCGTCAAAAACCCGCCCCGCGGGGTTTTCCGACACAGGCGAGCGTTCACCGCCAAACAACATATaaaaac                                                                     >  1:2930847/1‑147 (MQ=12)
                                                                                 gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAATCAACAGATAAAACGCAAg                                                                 <  2:3986617‑M1/134‑1 (MQ=255)
                                                                                 agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                 <  2:2563831/149‑1 (MQ=255)
                                                                                   tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                               <  2:3227989‑M1/136‑1 (MQ=255)
                                                                                     ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                             <  2:2163346/149‑1 (MQ=255)
                                                                                       ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCACCGGGCCGTTGCTTCGCAGAGTTCAAATCCGCTCCCGGCGGAATTGTCCTTCGCAGGAGAGGGTTAAACGACCAACCACAGATGAATCGGAAGGccccg                                                           >  2:3513493/1‑147 (MQ=14)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                        <  2:3445986/149‑1 (MQ=38)
                                                                                           gaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAAAACCCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGGTGACCGACACACAACACATAAAACGAAAGGCCCAGGCtt                                                       <  1:2471221‑M1/144‑1 (MQ=255)
                                                                                            ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                      <  2:3806290/149‑1 (MQ=35)
                                                                                            ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                      >  2:1361268/1‑149 (MQ=35)
                                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACCTTCAAATACGCTCCCCGCGGATTTGTCCTACTCAGGGGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTcgg                                                   <  2:1053608‑M1/148‑2 (MQ=255)
                                                                                                                         gccggCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGcc                         <  1:1159018/145‑1 (MQ=1)
                                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                    <  1:1949514/149‑1 (MQ=31)
                                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                    <  2:2337510/149‑1 (MQ=14)
                                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                   <  2:1760565/149‑1 (MQ=17)
                                                                                                                                cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                  <  1:46786/149‑1 (MQ=32)
                                                                                                                                  ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                <  1:2229248/149‑1 (MQ=34)
                                                                                                                                     cccTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACCAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg             <  1:1045697/149‑1 (MQ=255)
                                                                                                                                             gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt     <  2:2987805/149‑1 (MQ=255)
                                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt  <  1:3092572/149‑1 (MQ=255)
                                                                                                                                                    cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                  >  2:1943641/1‑129 (MQ=33)
                                                                                                                                                    |                                                                                                                                                
GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTT  >  CP009273/360008‑360300

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: