Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 15.6% intergenic (‑91/+504) 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,1570TG15.6% 311.6 / 14.0 174intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (51/91);  new base G (22/10);  total (73/101)
Fisher's exact test for biased strand distribution p-value = 1.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

ACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360016‑360302
                                                                                                                                             |                                                                                                                                                 
aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCaa                                                                                                                                             <  2:2523501‑M1/60‑1 (MQ=255)
 ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCCAc                                                                                                                                            >  1:104682/1‑148 (MQ=255)
  cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACg                                                                                                                                           >  1:56290/1‑148 (MQ=38)
    cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         <  2:2151021/148‑1 (MQ=255)
    cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         >  2:246188/1‑148 (MQ=255)
      tGCCAGCGGCTTACCATCCCGCGCCACCATCCAGGGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                       <  1:1185586/148‑1 (MQ=255)
       gcccggggCTTCCCACCCCGCGCCACCATCACGTGCTGGAGCTCCGTTTATACGCTGTAAAACATGAGATTTCCGAGCCCCATTGCGGGGCTTTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  1:1702224/142‑1 (MQ=12)
       gCCAGCGGCTAACCCACCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTAAAACATGATATTTCCGATCCCCATTTCCTGGCTTTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  1:1978609/148‑1 (MQ=25)
       cgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  2:2352011‑M1/67‑1 (MQ=255)
        ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                     <  1:3555071/148‑1 (MQ=255)
        ccAGCGGCTTACCAACCAGCGCCCCCATCCAGTGCAGGAGCTCCGATGATAAGCTTTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                     <  1:3555036/148‑1 (MQ=255)
          tcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAata                                                                                                                                   <  1:827398‑M1/70‑2 (MQ=255)
              cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                               <  2:2570966‑M1/74‑1 (MQ=255)
               cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                              >  2:1340283/1‑148 (MQ=255)
               cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGcgc                                                                                                                              >  1:2579595/1‑146 (MQ=255)
                acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           >  2:2661278‑M1/73‑134 (MQ=255)
                acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           >  2:2661281‑M1/73‑134 (MQ=255)
                acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           <  1:2661281‑M1/62‑1 (MQ=255)
                acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACg                                                                                                                                           <  1:2661278‑M1/62‑1 (MQ=255)
                acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATccgcgcc                                                                                                                             <  1:2867577‑M1/76‑4 (MQ=255)
                 tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGGCGTGGGTTCGCACCGGGCAAATCCGCTccc                                                                                                                            >  1:1056232/1‑148 (MQ=37)
                  aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattTATTCCTGGCAGTTTATGGCGGGCGTGCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTGAAATGGGCTCCcg                                                                                                                           <  2:945614‑M1/78‑1 (MQ=255)
                  aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                           <  2:1495498‑M1/78‑1 (MQ=255)
                     caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAACGTTCAAATCCGCTCCcggcg                                                                                                                        <  1:2016823‑M1/81‑1 (MQ=255)
                     caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                        >  1:2352011‑M1/68‑148 (MQ=255)
                      tCCACCCCCACCATCCAGTGCAGGAGCGCCGATTATAAGCTGAAAAACATGATAATTCCGATCCCCACTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                       <  1:859646/148‑1 (MQ=21)
                       ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  2:3213774/148‑1 (MQ=255)
                       acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  1:69619‑M1/83‑1 (MQ=255)
                       acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGTGGCAACGTTGAAATGCGcgcccgggggc                                                                                                                      <  1:1922557‑M1/83‑5 (MQ=255)
                        cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGGGGCCACGTTGAAATACGCGGCcggcggga                                                                                                                     >  1:1241140/1‑146 (MQ=37)
                         aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGAtt                                                                                                                    >  1:2967113/1‑148 (MQ=255)
                          gataacaccacaggcccagtcgttcgactgagcctgtcgatttatttgaggcctggggaagtAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                   >  2:1840840‑M1/63‑148 (MQ=255)
                           ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                  <  2:2595035‑M1/87‑1 (MQ=255)
                             ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                <  2:531767/148‑1 (MQ=255)
                              caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                               <  2:578001/148‑1 (MQ=255)
                               accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGTCAACGATCAAAAACCCACCCGGCGGCtttttgcc                                                                                                              >  1:2598268/1‑143 (MQ=25)
                                ccaTCCCGGGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATGCCGATCCCCATGTCCGGGCAGTTTAGGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                             <  1:3393060/148‑1 (MQ=18)
                                  aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                           <  2:3132497/148‑1 (MQ=255)
                                  aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                           <  2:100725/148‑1 (MQ=255)
                                   tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACtg                                                                                                          <  2:2804614/148‑2 (MQ=255)
                                   tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACtg                                                                                                          <  2:1032750/148‑2 (MQ=255)
                                        tGCAGGAGCTCCGATGATAAGCTGTAAACCATGAGAATTCCGAGCCCCAATTCCGGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                     <  1:3141219/148‑1 (MQ=18)
                                        cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTcaggtg                                                                                                     <  2:419431‑M1/100‑3 (MQ=255)
                                         gCAGGAGCTCCGATGATAAGCTGTCAAACATGACAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                    <  1:1719856/148‑1 (MQ=255)
                                          cAGGCGCTCCTATGATACGCTGACAAACATGACAATGCCCAGCCCAAATTCCTGGCTGTTTAGGGCGGGCGTCCGCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                   <  1:2103819/148‑1 (MQ=17)
                                           agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATGTGTGCTACTCCGGCGAGc                                                                                                  <  1:2024324‑M1/103‑1 (MQ=255)
                                            ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                 >  1:3366001/1‑148 (MQ=255)
                                               gCTCCGATGAGACGCTGTCAAACATCACAATGCCGACCCCCAATGCCGGGCAGTTTCGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                              <  1:400536/148‑1 (MQ=14)
                                                cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:172699/148‑1 (MQ=255)
                                                  ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATCCCCTCCCCGCGGATTTGTCCTAATCCGGAGAGCGTTAAcc                                                                                           >  2:910435/1‑148 (MQ=18)
                                                   cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                          <  2:763252/148‑1 (MQ=255)
                                                     ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACCTTCAAAAACGCTCCCGGCGGATTTGTCCTAATCAGGAGAGCGTTCACCGac                                                                                        <  1:2576470‑M1/113‑1 (MQ=255)
                                                      tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                       <  2:1694893‑M1/114‑1 (MQ=255)
                                                       gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGTGGAAATCCGCGCCCCGGGGTTTTGTGATACTAAGGAGACCGTTGACCGacca                                                                                      >  1:2594975/1‑146 (MQ=18)
                                                        atAAGCTGTCAAACGTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                     <  2:888125/148‑1 (MQ=255)
                                                        atAAGCTGTCAAACGTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                     <  2:2669911/148‑1 (MQ=255)
                                                        atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGGTCAAAAACGCACCCCCCGGCATAGTGCTACAAAGGTGAGCGTTCACCGacaaa                                                                                     >  1:2654892/1‑148 (MQ=18)
                                                            tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCACCCCGCGGAATTTTACAAAGCAGGAGAGCGGTGAACGACAaacaac                                                                                 <  1:3512469‑M1/120‑1 (MQ=255)
                                                            gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAACCTTAAAATACGCTCCCGGCGGATTTGTGCTACTAAGGAGAGCGGTCACCGACAaaaaac                                                                                 >  1:543110/1‑148 (MQ=17)
                                                            gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                 <  2:346117/148‑1 (MQ=255)
                                                             cTGTAAAACATGAGATTTCCGAGCCCCATTTCCTGGATGTTTAGGGCGGGCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                <  1:3027736/148‑1 (MQ=14)
                                                              tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                               <  2:3143520/148‑1 (MQ=255)
                                                                 ttttatttgttggctggcgtattAATTGCTGGCAGTTTATGTCGGGCGTCCTGCGCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCTCCGTCAAACAACAGATa                                                                            >  2:2654892‑M1/24‑148 (MQ=255)
                                                                  aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCAAAGTGCAAATCCGCGCCCGGCGGCTTGGTGCTACTCCGGGGCTGGTTGACCGCCCAAAAACAGAtta                                                                           >  1:2466093/1‑146 (MQ=17)
                                                                   aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTGGTGGCCACCTTCAAAAAAGCACCCCGCGGCTTTGTCCGACGCAGGCGTGGGGTTACCGCACAACAACACATaaa                                                                          >  1:1272402/1‑148 (MQ=12)
                                                                   aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGGAGCCACCTGAAAAAACCCGCCCCGCGGCTTTGGGCAACACAGGAGAGCGGGCACCCCCCAACAACACATaaa                                                                          >  1:1272338/1‑148 (MQ=12)
                                                                      tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCTACCGGCCGAATTTTGCGACAAAGGAGAGCGTTGACCCACAAACAACagaaaaaaag                                                                       <  1:3236504‑M1/130‑3 (MQ=255)
                                                                      aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAAACCGCTCCCCGCGGCGTTGTGCTACTCCGGAGAGCGTTCGCCGACACACAACAGATAAAACg                                                                       >  2:2688240/1‑148 (MQ=14)
                                                                       ttgtaggctggtgaatgAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAAGAGATAAAACGa                                                                      >  2:2967113‑M1/18‑148 (MQ=255)
                                                                       tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCGACCGTGAAAAACGCGCCCCGCCGCTGTTTCCGACTCAGGGGAGCGGTGAACGACACACAACACAAAAAAcca                                                                      >  1:3383868/1‑146 (MQ=9)
                                                                       tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCGTAAAAAACCCGCCCGGCGGCGTTGTGCTACAAAGGGGGGGGGTAACCCACAAAAAACCAATAAAAAGa                                                                      >  1:2827843/1‑148 (MQ=9)
                                                                       tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATACTACGa                                                                      >  2:2731927/1‑148 (MQ=255)
                                                                        gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                     <  2:1355560/148‑1 (MQ=255)
                                                                           tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATAAAACGAAAgg                                                                  <  2:2709742‑M1/135‑1 (MQ=255)
                                                                           tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGTAAgg                                                                  <  2:2709827‑M1/135‑1 (MQ=255)
                                                                            aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCGTCAAATACGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGGTCACCGACAAACAACACATAAAAACAAAgtc                                                                 >  1:550270/1‑146 (MQ=14)
                                                                             ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATACGCTCCCGGCGGGTTTGTGCGACTCAGGCGTGCGTTGACCCACCAACAACAGATAAAAAGAAATGcc                                                                >  1:1746952/1‑148 (MQ=14)
                                                                                 gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            >  2:160661‑M1/8‑148 (MQ=255)
                                                                                 gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            <  2:2810442/148‑1 (MQ=38)
                                                                                 gATCCCCAATTCCTGGCAGTTTATGGCGGGCCTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCAGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            <  2:645144/148‑1 (MQ=255)
                                                                                  ggaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGGCGCAACGGTCCAATCCGGTCCCGGCGGATTTGGGCTACTCCGGCGAGCGTTCCCCGACCCACCCCAGCTAAATCGAACGGCCCCGtg                                                           >  2:2396201/7‑147 (MQ=2)
                                                                                  ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCACCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAAAACATAAAACGAAAGGcccactc                                                           <  2:2394705‑M1/142‑4 (MQ=255)
                                                                                  aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                           <  2:1850625/148‑1 (MQ=38)
                                                                                     cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                        <  1:569562/148‑1 (MQ=35)
                                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                     <  2:2779504/148‑1 (MQ=32)
                                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                     <  2:370804/148‑1 (MQ=32)
                                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                     >  2:1154023/1‑148 (MQ=17)
                                                                                           tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                  <  2:3366001/148‑1 (MQ=11)
                                                                                           tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACCAACAACAGATAAAACCAAAGGCCCCGTCTTTCGACTg                                                  >  2:2280371/1‑148 (MQ=2)
                                                                                                      tagggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGtt                                       <  2:2139919/145‑1 (MQ=1)
                                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGTGCCTTTCGTTTTATTTGTTGCCTgggg                      >  2:583525/1‑148 (MQ=9)
                                                                                                                       ccGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTGCGTTTTATTTGATGCCTgggg                      >  2:583490/1‑148 (MQ=9)
                                                                                                                         gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  1:1495498/148‑1 (MQ=17)
                                                                                                                          ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   <  1:2280371/148‑1 (MQ=33)
                                                                                                                           cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                  >  2:3079892/1‑148 (MQ=34)
                                                                                                                               cTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGCCTGAGCCGTTCGTTTTATTTGATGGCGGgggaattggga              >  2:1111089/1‑148 (MQ=11)
                                                                                                                                     gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTTATGCCTGGGGAATTGGGaattaat        >  2:1375124/1‑148 (MQ=255)
                                                                                                                                        gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt     >  2:1304029/1‑148 (MQ=255)
                                                                                                                                          tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTg   >  2:2499496/1‑148 (MQ=255)
                                                                                                                                           gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                <  1:250416/134‑1 (MQ=37)
                                                                                                                                             |                                                                                                                                                 
ACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360016‑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: