Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,915,555 A→C 14.4% 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,5550AC14.4% 791.5 / 133.4 296Q359P (CAG→CCG) rsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (101/152);  new base C (27/16);  total (128/168)
Fisher's exact test for biased strand distribution p-value = 7.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

AAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCC  >  CP009273/1915407‑1915703
                                                                                                                                                    |                                                                                                                                                    
aaaaCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATCCAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCa                                                                                                                                                      <  2:4052745/149‑1 (MQ=255)
   aCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCg                                                                                                                                                   >  2:1571581/1‑149 (MQ=255)
    cTCAGGCGATTCCCGCCTTACCCGCCCCCAAATCCAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGt                                                                                                                                                  <  2:3068561/149‑1 (MQ=255)
     tCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTg                                                                                                                                                 <  2:1324932/149‑1 (MQ=255)
      cAGGCGATTCCCGCCTTACCCGCCCCCAAATCCAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGa                                                                                                                                                <  2:3052285/149‑1 (MQ=255)
      cAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGa                                                                                                                                                <  1:4299672/149‑1 (MQ=255)
       aGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGAc                                                                                                                                               >  1:601672/1‑149 (MQ=255)
            aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGa                                                                                                                                          >  2:464552/1‑149 (MQ=255)
                ccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGACAACCTGCGCCTCTGGCAGCGTGACACAGAACtg                                                                                                                                      >  2:3949839/1‑149 (MQ=255)
                ccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtg                                                                                                                                      >  2:4041585/1‑149 (MQ=255)
                ccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtg                                                                                                                                      <  2:2672172/149‑1 (MQ=255)
                  gCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtgtg                                                                                                                                    <  1:1807121/149‑1 (MQ=255)
                    cTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgt                                                                                                                                  >  2:25179/1‑149 (MQ=255)
                     ttACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCATCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgtt                                                                                                                                 >  1:3198497/1‑149 (MQ=255)
                     ttACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgtt                                                                                                                                 <  1:2981112/149‑1 (MQ=255)
                      tACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgttg                                                                                                                                >  1:3858245/1‑149 (MQ=255)
                        cccGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgttgtt                                                                                                                              <  2:407315/149‑1 (MQ=255)
                         ccGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTc                                                                                                                             >  1:4139459/1‑149 (MQ=255)
                         ccGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTc                                                                                                                             >  2:1001357/1‑149 (MQ=255)
                          cGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTcc                                                                                                                            <  2:1446438/149‑1 (MQ=255)
                           gCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTccc                                                                                                                           <  1:1645979/149‑1 (MQ=255)
                            cccccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCg                                                                                                                          >  1:2253711/1‑149 (MQ=255)
                              cccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGt                                                                                                                        <  2:3403814/149‑1 (MQ=255)
                                cAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTgg                                                                                                                      <  1:4599780/149‑1 (MQ=255)
                                 aaaTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTggg                                                                                                                     >  1:5158305/1‑149 (MQ=255)
                                   aTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCa                                                                                                                   >  2:4557413/1‑149 (MQ=255)
                                       aaaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGa                                                                                                               <  2:2118305/149‑1 (MQ=255)
                                        aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATTGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGaa                                                                                                              >  2:257999/1‑149 (MQ=255)
                                         aGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAg                                                                                                             <  1:1162393/149‑1 (MQ=255)
                                             ggTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCt                                                                                                         >  2:3667692/1‑149 (MQ=255)
                                               tAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGa                                                                                                       <  1:4969806/149‑1 (MQ=255)
                                                  ttttCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCg                                                                                                    >  1:3133356/1‑149 (MQ=255)
                                                  ttttCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCg                                                                                                    >  1:1287042/1‑149 (MQ=255)
                                                     tCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTa                                                                                                 <  2:2047181/149‑1 (MQ=255)
                                                       cGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTaaa                                                                                               <  1:1310796/149‑1 (MQ=255)
                                                          tCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTc                                                                                            <  1:3686946/149‑1 (MQ=255)
                                                            aGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCg                                                                                          >  2:2918374/1‑149 (MQ=255)
                                                              cccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAt                                                                                        <  2:447610/149‑1 (MQ=255)
                                                              cccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAt                                                                                        <  1:2408969/149‑1 (MQ=255)
                                                                cGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAttt                                                                                      >  2:1382660/1‑149 (MQ=255)
                                                                  gTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtc                                                                                    <  2:2384424/149‑1 (MQ=255)
                                                                   tGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtct                                                                                   >  1:4753593/1‑149 (MQ=255)
                                                                    gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtctc                                                                                  <  2:2819155/149‑1 (MQ=255)
                                                                    gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtctc                                                                                  <  2:45334/149‑1 (MQ=255)
                                                                       aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                               >  2:2227839/1‑149 (MQ=255)
                                                                         aTCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTg                                                                             <  1:1736803/149‑1 (MQ=255)
                                                                              gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAt                                                                        <  2:4274247/149‑1 (MQ=255)
                                                                                   tGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACACAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAAc                                                                   <  1:4847788/149‑1 (MQ=255)
                                                                                   tGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAAc                                                                   >  2:1883211/1‑149 (MQ=255)
                                                                                   tGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAAc                                                                   <  1:2665299/149‑1 (MQ=255)
                                                                                       cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGc                                                               <  1:3341220/149‑1 (MQ=255)
                                                                                       cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGc                                                               >  1:4664468/1‑149 (MQ=255)
                                                                                         aaTTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCg                                                             >  2:3027446/1‑149 (MQ=255)
                                                                                          aTTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCga                                                            >  2:3149559/1‑149 (MQ=255)
                                                                                          aTTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCga                                                            >  2:4466000/1‑149 (MQ=255)
                                                                                           ttCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgag                                                           <  1:783301/149‑1 (MQ=255)
                                                                                            tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAATCACTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                          >  2:734638/1‑149 (MQ=255)
                                                                                            tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                          <  2:3771215/149‑1 (MQ=255)
                                                                                            tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                          >  1:319275/1‑149 (MQ=255)
                                                                                               tCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGc                                                       >  1:2142295/1‑149 (MQ=255)
                                                                                               tCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGc                                                       <  1:5046298/149‑1 (MQ=255)
                                                                                                     ggcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCacaaca                                                 >  2:3077598/1‑149 (MQ=255)
                                                                                                     ggcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCacaaca                                                 <  2:1748781/149‑1 (MQ=255)
                                                                                                       cTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACaaa                                               <  2:1034799/149‑1 (MQ=255)
                                                                                                        tACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAg                                              <  1:3457381/149‑1 (MQ=255)
                                                                                                         aCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAgg                                             <  1:3645038/149‑1 (MQ=255)
                                                                                                         aCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAgg                                             >  1:1448956/1‑149 (MQ=255)
                                                                                                          cAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGt                                            <  2:2631187/149‑1 (MQ=255)
                                                                                                             gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTAt                                         <  2:471659/149‑1 (MQ=255)
                                                                                                             gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTAt                                         <  2:4808235/149‑1 (MQ=255)
                                                                                                             gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTAt                                         <  2:3359179/149‑1 (MQ=255)
                                                                                                              tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                        <  2:135372/149‑1 (MQ=255)
                                                                                                                ttGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGc                                      <  1:3564286/149‑1 (MQ=255)
                                                                                                                 tGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCt                                     <  2:1512691/149‑1 (MQ=255)
                                                                                                                  ggCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTg                                    <  2:2323306/149‑1 (MQ=255)
                                                                                                                    cTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGgc                                  <  2:849613/149‑1 (MQ=255)
                                                                                                                    cTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGgc                                  >  1:3953070/1‑149 (MQ=255)
                                                                                                                         aCTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCAt                             <  1:2918374/149‑1 (MQ=255)
                                                                                                                          cTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATg                            >  1:2884518/1‑149 (MQ=255)
                                                                                                                             ggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                         >  2:706618/1‑149 (MQ=255)
                                                                                                                               aTGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCg                       >  2:5104081/1‑149 (MQ=255)
                                                                                                                                 gAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGt                     <  1:5242102/149‑1 (MQ=255)
                                                                                                                                 gAAAACCTGCGCCTCTGGCAGCGTAACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGt                     <  2:4952605/149‑1 (MQ=255)
                                                                                                                                  aaaaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGtt                    <  2:1651653/149‑1 (MQ=255)
                                                                                                                                   aaaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                   <  2:5208774/149‑1 (MQ=255)
                                                                                                                                   aaaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                   <  2:4719369/149‑1 (MQ=255)
                                                                                                                                     aCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAtt                 <  1:1867607/149‑1 (MQ=255)
                                                                                                                                       cTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGc               >  1:2554889/1‑149 (MQ=255)
                                                                                                                                       cTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGc               <  1:1001357/149‑1 (MQ=255)
                                                                                                                                        tGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGct              >  1:250076/1‑149 (MQ=255)
                                                                                                                                         gcgcCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctc             <  1:3237181/149‑1 (MQ=255)
                                                                                                                                          cgcCTCTGTCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctct            <  2:3548439/149‑1 (MQ=255)
                                                                                                                                          cgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCGGCATGAAGCGGTTATTGctct            <  2:4817905/149‑1 (MQ=255)
                                                                                                                                          cgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctct            >  2:1262478/1‑149 (MQ=255)
                                                                                                                                            ccTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTg          >  1:1106166/1‑149 (MQ=255)
                                                                                                                                             ctctGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGc         <  2:2446440/149‑1 (MQ=255)
                                                                                                                                              tctGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTCAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGcc        >  2:3271303/1‑149 (MQ=255)
                                                                                                                                                tGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTc      <  2:3817620/149‑1 (MQ=255)
                                                                                                                                                    aGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAccc  <  1:1968652/149‑1 (MQ=255)
                                                                                                                                                    aGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAccc  <  2:349179/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCC  >  CP009273/1915407‑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: