Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,915,555 A→C 16.5% Q359P (CAG→CCG)  rsmF → 16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,915,5550AC16.5% 742.6 / 162.0 301Q359P (CAG→CCG) rsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (121/130);  new base C (20/30);  total (141/160)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 307 total aligned reads displayed.

AAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCC  >  CP009273/1915408‑1915703
                                                                                                                                                   |                                                                                                                                                    
aaaCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAg                                                                                                                                                     >  1:2337554/1‑149 (MQ=255)
  aCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCg                                                                                                                                                   >  1:3597165/1‑149 (MQ=255)
    tCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAACGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTg                                                                                                                                                 <  2:1801678/149‑1 (MQ=255)
     cAGGCGAGTCCCGCCTTACCCGCCCCCAAAGACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGa                                                                                                                                                <  2:2743623/149‑1 (MQ=255)
      aGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGAc                                                                                                                                               <  2:1932562/149‑1 (MQ=255)
      aGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGAc                                                                                                                                               <  2:3343049/149‑1 (MQ=255)
             tCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAAc                                                                                                                                        >  1:1817391/1‑149 (MQ=255)
              cccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACt                                                                                                                                       <  2:3910315/149‑1 (MQ=255)
              cccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACt                                                                                                                                       <  2:1798824/149‑1 (MQ=255)
               ccGCCTTACCCGCCCCCAAATACACAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTCCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtg                                                                                                                                      >  2:1941291/1‑149 (MQ=255)
                  ccTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTgg                                                                                                                                   <  2:1161193/149‑1 (MQ=255)
                   cTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgt                                                                                                                                  <  2:2465475/149‑1 (MQ=255)
                   cTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgt                                                                                                                                  <  1:4290578/149‑1 (MQ=255)
                    ttACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgtt                                                                                                                                 >  2:3514835/1‑149 (MQ=255)
                    ttACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgtt                                                                                                                                 <  2:1840202/149‑1 (MQ=255)
                     tACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgttg                                                                                                                                >  1:2540482/1‑149 (MQ=255)
                         cGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTcc                                                                                                                            <  1:1340114/149‑1 (MQ=255)
                         cGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTcc                                                                                                                            <  1:520072/149‑1 (MQ=255)
                         cGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAAACTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTcc                                                                                                                            <  1:2965072/149‑1 (MQ=255)
                          gCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTccc                                                                                                                           <  2:1716795/149‑1 (MQ=255)
                             cccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGt                                                                                                                        >  2:2312833/1‑149 (MQ=255)
                               cAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTgg                                                                                                                      >  1:1769145/1‑149 (MQ=255)
                                 aaTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGc                                                                                                                    <  1:2078559/149‑1 (MQ=255)
                                 aaTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGc                                                                                                                    >  2:708926/1‑149 (MQ=255)
                                  aTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCa                                                                                                                   <  1:3401796/149‑1 (MQ=255)
                                    aCAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAtt                                                                                                                 >  1:546581/1‑149 (MQ=255)
                                    aCAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAtt                                                                                                                 >  1:3709691/1‑149 (MQ=255)
                                      aaaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGa                                                                                                               <  1:4656012/149‑1 (MQ=255)
                                       aaGTCGGTAATTTTCCGGTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGaa                                                                                                              <  2:2934681/149‑1 (MQ=255)
                                             gTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTg                                                                                                        >  2:2828350/1‑149 (MQ=255)
                                              tAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGa                                                                                                       <  2:2125005/149‑1 (MQ=255)
                                               aaTTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGAt                                                                                                      <  2:1524516/149‑1 (MQ=255)
                                               aaTTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGAt                                                                                                      <  1:1009278/149‑1 (MQ=255)
                                                 ttttCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCg                                                                                                    <  1:1144829/149‑1 (MQ=255)
                                                  tttCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCgg                                                                                                   >  1:618803/1‑149 (MQ=255)
                                                  tttCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCgg                                                                                                   <  1:3271287/149‑1 (MQ=255)
                                                      cGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTaaa                                                                                               <  1:4393987/149‑1 (MQ=255)
                                                       gTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAg                                                                                              >  1:2752541/1‑149 (MQ=255)
                                                        ttCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCTGGTGGGCATTGAAGCCCTGATCGGTAAAGt                                                                                             >  2:2208112/1‑149 (MQ=255)
                                                           aGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCg                                                                                          <  2:1464184/149‑1 (MQ=255)
                                                               cGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAttt                                                                                      >  2:4262188/1‑149 (MQ=255)
                                                                  tGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtct                                                                                   <  1:2705021/149‑1 (MQ=255)
                                                                     aaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCgg                                                                                >  2:2059188/1‑149 (MQ=255)
                                                                     aaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCgg                                                                                >  1:3724817/1‑149 (MQ=255)
                                                                      aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                               >  2:1028197/1‑149 (MQ=255)
                                                                      aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                               >  1:2537104/1‑149 (MQ=255)
                                                                       gATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGtt                                                                              <  1:1076594/149‑1 (MQ=255)
                                                                        aTCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTg                                                                             >  1:2035609/1‑149 (MQ=255)
                                                                           gcgAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTgggg                                                                          >  1:540867/1‑149 (MQ=255)
                                                                           gcgAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTgggg                                                                          >  2:4066582/1‑149 (MQ=255)
                                                                             gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAt                                                                        <  2:4053852/149‑1 (MQ=255)
                                                                             gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAt                                                                        <  1:2569409/149‑1 (MQ=255)
                                                                             gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAt                                                                        <  1:3526988/149‑1 (MQ=255)
                                                                                  tGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAAc                                                                   >  2:4465922/1‑149 (MQ=255)
                                                                                   ggACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACt                                                                  >  1:2256608/1‑149 (MQ=255)
                                                                                   ggACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACt                                                                  >  1:2842067/1‑149 (MQ=255)
                                                                                         aTTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCga                                                            <  1:3090629/149‑1 (MQ=255)
                                                                                          ttCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgag                                                           <  1:368500/149‑1 (MQ=255)
                                                                                           tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                          <  2:577529/149‑1 (MQ=255)
                                                                                           tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                          <  2:3263176/149‑1 (MQ=255)
                                                                                               cAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGca                                                      <  1:3492323/149‑1 (MQ=255)
                                                                                                aGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGcac                                                     <  1:3328069/149‑1 (MQ=255)
                                                                                                aGGCGGATACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGcac                                                     <  2:1793201/149‑1 (MQ=255)
                                                                                                   cggcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCacaac                                                  <  2:634435/149‑1 (MQ=255)
                                                                                                     gcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCAcaacaa                                                >  2:3027194/1‑149 (MQ=255)
                                                                                                       tACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAg                                              >  2:3301911/1‑149 (MQ=255)
                                                                                                        aCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAgg                                             >  1:1225906/1‑149 (MQ=255)
                                                                                                         cAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGATGGGGATTAAACTTGCCGAGACGCACAACAAAGGt                                            >  1:1948825/1‑149 (MQ=255)
                                                                                                          aGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGtt                                           <  1:563583/149‑1 (MQ=255)
                                                                                                           ggTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTa                                          <  2:1654586/149‑1 (MQ=255)
                                                                                                           ggTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTa                                          >  2:3578893/1‑149 (MQ=255)
                                                                                                            gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTAt                                         >  2:535204/1‑149 (MQ=255)
                                                                                                             tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                        <  1:3580309/149‑1 (MQ=255)
                                                                                                               ttGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGc                                      <  1:3908197/149‑1 (MQ=255)
                                                                                                               ttGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGc                                      <  2:1956794/149‑1 (MQ=255)
                                                                                                                  gCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTgg                                   <  2:4171107/149‑1 (MQ=255)
                                                                                                                   cTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGgc                                  <  1:4065810/149‑1 (MQ=255)
                                                                                                                       aaCTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGTCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAATGTTATCGCTGgcagca                              <  1:3601958/149‑1 (MQ=255)
                                                                                                                        aCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCAt                             >  2:3112314/1‑149 (MQ=255)
                                                                                                                        aCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCAt                             <  2:1989885/149‑1 (MQ=255)
                                                                                                                           gggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGaa                          >  1:2392196/1‑149 (MQ=255)
                                                                                                                           gggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGaa                          >  1:2956398/1‑149 (MQ=255)
                                                                                                                            ggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                         >  1:1577602/1‑149 (MQ=255)
                                                                                                                            ggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                         >  1:4477077/1‑149 (MQ=255)
                                                                                                                                  aaaCCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                   <  2:4485849/149‑1 (MQ=255)
                                                                                                                                     ccTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTg                >  2:3743819/1‑149 (MQ=255)
                                                                                                                                      cTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGc               >  2:117827/1‑149 (MQ=255)
                                                                                                                                        gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctc             <  1:3027194/149‑1 (MQ=255)
                                                                                                                                        gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctc             <  2:2084514/149‑1 (MQ=255)
                                                                                                                                        gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctc             <  1:688105/149‑1 (MQ=255)
                                                                                                                                         cgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctct            >  2:4675818/1‑149 (MQ=255)
                                                                                                                                           ccTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTg          >  1:948549/1‑149 (MQ=255)
                                                                                                                                              ctGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCt       <  2:4603957/149‑1 (MQ=255)
                                                                                                                                               tGGCAGCGTGACAAAGATCTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTc      <  1:2641335/149‑1 (MQ=255)
                                                                                                                                               tGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTc      <  1:1693971/149‑1 (MQ=255)
                                                                                                                                                ggCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCa     >  1:4094369/1‑149 (MQ=255)
                                                                                                                                                 gCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAc    >  2:136988/1‑149 (MQ=255)
                                                                                                                                                   aGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAccc  >  2:2559196/1‑149 (MQ=255)
                                                                                                                                                   aGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAccc  >  2:795387/1‑149 (MQ=255)
                                                                                                                                                   aGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATGTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAccc  <  2:2633692/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCC  >  CP009273/1915408‑1915703

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: