Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,915,555 A→C 15.9% 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,5550AC15.9% 898.1 / 180.4 343Q359P (CAG→CCG) rsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (131/157);  new base C (33/22);  total (164/179)
Fisher's exact test for biased strand distribution p-value = 5.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 350 total aligned reads displayed.

AAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCC  >  CP009273/1915407‑1915703
                                                                                                                                                    |                                                                                                                                                    
aaaaCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACGGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCa                                                                                                                                                      <  1:3874013/149‑1 (MQ=255)
aaaaCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCa                                                                                                                                                      >  1:2877122/1‑149 (MQ=255)
  aaCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGc                                                                                                                                                    >  2:2102003/1‑149 (MQ=255)
      cAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGa                                                                                                                                                <  2:2420107/149‑1 (MQ=255)
           gATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAg                                                                                                                                           >  1:913183/1‑149 (MQ=255)
            aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGa                                                                                                                                          <  2:1643363/149‑1 (MQ=255)
            aTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGa                                                                                                                                          >  1:1661226/1‑149 (MQ=255)
                 cGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtgt                                                                                                                                     <  1:396010/149‑1 (MQ=255)
                  gCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtgtg                                                                                                                                    >  2:4597143/1‑149 (MQ=255)
                  gCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACtgtg                                                                                                                                    >  1:4811662/1‑149 (MQ=255)
                   ccTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTgg                                                                                                                                   <  2:4998550/149‑1 (MQ=255)
                    cTTACCCGCCCCCCAATACACAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAAGTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGgt                                                                                                                                  <  2:4842453/149‑1 (MQ=255)
                         ccGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTCCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTc                                                                                                                             >  2:4454707/1‑149 (MQ=255)
                             ccccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCgg                                                                                                                         >  1:4476051/1‑149 (MQ=255)
                              cccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGt                                                                                                                        >  1:5172645/1‑149 (MQ=255)
                               ccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTg                                                                                                                       <  2:29201/149‑1 (MQ=255)
                               ccAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTg                                                                                                                       >  1:3313365/1‑149 (MQ=255)
                                    tACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAt                                                                                                                  >  2:3565074/1‑149 (MQ=255)
                                    tACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACATGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAt                                                                                                                  >  2:1164087/1‑149 (MQ=255)
                                     aCAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCAtt                                                                                                                 <  1:1513500/149‑1 (MQ=255)
                                      cAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTg                                                                                                                >  1:728806/1‑149 (MQ=255)
                                       acaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGa                                                                                                               <  1:3954918/147‑1 (MQ=255)
                                       aaaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGa                                                                                                               >  2:405912/1‑149 (MQ=255)
                                        aaGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGaa                                                                                                              >  1:313535/1‑149 (MQ=255)
                                                 aTTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATc                                                                                                     <  1:3787356/149‑1 (MQ=255)
                                                     tCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTa                                                                                                 <  1:4592509/149‑1 (MQ=255)
                                                      ccGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTaa                                                                                                >  1:21431/1‑149 (MQ=255)
                                                       cGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTaaa                                                                                               <  2:2091342/149‑1 (MQ=255)
                                                          tCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTc                                                                                            >  1:3771610/1‑149 (MQ=255)
                                                            aGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCg                                                                                          <  2:3725717/149‑1 (MQ=255)
                                                             gCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGa                                                                                         <  1:2080552/149‑1 (MQ=255)
                                                              cccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAt                                                                                        <  2:4989483/149‑1 (MQ=255)
                                                               ccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAtt                                                                                       <  2:1578292/149‑1 (MQ=255)
                                                               ccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAtt                                                                                       <  1:1813004/149‑1 (MQ=255)
                                                               ccGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAtt                                                                                       >  2:1103320/1‑149 (MQ=255)
                                                                cGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGAttt                                                                                      <  1:1298478/149‑1 (MQ=255)
                                                                  gTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGAGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtc                                                                                    >  2:2239230/1‑149 (MQ=255)
                                                                   tGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTtct                                                                                   >  2:3013115/1‑149 (MQ=255)
                                                                     aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCg                                                                                 >  1:22742/1‑149 (MQ=255)
                                                                     aaaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCg                                                                                 >  1:2626619/1‑149 (MQ=255)
                                                                      aaGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCgg                                                                                <  2:3189866/149‑1 (MQ=255)
                                                                       aGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGt                                                                               >  2:4644315/1‑149 (MQ=255)
                                                                        gATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGtt                                                                              >  2:2980188/1‑149 (MQ=255)
                                                                          tCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTgg                                                                            >  1:433571/1‑149 (MQ=255)
                                                                             cgAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGa                                                                         <  1:1062897/149‑1 (MQ=255)
                                                                              gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAt                                                                        >  1:4995361/1‑149 (MQ=255)
                                                                              gAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGAt                                                                        >  2:794455/1‑149 (MQ=255)
                                                                                 gCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTaa                                                                     <  1:1237703/149‑1 (MQ=255)
                                                                                   ttgACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAAc                                                                   <  2:5070311/147‑1 (MQ=255)
                                                                                    ggACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACt                                                                  <  1:2102267/149‑1 (MQ=255)
                                                                                    ggACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACt                                                                  <  1:4768320/149‑1 (MQ=255)
                                                                                    ggACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACt                                                                  >  2:1936825/1‑149 (MQ=255)
                                                                                    ggACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACt                                                                  <  1:2524177/149‑1 (MQ=255)
                                                                                     gACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACtt                                                                 <  2:1788173/149‑1 (MQ=255)
                                                                                     gACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACtt                                                                 <  2:3849796/149‑1 (MQ=255)
                                                                                       cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAATCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGc                                                               >  2:1360494/1‑149 (MQ=255)
                                                                                       cAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGc                                                               <  2:1879085/149‑1 (MQ=255)
                                                                                           ttCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgag                                                           >  2:2274010/1‑149 (MQ=255)
                                                                                            tCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCgaga                                                          <  1:3424337/149‑1 (MQ=255)
                                                                                              gTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACg                                                        >  1:851740/1‑149 (MQ=255)
                                                                                               tCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGc                                                       >  1:2353993/1‑149 (MQ=255)
                                                                                               tCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGc                                                       >  1:1457898/1‑149 (MQ=255)
                                                                                                 aGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGcac                                                     <  2:5140757/149‑1 (MQ=255)
                                                                                                 aGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGcac                                                     <  2:4135346/149‑1 (MQ=255)
                                                                                                   gcggcTACAGGTGTTGGTTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCacaa                                                   <  2:873700/149‑1 (MQ=255)
                                                                                                   gcggcTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCacaa                                                   <  2:2232167/149‑1 (MQ=255)
                                                                                                       cTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACaaa                                               >  1:4088598/1‑149 (MQ=255)
                                                                                                        tACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGGTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAg                                              >  1:3655234/1‑149 (MQ=255)
                                                                                                         aCAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAgg                                             >  1:322419/1‑149 (MQ=255)
                                                                                                          cAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGt                                            >  1:1480970/1‑149 (MQ=255)
                                                                                                           aGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGtt                                           <  1:2559262/149‑1 (MQ=255)
                                                                                                           aGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGtt                                           >  1:4326913/1‑149 (MQ=255)
                                                                                                              tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                        <  2:2656277/149‑1 (MQ=255)
                                                                                                              tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                        <  2:1229825/149‑1 (MQ=255)
                                                                                                              tgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATc                                        <  2:420274/149‑1 (MQ=255)
                                                                                                               gtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCg                                       >  1:3721575/1‑149 (MQ=255)
                                                                                                                ttGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGc                                      >  2:4374834/1‑149 (MQ=255)
                                                                                                                ttGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGc                                      <  2:633497/149‑1 (MQ=255)
                                                                                                                        aaCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGgcagca                              >  1:4150509/1‑149 (MQ=255)
                                                                                                                        aaCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCAGAGACGCACAACAAAGGTTATCGCTGgcagca                              >  1:1021554/1‑149 (MQ=255)
                                                                                                                         aCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCAt                             <  2:2955600/149‑1 (MQ=255)
                                                                                                                         aCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCAt                             <  2:4811662/149‑1 (MQ=255)
                                                                                                                             ggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                         <  2:1328500/149‑1 (MQ=255)
                                                                                                                              gATGAAAACCTGCGCCTCTGGCCGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGc                        <  1:2966512/149‑1 (MQ=255)
                                                                                                                                 gAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGt                     <  1:2661048/149‑1 (MQ=255)
                                                                                                                                  aaaaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGtt                    >  2:3264924/1‑149 (MQ=255)
                                                                                                                                   aaaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                   <  2:569962/149‑1 (MQ=255)
                                                                                                                                    aaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAt                  <  2:3428121/149‑1 (MQ=255)
                                                                                                                                    aaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAt                  >  1:602179/1‑149 (MQ=255)
                                                                                                                                    aaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAt                  >  2:2669501/1‑149 (MQ=255)
                                                                                                                                    aaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAt                  >  1:3045716/1‑149 (MQ=255)
                                                                                                                                      ccTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTg                >  1:300731/1‑149 (MQ=255)
                                                                                                                                        tGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGct              >  2:103668/1‑149 (MQ=255)
                                                                                                                                        tGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGct              >  2:1804448/1‑149 (MQ=255)
                                                                                                                                         gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGctc             <  1:2789574/149‑1 (MQ=255)
                                                                                                                                            ccTCTGGCAGAGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTg          <  1:676647/149‑1 (MQ=255)
                                                                                                                                             ctctGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGc         >  2:3166071/1‑149 (MQ=255)
                                                                                                                                              tctGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGcc        >  1:4698010/1‑149 (MQ=255)
                                                                                                                                                  gCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAc    >  2:4756699/1‑149 (MQ=255)
                                                                                                                                                    aGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAccc  <  2:72904/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: