Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,915,555 A→C 22.6% 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,5550AC22.6% 614.8 / 236.7 319Q359P (CAG→CCG) rsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (121/125);  new base C (33/39);  total (154/165)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 320 total aligned reads displayed.

AAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTC  >  CP009273/1915408‑1915699
                                                                                                                                                   |                                                                                                                                                
aaaCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGc                                                                                                                                                <  1:3003979/150‑1 (MQ=255)
  aCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGt                                                                                                                                              >  1:4242324/1‑150 (MQ=255)
   cTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTg                                                                                                                                             <  1:1245490/150‑1 (MQ=255)
   cTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTg                                                                                                                                             >  1:1538068/1‑150 (MQ=255)
   cTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTg                                                                                                                                             >  1:159558/1‑150 (MQ=255)
     cAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGAc                                                                                                                                           >  2:4316697/1‑150 (MQ=255)
      aGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACa                                                                                                                                          <  1:3484052/150‑1 (MQ=255)
      aGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACa                                                                                                                                          >  1:4194854/1‑150 (MQ=255)
           aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGaa                                                                                                                                     >  1:812676/1‑150 (MQ=255)
           aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGaa                                                                                                                                     >  1:3868757/1‑150 (MQ=255)
           aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGaa                                                                                                                                     <  2:1938527/150‑1 (MQ=255)
           aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGaa                                                                                                                                     <  2:5073333/150‑1 (MQ=255)
           aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGCCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACCAAGaa                                                                                                                                     >  2:3019431/1‑150 (MQ=255)
             tCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACt                                                                                                                                   <  2:655191/150‑1 (MQ=255)
               ccGCCTTACCCGCCCCCAAATACACAGTCGGTAATTGTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtgt                                                                                                                                 <  2:5024874/150‑1 (MQ=255)
               ccGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtgt                                                                                                                                 <  1:1253706/150‑1 (MQ=255)
                 gCCTTACCCGCCCCCACATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTgg                                                                                                                               <  1:4453960/150‑1 (MQ=255)
                    ttACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgttg                                                                                                                            <  2:1746219/150‑1 (MQ=255)
                     tACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgttgt                                                                                                                           >  2:2236541/1‑150 (MQ=255)
                      aCCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgttgtt                                                                                                                          <  1:4751081/150‑1 (MQ=255)
                         cGCCCCCAAAGACAATGTCGGTAATGTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTccc                                                                                                                       <  2:4034994/150‑1 (MQ=255)
                            ccccAAATACAAAGTCGGTAATTTTCCGGTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGt                                                                                                                    <  2:2426408/150‑1 (MQ=255)
                            caccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGAGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACACAGAACTGTGGTTGTTCCCGGt                                                                                                                    >  2:536376/3‑150 (MQ=255)
                             cccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTg                                                                                                                   >  2:2716186/1‑150 (MQ=255)
                             cccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTg                                                                                                                   <  2:286190/150‑1 (MQ=255)
                              ccAAATACAAACTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTgg                                                                                                                  <  2:4823456/150‑1 (MQ=255)
                                 aaTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCa                                                                                                               <  1:2367226/150‑1 (MQ=255)
                                  aTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAt                                                                                                              <  2:4474965/150‑1 (MQ=255)
                                  aTACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAt                                                                                                              <  2:2766921/150‑1 (MQ=255)
                                   tACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAtt                                                                                                             <  2:1615593/150‑1 (MQ=255)
                                      acaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGaa                                                                                                          <  1:946247/148‑1 (MQ=255)
                                      aaaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGaa                                                                                                          <  1:1768765/150‑1 (MQ=255)
                                       aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAg                                                                                                         >  1:589735/1‑150 (MQ=255)
                                        aGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGc                                                                                                        >  2:3676560/1‑150 (MQ=255)
                                           cGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCt                                                                                                     <  2:266073/150‑1 (MQ=255)
                                             gTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGa                                                                                                   >  2:3632892/1‑150 (MQ=255)
                                             gTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGa                                                                                                   <  2:3319858/150‑1 (MQ=255)
                                               aaTTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATc                                                                                                 <  1:4113065/150‑1 (MQ=255)
                                                      cGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAg                                                                                          >  1:3159891/1‑150 (MQ=255)
                                                      cGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAg                                                                                          <  1:1657500/150‑1 (MQ=255)
                                                       gTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGt                                                                                         <  1:3333673/150‑1 (MQ=255)
                                                         tCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTcc                                                                                       >  2:4225448/1‑150 (MQ=255)
                                                          cAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCg                                                                                      >  2:2190446/1‑150 (MQ=255)
                                                           aGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGa                                                                                     >  2:405925/1‑150 (MQ=255)
                                                            gCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAt                                                                                    <  2:3411278/150‑1 (MQ=255)
                                                            gCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAt                                                                                    <  2:786832/150‑1 (MQ=255)
                                                             cccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAtt                                                                                   <  2:2830978/150‑1 (MQ=255)
                                                              ccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAACACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAttt                                                                                  <  2:433222/150‑1 (MQ=255)
                                                               cGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAtttt                                                                                 <  1:2633887/150‑1 (MQ=255)
                                                                ggTGAAAGATCGCGAAGCTGGACACATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtc                                                                                <  2:560383/150‑1 (MQ=255)
                                                                ggTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtc                                                                                >  1:2431940/1‑150 (MQ=255)
                                                                 gTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtct                                                                               <  1:4491356/150‑1 (MQ=255)
                                                                   gAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCg                                                                             <  1:5055630/150‑1 (MQ=255)
                                                                    aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCgg                                                                            >  1:1362554/1‑150 (MQ=255)
                                                                    aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCgg                                                                            >  2:127457/1‑150 (MQ=255)
                                                                    aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCgg                                                                            <  1:3368308/150‑1 (MQ=255)
                                                                     aaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                           <  2:4640357/150‑1 (MQ=255)
                                                                     aaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                           >  2:1449191/1‑150 (MQ=255)
                                                                     aaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                           >  1:2036912/1‑150 (MQ=255)
                                                                      aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGtt                                                                          >  1:3209683/1‑150 (MQ=255)
                                                                         tCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTggg                                                                       <  2:919913/150‑1 (MQ=255)
                                                                          cgagAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTgggg                                                                      <  2:5234981/147‑1 (MQ=255)
                                                                             gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAtt                                                                   <  1:1301177/150‑1 (MQ=255)
                                                                              aaGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTa                                                                  >  1:5107685/1‑150 (MQ=255)
                                                                               aGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTaa                                                                 <  2:1195882/150‑1 (MQ=255)
                                                                                      cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGcc                                                          <  1:2896277/150‑1 (MQ=255)
                                                                                      cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGcc                                                          <  2:3848302/150‑1 (MQ=255)
                                                                                      cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGcc                                                          <  1:206622/150‑1 (MQ=255)
                                                                                          ttCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                      <  1:775755/150‑1 (MQ=255)
                                                                                           tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGAc                                                     <  1:3624817/150‑1 (MQ=255)
                                                                                             gTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGc                                                   >  1:3985527/1‑150 (MQ=255)
                                                                                                aGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGcaca                                                <  1:1045695/150‑1 (MQ=255)
                                                                                                   cggcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCacaaca                                             <  1:405025/150‑1 (MQ=255)
                                                                                                    ggcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCAcaacaa                                            <  1:879666/150‑1 (MQ=255)
                                                                                                     gcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACaaa                                           >  1:3945246/1‑150 (MQ=255)
                                                                                                     gcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACaaa                                           >  1:2951697/1‑150 (MQ=255)
                                                                                                      cTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAg                                          >  1:4394774/1‑150 (MQ=255)
                                                                                                        aCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGt                                        >  2:2207862/1‑150 (MQ=255)
                                                                                                            gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                    >  1:652477/1‑150 (MQ=255)
                                                                                                            gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                    <  1:5238855/150‑1 (MQ=255)
                                                                                                            gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                    <  1:3410531/150‑1 (MQ=255)
                                                                                                             tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCg                                   >  2:5141859/1‑150 (MQ=255)
                                                                                                             tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCg                                   >  2:36043/1‑150 (MQ=255)
                                                                                                             tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCg                                   <  1:2282246/150‑1 (MQ=255)
                                                                                                               ttGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCt                                 >  1:3895547/1‑150 (MQ=255)
                                                                                                                tGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTg                                >  2:2487977/1‑150 (MQ=255)
                                                                                                                       aaCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCAt                         >  2:3085511/1‑150 (MQ=255)
                                                                                                                           gggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                     <  2:5179763/150‑1 (MQ=255)
                                                                                                                               tGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGt                 <  1:1935909/150‑1 (MQ=255)
                                                                                                                                 aaaaCCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa               >  2:2925564/1‑150 (MQ=255)
                                                                                                                                 aaaaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGGTATCGCTGGCAGCATGAAGCGGTTa               >  2:1175885/1‑150 (MQ=255)
                                                                                                                                  aaaCCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAt              <  1:128057/150‑1 (MQ=255)
                                                                                                                                   aaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGGTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAtt             <  1:4553263/150‑1 (MQ=255)
                                                                                                                                     ccTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGc           >  1:4907396/1‑150 (MQ=255)
                                                                                                                                       tGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctc         >  2:2547242/1‑150 (MQ=255)
                                                                                                                                         cgcCTCTGGCAGCGTGACAACGCACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTACAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCtt       <  2:2890978/150‑1 (MQ=255)
                                                                                                                                         cgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCtt       <  1:2099730/150‑1 (MQ=255)
                                                                                                                                         cgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCtt       <  2:4242324/150‑1 (MQ=255)
                                                                                                                                         cgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCtt       >  2:1897318/1‑150 (MQ=255)
                                                                                                                                           ccTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGc     >  1:4016459/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTC  >  CP009273/1915408‑1915699

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: