Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,172 T→A 5.3% intergenic (‑106/+489) 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,1720TA 5.3% 626.0 / 2.8 209intergenic (‑106/+489)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (59/132);  new base A (12/4);  total (72/137)
Fisher's exact test for biased strand distribution p-value = 6.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCG  >  CP009273/360024‑360318
                                                                                                                                                    |                                                                                                                                                  
gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTATATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGGTGAAAt                                                                                                                                                    <  1:1215769‑M1/69‑1 (MQ=255)
 cAGCGGCTTACCAACCCGCGCCACCATCCCGTGCAGGAGCTCCTTTTATATGGTGAAAAACATGCTAATTCCGACCCCAAATCCCTGGCATTTTTGGGGGGGCGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                                   <  1:1296124/149‑1 (MQ=16)
    cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                                <  2:4265048/149‑1 (MQ=255)
        ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                            <  2:3199621/149‑1 (MQ=255)
        acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACTTTCAAATCCGCGccc                                                                                                                                            <  1:1736170‑M1/77‑1 (MQ=255)
          aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCAACcgg                                                                                                                                          >  1:305775/1‑149 (MQ=255)
           aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                                         <  2:3069253‑M1/80‑1 (MQ=255)
            cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCTGGGCCGTTGGTTCGCAACCTTCAAAACCCCTCCccgcg                                                                                                                                        >  1:3423158/1‑149 (MQ=37)
            cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCTGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                        >  1:3889490/1‑149 (MQ=255)
              tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                      <  1:1666525/149‑1 (MQ=255)
              tCCAGCGCCACCAGCCAGGGCCGGAGCTCCTATGATAAGCTGTCACCCATGAGATTTCCGCGCCCCAATGCCTGGCTGTTTATGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                      <  2:2969132/149‑1 (MQ=18)
                cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                    <  2:1057398/149‑1 (MQ=255)
                  gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTCAAAAACGCTCCCGGCGGCTTTg                                                                                                                                  >  1:4271106/1‑149 (MQ=25)
                   cgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                                 <  2:2363820/149‑1 (MQ=255)
                   ataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggacttAATTCCTGGCAGGTTATGGCGGGCGTCCGGCCCGCCCCCCTCCGGGCCGTTGCTTCGCCACGGTCACATCCGCTCCCGGCGGATgtgt                                                                                                                                 <  2:3167483‑M1/88‑1 (MQ=255)
                    taaaacgaatggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTGCGCTACGTGCCATTCCGCTCCcggcggagtggtc                                                                                                                                <  2:3227992‑M1/89‑7 (MQ=255)
                    gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGCTTTGTc                                                                                                                                >  2:3692810/1‑149 (MQ=255)
                     ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTGAAATCCGCTCCCGGCGGATTTGtgc                                                                                                                               >  1:2100904/1‑147 (MQ=255)
                     aaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCAACGTTAAAATACGCTCCCGGCGGATTtgtac                                                                                                                               <  1:3261095‑M1/90‑3 (MQ=255)
                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATACGCTCCCGGCGGTATTGTCCt                                                                                                                              >  1:4348999/1‑149 (MQ=255)
                        ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                            <  2:612541/149‑1 (MQ=255)
                        ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACCTTCAAATACGCTCCCCGCGGAGTTGTCCTAc                                                                                                                            >  1:1382203/1‑149 (MQ=25)
                          aTCCAGTGCAGGAGCTCCGTTGATACGCTGTCACACATGAGAATCCCGACCCCCAATTCCGGGCTTTTTTGGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                                          <  1:1502627/149‑1 (MQ=18)
                          aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCAGCTCCCGGCGGATTTGTCCTACTc                                                                                                                          <  2:1782705/149‑1 (MQ=255)
                           aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTAGTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCAACCGGGCCGTGGCTTAGCAACGTTAAATTGCGCTCCCGGGGGAGTTGTCctacgca                                                                                                                         <  2:253133‑M1/96‑4 (MQ=255)
                              aGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                      <  1:3882380/149‑1 (MQ=255)
                                tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                    <  2:2384372/149‑1 (MQ=255)
                                 gccggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACCTTCAAAAACGCACCCCGCGGATTTGTCCGACGCAGgaggg                                                                                                                   >  1:3866592/4‑147 (MQ=21)
                                    ggAGCTCCGATGATACGCTGTCAAACAGGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                <  1:306151/149‑1 (MQ=255)
                                    ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                >  2:645097/1‑149 (MQ=255)
                                      aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                              <  1:1415476/149‑1 (MQ=255)
                                          ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                          <  2:404807/149‑1 (MQ=255)
                                              tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                      <  1:681313/149‑1 (MQ=255)
                                              tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                      >  2:1723844‑M1/35‑149 (MQ=255)
                                               gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCGCCCGGCGGGTTTGTCCTACTCAGGAGAGCGTTGACCGccaaa                                                                                                     >  2:3023359/1‑149 (MQ=18)
                                               gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAAACGCTCCCCGCGGAGTTGTCCGACTCAGGAGAGCGTTCACCCacaaa                                                                                                     >  2:817732/1‑149 (MQ=18)
                                                      tcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                              >  1:3261056‑M1/27‑149 (MQ=255)
                                                      tcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGACCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGATCGTTCACCGACAAACAACAGa                                                                                              >  2:818787‑M1/27‑149 (MQ=255)
                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                              <  2:3091798/149‑1 (MQ=255)
                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                             <  2:3241591/149‑1 (MQ=255)
                                                         cAAACTTGAGATTTCCGATCCCCATTTCCTGGCTGTTTTGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                           <  1:2138529/149‑1 (MQ=14)
                                                          aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                          <  2:4378333/149‑1 (MQ=255)
                                                           aaCATGAGATTTCCGATCCCCACTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                         <  1:2406801/149‑1 (MQ=255)
                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                         <  2:3883327/149‑1 (MQ=255)
                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAAGTTCAAAAAAGCGCCCGGCCGGTTTTTGCTACACAGGAGAGGGTTCACCGACAAACAACAGATaaaaag                                                                                       >  1:1004222/1‑147 (MQ=17)
                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACGCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                       <  2:324616‑M1/130‑1 (MQ=255)
                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCGTCAAATCCCCTCCCCGCCGATTTGTCCGACACACGAGAGCGTTCACCCACAAACAACACAAAAAacc                                                                                       <  1:2912872‑M1/130‑2 (MQ=255)
                                                             atttgatgcctggggaattAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAACg                                                                                       <  2:1019065‑M1/130‑1 (MQ=255)
                                                              tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAgataaaaaga                                                                                      <  2:1949939‑M1/131‑4 (MQ=255)
                                                              aggagaAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                                      <  2:1300609/147‑1 (MQ=255)
                                                              aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                                      <  2:2719994/149‑1 (MQ=255)
                                                                tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACACAAAAAACGaaa                                                                                    <  2:1277382‑M1/133‑1 (MQ=255)
                                                                gttatttCCGATACAAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                    <  2:2100904/144‑1 (MQ=255)
                                                                 gttggagggggaattAATTCGTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                   >  1:3916785‑M1/16‑149 (MQ=255)
                                                                  gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGGTGGCCACGGTGAAAAAACCTACCCGCGGGTTTGTGCGAATGAGGCGAGCGTTGACCGAGCAAAAACACAAAAAACGgaatg                                                                                  >  1:3739242/1‑147 (MQ=9)
                                                                  atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGGTTTGTCCTACGCAGGGGCGCGGTCACCGACAACCAACAGATAAAACGAAAgg                                                                                  <  2:326174‑M1/135‑1 (MQ=255)
                                                                    aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                <  2:2786968/149‑1 (MQ=255)
                                                                     gatgtgtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                               >  2:1310773‑M1/12‑149 (MQ=255)
                                                                      tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                              <  2:1325889/149‑1 (MQ=40)
                                                                      tCCGAGCCCCAATTCCTGGCTGTTTTGGGGGGGCGACCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                              <  1:1887275/149‑1 (MQ=11)
                                                                        ggATCCCCAATTCCTGGCAGTTTTGGGCGGGCGTCCTGCCTGCCACCCTCCGGGCCGTTGCTTCGCAACTTTCAAATCCGATCCGGGGTTATTTGTCCTACTATGGAGAGCTTTCACCGACAAACAACATATAAAACGAAAGGCCCAGt                                                                            <  1:3081560/148‑1 (MQ=11)
                                                                           cccccAATTCCTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                         <  1:913294/148‑1 (MQ=11)
                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTGCACCGCCACACAACAGATAAAACGAAAGGCCCCGTCttt                                                                        <  2:3953815‑M1/145‑1 (MQ=255)
                                                                             cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                       <  2:343007/149‑1 (MQ=35)
                                                                             attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                       >  1:325217‑M1/4‑149 (MQ=255)
                                                                              ttaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACCCACCCCGCGGATTTTTACTACAAAGGAGAGCGTTAACCGACAAAAAAAAAATAAAACTAAAGGCCCCGTCTTTCg                                                                      >  1:2178522/3‑149 (MQ=11)
                                                                              ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                      <  2:4373831‑M1/147‑1 (MQ=255)
                                                                              ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAAAAACGCTCCCGGCGGGTTTGTCCAAAAAAGGAGAGCGTTGAACGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCg                                                                      >  1:305100/1‑149 (MQ=11)
                                                                               tgaggcgtGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                     <  1:2706809/142‑1 (MQ=11)
                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCCGTCTTTCGa                                                                     <  2:2937311‑M1/148‑1 (MQ=255)
                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTACTACTCAGGAGAGCGTTAACCCACAAACAACAAAAAAAACGAAAGGCCCAGTCTTTAGAc                                                                    >  1:2849131/1‑149 (MQ=11)
                                                                                aaTTCCGGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                    <  2:2615702/149‑1 (MQ=17)
                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTAGACTg                                                                  >  2:1407629/1‑149 (MQ=14)
                                                                                   tCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                 <  2:1185296/149‑1 (MQ=11)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                    <  1:3191435/130‑1 (MQ=11)
                                                                                                gggggggcgTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                                    <  1:1268364/145‑1 (MQ=0)
                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                       <  2:1480688/149‑1 (MQ=11)
                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCAACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                      <  1:4284605/149‑1 (MQ=14)
                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                     <  2:2452853/149‑1 (MQ=32)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                   <  2:908274/149‑1 (MQ=33)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                   <  1:2758972/149‑1 (MQ=33)
                                                                                                                  ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                  <  2:1895108/149‑1 (MQ=34)
                                                                                                                   cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                 <  2:540351/149‑1 (MQ=35)
                                                                                                                     cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg                               >  2:2047182/1‑149 (MQ=37)
                                                                                                                           gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAATGCCCCGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGacttaa                         >  2:1038065/1‑149 (MQ=255)
                                                                                                                            ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat                        <  1:1407629/149‑1 (MQ=255)
                                                                                                                             gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt                       <  1:4373831/149‑1 (MQ=255)
                                                                                                                              ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTGGATTTGATGCCGGGGGCAGTGGGCATTAATTc                      >  1:836115/1‑149 (MQ=11)
                                                                                                                                ggtGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTACACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                    <  1:2291214/147‑1 (MQ=255)
                                                                                                                                ggagCTTCGCAACGTTCACATCCGCTCCCGGCGGATTTGTCCTCCTCCGGAGAGCGTTCACCGACAAACAACAGATCAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                    <  1:253133/146‑1 (MQ=14)
                                                                                                                                   gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGTGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                 >  2:2835295/1‑149 (MQ=255)
                                                                                                                                        gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGa            <  1:2814277/149‑1 (MQ=255)
                                                                                                                                         cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaa           >  1:728570/1‑149 (MQ=255)
                                                                                                                                          aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTGTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa          >  2:178535/1‑149 (MQ=255)
                                                                                                                                           aCGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg         >  2:1626687/1‑149 (MQ=255)
                                                                                                                                             gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAggg       <  1:1317752/149‑1 (MQ=255)
                                                                                                                                             gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAggg       <  1:2047182/149‑1 (MQ=255)
                                                                                                                                              ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc      >  1:3463933/1‑149 (MQ=255)
                                                                                                                                              ttCAAATCCGCTCCCGGCGGATGGGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc      <  2:3758154/149‑1 (MQ=255)
                                                                                                                                                cAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt    >  1:1815307/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCG  >  CP009273/360024‑360318

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: