Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,181 G→C 7.7% intergenic (‑115/+480) 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,1810GC 7.7% 595.7 / 6.2 204intergenic (‑115/+480)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (52/134);  new base C (10/8);  total (62/142)
Fisher's exact test for biased strand distribution p-value = 2.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACG  >  CP009273/360033‑360325
                                                                                                                                                     |                                                                                                                                                
tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTGCAAATCCGCGCCcg                                                                                                                                                  >  2:3533708/1‑149 (MQ=37)
  aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattTATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCTACGTTCAAATCCGCTCCcggc                                                                                                                                                <  2:2696411‑M1/80‑1 (MQ=255)
    caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCATGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCGcggcgg                                                                                                                                              <  2:274263‑M1/82‑1 (MQ=255)
    caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttagttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                              <  1:4351098‑M1/82‑1 (MQ=255)
     tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                             <  2:3032783/149‑1 (MQ=255)
     tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCACCcggcggc                                                                                                                                             >  1:4381440/1‑148 (MQ=255)
      ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                            <  2:767952/149‑1 (MQ=255)
      acagataaaacgaaaggcccagtctttcgactgagcctttcgtgtgatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGGTCCAATCCGCGcccggcggct                                                                                                                                            <  2:4412122‑M1/84‑3 (MQ=255)
       cagataaaacgaccggccccgtctttcgactgagcctattgttttttttgatgcgtgggtaattAATTCCTGTATGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                           >  1:67071‑M1/65‑149 (MQ=255)
        aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                          >  1:2538986/1‑149 (MQ=255)
           gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACCTGAGAATTCCGATCCCCCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                       <  2:2517679/149‑1 (MQ=255)
            ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                      <  2:2360267/149‑1 (MQ=255)
             caccaTCCAGTGCAGGCGCTCCGATTATAAGCTGTCAAAATTGAGAATTCCGATCCCAAATGCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                     <  1:167309/149‑1 (MQ=18)
              aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                    <  2:3846510‑M1/92‑1 (MQ=255)
               acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCTCTCCCGGGGGCTTTGTCCTAc                                                                                                                                   <  1:1302424‑M1/93‑1 (MQ=255)
               acgaaaggcccagtctttcgactgagcctttcgttgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGAGGTGTCCTAc                                                                                                                                   <  2:788832‑M1/93‑1 (MQ=255)
                caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                                  <  1:2739514/149‑1 (MQ=255)
                 gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATCCGCTCCCGGCGGAGTTGGCCTACTc                                                                                                                                 <  2:484454‑M1/95‑1 (MQ=255)
                  acctgggCAGGAGCTCCGATTATAAGCTGACAAAAATTACAATTCCGACCCCCACTTCCTGGCAGTTTTTGGGGGGCGTCGTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                                <  1:2816468/143‑1 (MQ=17)
                   ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                               <  2:1495115/149‑1 (MQ=255)
                       tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                           <  2:2934794/149‑1 (MQ=255)
                        ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                          <  2:1016148‑M1/102‑1 (MQ=255)
                         cAGGAGCTCCGATGAGAAGCTGTCAAACATGTGACTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                         <  2:4524250/149‑1 (MQ=255)
                           ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                       <  1:326007/149‑1 (MQ=255)
                            tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGCGAGCGtt                                                                                                                      <  2:2826236‑M1/106‑1 (MQ=255)
                            tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                      >  1:3846510‑M1/44‑149 (MQ=255)
                            gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                      >  2:2293967/1‑149 (MQ=255)
                            gAGCTCCGATGAGAAGCTGTCAACCATGAGACTTCCGAGCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                      <  2:4678573/149‑1 (MQ=255)
                             ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGTCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTc                                                                                                                     <  2:484755‑M1/107‑1 (MQ=255)
                             ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGTCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTATGTCATACTCAGGAGagcgatc                                                                                                                     <  2:1598599‑M1/107‑4 (MQ=255)
                             aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTAGCAACGTTCAAATACGCTCCCGGCGGATTTGTACTACTGAGGAGAACGtta                                                                                                                     >  1:3102378/1‑148 (MQ=21)
                              tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                    >  2:2722865‑M1/42‑149 (MQ=255)
                              gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCACCGGGCAAATCCGCGCCCCGCGGCTTTGTCCTACTCAGGAGAGCGtgag                                                                                                                    >  1:881977/1‑146 (MQ=18)
                                  cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCGTAAAATCCGCAACCCGCCGCTTTGTCCAACTCAGGTGAGCCTTCACCGa                                                                                                                >  1:4094301/1‑149 (MQ=18)
                                   gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                               <  1:2618312/149‑1 (MQ=255)
                                   actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGGTCAAAACCCCTCCCCGCGGATTTTTCCAACACAGGCTAGCGTTAAccgcc                                                                                                               <  1:4704755‑M1/113‑3 (MQ=255)
                                     tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTACAACTCAGGAGAGCGTTCACCGacaa                                                                                                             <  2:1804288‑M1/115‑1 (MQ=255)
                                      gagcctttcgatctacatgatgcccggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCATCCGGTCCGTTGCTTAGCAACTATCAAATCCGCTCACGTCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                            >  2:3530260‑M1/34‑149 (MQ=255)
                                       atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                           <  2:2018676/149‑1 (MQ=255)
                                       atAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                           <  2:2756262/149‑1 (MQ=255)
                                       agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                           <  2:4263171‑M1/117‑1 (MQ=255)
                                          aGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                                        <  2:4104528/149‑1 (MQ=255)
                                           tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                       >  2:3458560‑M1/29‑149 (MQ=255)
                                           gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                       <  2:4427911/149‑1 (MQ=255)
                                             tGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                     <  1:1862134/149‑1 (MQ=255)
                                                 aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                                 <  2:4134254/149‑1 (MQ=255)
                                                  aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                <  2:1890509/149‑1 (MQ=255)
                                                   tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                               <  2:519941‑M1/129‑1 (MQ=255)
                                                    atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAAAACCCAACCGGCGGATTTGTCCTACACAGGAGATCGTTCACCCCCAAACAAAAGataaaaac                                                                                              <  1:4641480‑M1/130‑3 (MQ=255)
                                                    atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACAACAAATAAAACg                                                                                              <  1:1137990‑M1/130‑1 (MQ=255)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                            <  2:850524/149‑1 (MQ=255)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAACTTCAAATACGCTCCCCGCGGATTTTTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGaa                                                                                            >  1:1020639/1‑149 (MQ=14)
                                                      tGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                            <  2:2880206/149‑1 (MQ=255)
                                                       gagaATTCCGATCCCCATTTCGTGGCTGTTTTTGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                           <  1:1786042/149‑1 (MQ=14)
                                                        agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTGGCAACGGTGCAAACCCCGCCCCGCCGCATTTTCCGAAACAGGGGAGCGGTCACCGAACCACAACAGCTaaaaacaaaa                                                                                          >  1:883186/1‑143 (MQ=9)
                                                           aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                       >  2:3504501/1‑149 (MQ=255)
                                                            ttCCGATCCCGAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                      <  2:2538986/149‑1 (MQ=255)
                                                             ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                     <  2:1173793‑M1/139‑1 (MQ=255)
                                                             ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                     <  2:4663357‑M1/139‑1 (MQ=255)
                                                              cctaTCCCAAATTCGTGGCTGTTTTGGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTATTTGTCCTACTCAGGTGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                    <  1:1757393/146‑1 (MQ=11)
                                                              ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                    >  2:86147/1‑149 (MQ=40)
                                                               cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTGCCCCGACAAACAACAGATAAAACGTAAGGCCCAGt                                                                                   >  2:2280625/1‑149 (MQ=255)
                                                                 ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTACCGGCGGATTTTTCCTACACAGGAGAGCGTTCACCGACAAACAACAGATAAAACCAAATGCCCAGTCt                                                                                 <  1:3162039‑M1/143‑1 (MQ=255)
                                                                 aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                 <  2:3590867/149‑1 (MQ=38)
                                                                  gaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACCCTCCCGGCGGATTTTTCCTACTCAGGCGAGCGTTCACCCACAAACAACACATAAAAACAAATGGCCACTCtt                                                                                <  1:3835019‑M1/144‑1 (MQ=255)
                                                                   aattAATTCCTGGCAGTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                               >  1:788832‑M1/5‑149 (MQ=255)
                                                                     ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                                 >  1:3528965‑M1/3‑113 (MQ=255)
                                                                     ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                             <  2:3951133‑M1/147‑1 (MQ=255)
                                                                     ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACCTTAAAATACGCACCCGGCGGATTTTTACTACACAGGAGAGCGTTCACCGACCAAAAACAGAAAAAAACAAAGGCCCCCTCTTTCg                                                                             >  1:268401/1‑149 (MQ=11)
                                                                     ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTACGAATCAGGAGAGCGTTGACCGACAAACAACAGATAAAACGAAAGGCCCCGTCgtttg                                                                             >  1:1179310/1‑147 (MQ=11)
                                                                       aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAATGCCCCGTCTTTCGAc                                                                           >  1:1550023/1‑149 (MQ=11)
                                                                       aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGGGGATTTGTCCGACTCAGGAGAGCGTTCACCGGCAAACAACAGATAAAACGAAAGGCCCAGTGTTTGGAc                                                                           >  2:2107195/1‑149 (MQ=11)
                                                                        aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                          <  2:309379/149‑1 (MQ=32)
                                                                          tCCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                        <  1:124110/149‑1 (MQ=11)
                                                                          tCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                        <  1:4774920/149‑1 (MQ=11)
                                                                                                    gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                              <  1:4738670/149‑1 (MQ=11)
                                                                                                    gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                              <  1:484676/149‑1 (MQ=11)
                                                                                                       cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                           <  2:1206818/149‑1 (MQ=17)
                                                                                                       cGCCACCCTCCGGGCCGTTGCTTCGCAACGTGCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                           <  2:1918684/149‑1 (MQ=17)
                                                                                                        gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                          <  1:2940403/149‑1 (MQ=33)
                                                                                                        gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                          <  1:2100977/149‑1 (MQ=33)
                                                                                                         ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                         <  2:1958703/149‑1 (MQ=21)
                                                                                                         ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                         <  1:765286/149‑1 (MQ=34)
                                                                                                         ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                         <  2:2463399/149‑1 (MQ=34)
                                                                                                         ccACCCTCCGGGCCGTTGCGTAGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCAGGCGGGCGGTCACCCCCACACAACAGATAAAACGAAAGGCCCACTCTTTCGACGGTGCCGGTCGGGGTGTTTGATGCCTgggggtgg                                         >  2:2898756/1‑145 (MQ=2)
                                                                                                           aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGCTAAAACGACAGGCCCAGTCTTTCGTCTGAGCCGTTCGTTTTTTTTGCTGCGTGggggattgg                                       >  2:2616023/1‑149 (MQ=11)
                                                                                                           aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCAACCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCACTCTGTCGACTGAGCCGGTCGTTGTATTTGATGCCTGgggaattgg                                       >  2:1305550/1‑149 (MQ=11)
                                                                                                            cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattttt                                      >  2:678469/1‑146 (MQ=21)
                                                                                                              cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                        >  2:2781088/1‑145 (MQ=35)
                                                                                                                       gTTGCTTCGCAACTTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                           <  1:2004893/149‑1 (MQ=255)
                                                                                                                       gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCCCCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGGTTTCTTTGATGCCGGGGGACTTGGGAATTAAGTCtt                           >  2:2578369/1‑149 (MQ=14)
                                                                                                                          gCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                        <  1:2227489/149‑1 (MQ=255)
                                                                                                                                 aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa                 <  1:93457/149‑1 (MQ=255)
                                                                                                                                 aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa                 >  2:194937/1‑149 (MQ=255)
                                                                                                                                  aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg                >  1:3572154/1‑149 (MQ=255)
                                                                                                                                   cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                           <  1:2373103/121‑1 (MQ=33)
                                                                                                                                              cGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCCCCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAt     <  1:2445230/149‑1 (MQ=255)
                                                                                                                                              cGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACCGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAATTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAt     <  2:1378236/149‑1 (MQ=255)
                                                                                                                                               gCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa    >  2:4551927/1‑149 (MQ=255)
                                                                                                                                                 tCCCGGCGGAGTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGCTGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACg  >  1:588535/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACG  >  CP009273/360033‑360325

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: