Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 10.1% intergenic (‑107/+488) 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,1730CA10.1% 388.4 / 2.4 158intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (56/74);  new base A (14/10);  total (72/86)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

GCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTG  >  CP009273/360027‑360320
                                                                                                                                                  |                                                                                                                                                   
gCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATCAGCTGTCAAACATGAGAATTCCGATCCCCCATTCCTGGCAGGTTATGGCGGGCGTCCTGCCCGCCACCCTCCGTGCCGTGGCGTCGCAACGTTCAAATcc                                                                                                                                                    >  1:1764447/1‑148 (MQ=37)
   cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                 <  2:2570966‑M1/74‑1 (MQ=255)
    gacaaccaccagatcaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                >  2:2664637‑M1/74‑148 (MQ=255)
    cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGGAGCAACGTTAAAATCCGcac                                                                                                                                                >  1:2579581/1‑146 (MQ=37)
    cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                >  2:1340227/1‑148 (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)
           tCCACCCCCACCATCCAGTGCAGGAGCGCCGATTATAAGCTGAAAAACATGATAATTCCGATCCCCACTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                         <  1:859646/148‑1 (MQ=21)
            acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                        <  1:69619‑M1/83‑1 (MQ=255)
             cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGGGGCCACGTTGAAATACGCGGCcggcggga                                                                                                                                       >  1:1241140/1‑146 (MQ=37)
              aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGAtt                                                                                                                                      >  1:2967113/1‑148 (MQ=255)
               gataacaccacaggcccagtcgttcgactgagcctgtcgatttatttgaggcctggggaagtAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                     >  2:1840840‑M1/63‑148 (MQ=255)
               gataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                     <  2:248567‑M1/86‑1 (MQ=255)
                ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                    <  2:2595035‑M1/87‑1 (MQ=255)
                  ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                  <  2:531767/148‑1 (MQ=255)
                   caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                 <  2:578001/148‑1 (MQ=255)
                     ccaTCCCGGGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATGCCGATCCCCATGTCCGGGCAGTTTAGGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                               <  1:3393060/148‑1 (MQ=18)
                       aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                             <  2:3132497/148‑1 (MQ=255)
                       aTCCAGTGCAGGAGCGCCGATGAAAAGCTGACAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                             <  2:100744/148‑1 (MQ=255)
                        tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACtg                                                                                                                            <  2:1032750/148‑2 (MQ=255)
                             tGCAGGAGCTCCGATGATAAGCTGTAAACCATGAGAATTCCGAGCCCCAATTCCGGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                       <  1:3141219/148‑1 (MQ=18)
                             cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                       >  1:250778‑M1/49‑148 (MQ=255)
                               cAGGCGCTCCTATGATACGCTGACAAACATGACAATGCCCAGCCCAAATTCCTGGCTGTTTAGGGCGGGCGTCCGCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                     <  1:2103819/148‑1 (MQ=17)
                                agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATGTGTGCTACTCCGGCGAGc                                                                                                                    <  1:2024324‑M1/103‑1 (MQ=255)
                                 ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                   >  1:3366001/1‑148 (MQ=255)
                                  gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                  <  1:564782/148‑1 (MQ=255)
                                  gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                  <  2:1342610/148‑1 (MQ=255)
                                  gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCCACGTTGAAATGAACGCCACGCGGAGTTGTGCGTCTCAGGAGTGCGt                                                                                                                  >  1:964206/1‑148 (MQ=21)
                                   aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCAGCACCCGGCGGATTTGTCCTACTCAGGCGAGCGtt                                                                                                                 >  2:375750/1‑148 (MQ=255)
                                    gCTCCGATGAGACGCTGTCAAACATCACAATGCCGACCCCCAATGCCGGGCAGTTTCGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                <  1:400536/148‑1 (MQ=14)
                                       ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATCCCCTCCCCGCGGATTTGTCCTAATCCGGAGAGCGTTAAcc                                                                                                             >  2:910435/1‑148 (MQ=18)
                                         tatcttAAGCTGTCAAAATTTAGTCTTCCGCTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                                           <  2:3037884/143‑1 (MQ=14)
                                          ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACCCTCCCGGCGGCTGTGTCCTACTCAGGAGAGCGGTCACCGac                                                                                                          <  2:2685424‑M1/113‑1 (MQ=255)
                                           tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                                         <  2:1694893‑M1/114‑1 (MQ=255)
                                            gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGTGGAAATCCGCGCCCCGGGGTTTTGTGATACTAAGGAGACCGTTGACCGacca                                                                                                        >  1:2594975/1‑146 (MQ=18)
                                            gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACCGTCAAAAACCCTCCCCGCCGATTTGTACAACTCAGGAAAGCGGTAACCCacaa                                                                                                        <  2:3273557‑M1/115‑1 (MQ=255)
                                             atAAGCTGTCAAACGTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                       <  2:888125/148‑1 (MQ=255)
                                             atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGGTCAAAAACGCACCCCCCGGCATAGTGCTACAAAGGTGAGCGTTCACCGacaaa                                                                                                       >  1:2654892/1‑148 (MQ=18)
                                             atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACGCACCCGGCGGAAGTGTGCTACTGAGGAGAGCGTTCACCGaccaa                                                                                                       >  1:2654307/1‑148 (MQ=18)
                                                aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACCTTCAAATCCGCGCCCCGCGGCTTTGTGCTACGCAGGAGAGCGGTGACCCacaaacaa                                                                                                    >  1:3037884/1‑148 (MQ=17)
                                                 tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCACCCCGCGGAATTTTACAAAGCAGGAGAGCGGTGAACGACAaacaac                                                                                                   <  1:3512469‑M1/120‑1 (MQ=255)
                                                 gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAACCTTAAAATACGCTCCCGGCGGATTTGTGCTACTAAGGAGAGCGGTCACCGACAaaaaac                                                                                                   >  1:543110/1‑148 (MQ=17)
                                                 gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                                   <  2:346117/148‑1 (MQ=255)
                                                  cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTGAGGAGAGCGTTCACCGACAaaccacc                                                                                                  >  2:3207991/1‑147 (MQ=255)
                                                      ttttatttgttggctggcgtattAATTGCTGGCAGTTTATGTCGGGCGTCCTGCGCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCTCCGTCAAACAACAGATa                                                                                              >  2:2654892‑M1/24‑148 (MQ=255)
                                                       tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                             >  1:2595035‑M1/23‑148 (MQ=255)
                                                        ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                            >  2:69619‑M1/22‑148 (MQ=255)
                                                        aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTGGTGGCCACCTTCAAAAAAGCACCCCGCGGCTTTGTCCGACGCAGGCGTGGGGTTACCGCACAACAACACATaaa                                                                                            >  1:1272402/1‑148 (MQ=12)
                                                        aaCATGAGAATCCCGACCCCCAATTCCTGGATGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                            <  1:127176/148‑1 (MQ=14)
                                                           tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCTACCGGCCGAATTTTGCGACAAAGGAGAGCGTTGACCCACAAACAACagaaaaaaag                                                                                         <  1:3236504‑M1/130‑3 (MQ=255)
                                                           tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                          >  1:1490111‑M1/19‑147 (MQ=255)
                                                            tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCGACCGTGAAAAACGCGCCCCGCCGCTGTTTCCGACTCAGGGGAGCGGTGAACGACACACAACACAAAAAAcca                                                                                        >  1:3383868/1‑146 (MQ=9)
                                                             gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                       <  2:1355560/148‑1 (MQ=255)
                                                             gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                       <  2:2759330/148‑1 (MQ=255)
                                                             gagaAGTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                       <  2:2374299/148‑1 (MQ=255)
                                                               atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                     <  2:120303‑M1/134‑1 (MQ=255)
                                                                tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATAAAACGAAAgg                                                                                    <  2:2709742‑M1/135‑1 (MQ=255)
                                                                  ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATACGCTCCCGGCGGGTTTGTGCGACTCAGGCGTGCGTTGACCCACCAACAACAGATAAAAAGAAATGcc                                                                                  >  1:1746952/1‑148 (MQ=14)
                                                                      gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                              <  2:1728108‑M1/141‑1 (MQ=255)
                                                                      gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                              <  2:2810442/148‑1 (MQ=38)
                                                                      gATCCCCAATTCCTGGCAGTTTATGGCGGGCCTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCAGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                              <  2:645144/148‑1 (MQ=255)
                                                                       ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCACCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAAAACATAAAACGAAAGGcccactc                                                                             <  2:2394705‑M1/142‑4 (MQ=255)
                                                                       aTCCCCATTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                             <  1:574876/148‑1 (MQ=25)
                                                                       aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                             <  2:1850625/148‑1 (MQ=38)
                                                                           ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                            <  1:2323430/113‑1 (MQ=35)
                                                                             aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                       <  2:370804/148‑1 (MQ=32)
                                                                             aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                       <  2:2779504/148‑1 (MQ=32)
                                                                               ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                     >  2:2717495/1‑148 (MQ=31)
                                                                                tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                    <  2:3366001/148‑1 (MQ=11)
                                                                                tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACCAACAACAGATAAAACCAAAGGCCCCGTCTTTCGACTg                                                                    >  2:2280371/1‑148 (MQ=2)
                                                                                                           cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                         <  1:196716/148‑1 (MQ=11)
                                                                                                            ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGTGCCTTTCGTTTTATTTGTTGCCTgggg                                        >  2:583525/1‑148 (MQ=9)
                                                                                                            ccGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTGCGTTTTATTTGATGCCTgggg                                        >  2:583490/1‑148 (MQ=9)
                                                                                                              gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                      <  2:1496404/148‑1 (MQ=17)
                                                                                                              gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                      <  1:1495498/148‑1 (MQ=17)
                                                                                                               ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                     <  1:2280371/148‑1 (MQ=33)
                                                                                                                 aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAAAAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGGTTTATTTGTTGCCTGgggcattg                                   >  2:3090827/1‑148 (MQ=11)
                                                                                                                    cTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGCCTGAGCCGTTCGTTTTATTTGATGGCGGgggaattggga                                >  2:1111089/1‑148 (MQ=11)
                                                                                                                      ccGGGCCGTTGCTTCGCAACGTTCATATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat                              <  1:2837453/148‑1 (MQ=255)
                                                                                                                          gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat                          >  2:270156/1‑148 (MQ=255)
                                                                                                                            cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACACCAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc                        <  1:3544542/148‑1 (MQ=255)
                                                                                                                                gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                  <  1:250416/134‑1 (MQ=37)
                                                                                                                                gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                  >  2:250416/1‑134 (MQ=37)
                                                                                                                                   tCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                 <  1:875616/148‑1 (MQ=255)
                                                                                                                                       aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCGTGGGGAATTGGGAATTAATTCTTGAAGACGaa             >  2:2840428/1‑148 (MQ=255)
                                                                                                                                       aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaa             >  2:2840733/1‑148 (MQ=255)
                                                                                                                                        aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGaaa            >  1:2372335/1‑148 (MQ=255)
                                                                                                                                         cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                 <  1:2803338/142‑1 (MQ=255)
                                                                                                                                         cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                 >  2:2803338/1‑142 (MQ=255)
                                                                                                                                            tCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc        >  1:529741/1‑148 (MQ=255)
                                                                                                                                            tCAAATCCGCTCCCGGCGGAATTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc        <  1:1075230/148‑1 (MQ=255)
                                                                                                                                                aTCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg    <  2:3028715/148‑1 (MQ=255)
                                                                                                                                                aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg    >  2:147772/1‑148 (MQ=255)
                                                                                                                                                aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg    <  2:1257888/148‑1 (MQ=255)
                                                                                                                                                  ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg  >  1:171785/1‑148 (MQ=255)
                                                                                                                                                  ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg  >  2:901029/1‑148 (MQ=255)
                                                                                                                                                  ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAAAAAACGAAAGGCCCACTCTTTCGACTGAGACTATCTTTTTATATGATGAATGGGGAATAGGGAATTAATTATTTAATAATAAAGGGCCTCGTg  >  2:2219749/1‑148 (MQ=14)
                                                                                                                                                  |                                                                                                                                                   
GCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTG  >  CP009273/360027‑360320

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: