Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,639,351 C→A 9.8% Q203K (CAA→AAA)  rlmJ → 23S rRNA m(6)A2030 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,639,3510CA 9.8% 649.0 / 7.8 206Q203K (CAA→AAA) rlmJ23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (126/53);  new base A (2/25);  total (128/78)
Fisher's exact test for biased strand distribution p-value = 3.11e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CGGTTTCCCGCCGTGGTTTAATCCTTATCGACCCGCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAAC  >  CP009273/3639205‑3639498
                                                                                                                                                  |                                                                                                                                                   
  gtgtCCCGCCGTGGTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCaaa                                                                                                                                                  <  1:742692/145‑1 (MQ=255)
     tCCCGCCGTGGTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTa                                                                                                                                               <  2:2766096/148‑1 (MQ=255)
     tCCCGCCGTGGTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTa                                                                                                                                               <  1:853383/148‑1 (MQ=255)
             tGGTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATg                                                                                                                                       <  1:3430950/148‑1 (MQ=255)
             tGGTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATg                                                                                                                                       <  1:2143694/148‑1 (MQ=255)
             tGGTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATg                                                                                                                                       >  1:3217200/1‑148 (MQ=255)
              ggTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGa                                                                                                                                      >  2:2805472/1‑148 (MQ=255)
              ggTTTAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGa                                                                                                                                      <  2:1070812/148‑1 (MQ=255)
                 ttAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATcc                                                                                                                                   <  2:2567015/148‑1 (MQ=255)
                 ttAATCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATcc                                                                                                                                   >  1:3258137/1‑148 (MQ=255)
                    aTCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACg                                                                                                                                <  2:2898561/148‑1 (MQ=255)
                    aTCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACg                                                                                                                                >  2:2031518/1‑148 (MQ=255)
                     tCCTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGa                                                                                                                               <  1:3377015/148‑1 (MQ=255)
                      ccTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGAt                                                                                                                              <  2:3299733/148‑1 (MQ=255)
                      ccTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGAt                                                                                                                              <  2:1558155/148‑1 (MQ=255)
                      ccTTATCGACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGAt                                                                                                                              <  1:1692201/148‑1 (MQ=255)
                        ttATCGACCCTCCGTAAGAAATGAAAAATGACAATCAAGCGGTGGTCAGCGGGATAGCATAAGGTTACAAACGATTCGCAAATGTTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGAATGATCCACGATCt                                                                                                                            >  2:1132224/1‑148 (MQ=255)
                             gACCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAg                                                                                                                       <  1:1372246/148‑1 (MQ=255)
                              aCCCTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGTTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGAAAATTAAGCGCATGATCCACGATCTGGAAga                                                                                                                      >  1:1510943/1‑147 (MQ=255)
                               cccTCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCg                                                                                                                     >  2:1043752/1‑148 (MQ=255)
                                cctccgTATTAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGa                                                                                                                    <  2:1126792/145‑1 (MQ=255)
                                cctccgTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGa                                                                                                                    >  2:1960889/4‑148 (MQ=255)
                                cctccgTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGAAAATTAAGCGCATGATCCACGATCTGGAAGCGa                                                                                                                    <  2:2976543/145‑1 (MQ=255)
                                  tccgTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGAcc                                                                                                                  >  1:3502604/2‑148 (MQ=255)
                                         gaaatgaaaACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTc                                                                                                           >  2:1761778/1‑148 (MQ=255)
                                          aaatgaaaACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCg                                                                                                          >  2:3032180/1‑148 (MQ=255)
                                           aatgaaaACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGc                                                                                                         >  2:3472966/1‑148 (MQ=255)
                                            atgaaaACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCa                                                                                                        >  1:1889154/1‑148 (MQ=255)
                                            atgaaaACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCa                                                                                                        >  2:1316885/1‑148 (MQ=255)
                                              gaaaACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCaaa                                                                                                      >  1:1745483/1‑148 (MQ=255)
                                                   cTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCt                                                                                                 >  2:2347053/1‑148 (MQ=255)
                                                    tGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCGCTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTg                                                                                                >  1:1008463/1‑148 (MQ=255)
                                                       cTATCACGGGGGGGTCAGCGGTATCGCAGAAGGTTACAACCGTTGCGCCACGGGTTTTTCCGCACGGTGGTATCCGGTGGTGCTGGGTAAGCAAATTACGCGCATGATCCACGATGGGGACGCGACCGGTATTCGCAAAATTCTGCaa                                                                                             <  2:681453/148‑1 (MQ=255)
                                                         aTCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAAt                                                                                           >  1:1913143/1‑148 (MQ=255)
                                                          tCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAAtt                                                                                          >  2:3083852/1‑148 (MQ=255)
                                                          tCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAAtt                                                                                          >  1:1247979/1‑148 (MQ=255)
                                                             aGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGaa                                                                                       >  2:1879513/1‑148 (MQ=255)
                                                              gCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAAc                                                                                      >  2:3251678/1‑148 (MQ=255)
                                                               cGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACt                                                                                     >  1:916815/1‑148 (MQ=255)
                                                               cGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACt                                                                                     >  2:2208434/1‑148 (MQ=255)
                                                               cGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACt                                                                                     >  1:3393549/1‑148 (MQ=255)
                                                               cGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACt                                                                                     >  2:3016102/1‑148 (MQ=255)
                                                                 gtggtCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTgg                                                                                   >  1:2503526/1‑148 (MQ=255)
                                                                  tggtCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGc                                                                                  <  2:887012/148‑1 (MQ=255)
                                                                  tggtCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGc                                                                                  >  1:2697115/1‑148 (MQ=255)
                                                                  tggtCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGc                                                                                  >  1:3137221/1‑148 (MQ=255)
                                                                    gtCAGGGGTTTAGCTGATGTTTACAAACGTTTCGCCACTGGTATTTCGGCCGTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATTATCCACGATCTGGACGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCgg                                                                                <  2:570715/148‑1 (MQ=255)
                                                                    gtCAGGGGTATAGCAGACGGTTACAACCGTTTCGCCACTGGTATTTCCGCACTGTGGTATCGGGTGGTGCTGCTTCAGAAAATTAAGCGCATGATCCACGATCTGCAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCgg                                                                                <  1:561894/148‑1 (MQ=255)
                                                                    gtCAGCGGTATAGCAGTCGGTTACAACCGTTTCGCCACTGGTATTTCCGCAGTGTGGTATCGGTTGGTGGTGCGTCAGAAAATTAAGCGCTTGATCCCCGATCTGTAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCgg                                                                                <  2:3281195/148‑1 (MQ=255)
                                                                     tCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGt                                                                               >  1:2116066/1‑148 (MQ=255)
                                                                      cAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTa                                                                              >  2:2179501/1‑148 (MQ=255)
                                                                         cGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTg                                                                           >  2:2055876/1‑148 (MQ=255)
                                                                         cGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTg                                                                           >  2:2055046/1‑148 (MQ=255)
                                                                                gCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACa                                                                    >  2:3045039/1‑148 (MQ=255)
                                                                                 cAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAg                                                                   >  2:520912/1‑148 (MQ=255)
                                                                                   gAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCg                                                                 >  2:2233392/1‑148 (MQ=255)
                                                                                    atgtttACAAAGTTTTCGCCACTGTTATTTCCGCAGTGTGGTATCGGGTGGTGCTGCGTCAGAAAATTAAGCGCATTATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGa                                                                <  1:2031518/144‑1 (MQ=255)
                                                                                    aaGGTTACAACCGTTTCGCCACTGGTTTTTCCGCACGGTGTTATCCGGGGGGGCGGCGTCAGCAAATTAAGCGCATTATCCACGAGCTGGACGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGa                                                                <  2:3367984/148‑1 (MQ=255)
                                                                                    aaGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGa                                                                >  2:376702/1‑148 (MQ=255)
                                                                                     aGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGAt                                                               >  1:2703256/1‑148 (MQ=255)
                                                                                       gTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCg                                                             >  2:1213138/1‑148 (MQ=255)
                                                                                             accGTTTCGCCAGGGGTTTTTACGCAGTTTGGTTTCCGGTGGTGCTGCGTCAGAATATTAAGCGCATGATCCCCGATCTGGACGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTgg                                                       <  2:1503608/146‑1 (MQ=255)
                                                                                               ggTTTCGCCAGTGTTATTTCCGAAGTGTGTTATCCGGTGGTGCTGCTACACCAAATTACCCGCCTGACCCACGATCTGCAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCa                                                     <  1:439440/147‑1 (MQ=255)
                                                                                                  ttCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGa                                                  >  2:1104111/1‑148 (MQ=255)
                                                                                                  ttCGCCACTGGTATTTACGCACTGTGGTATCCGGGGGTGCTGCGTCAGAAAATTAAGCGCATGATCCACGATCTGGACGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGa                                                  <  2:3271427/148‑1 (MQ=255)
                                                                                                  ttCGCAACTGGTATTACCTCACGGGGTTATCCGGTGGTGCGGCGTCAGCAAATAAAGCGCATTATCCCCGATCTGCAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGa                                                  <  2:1842968/148‑1 (MQ=255)
                                                                                                  ttCGCAACTGGTATTACCGCCCTGTGTTATCCGGTGGTGCGGCGAAAGCAAATTAAGCGCATGATCCCCCCTCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGa                                                  <  2:1639281/148‑1 (MQ=255)
                                                                                                   tCGCCAGTGTTATTTCCGCTGTGTGGTATCCGGTGGTGCTGCGACAGAAAATTAAGCGCATTATCCACGATCTGCAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGAc                                                 <  2:1486435/148‑1 (MQ=255)
                                                                                                   tCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGAc                                                 >  1:2252251/1‑148 (MQ=255)
                                                                                                          tGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTcc                                          >  1:451300/1‑148 (MQ=255)
                                                                                                          tGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTcc                                          >  2:287367/1‑148 (MQ=255)
                                                                                                            gTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCgg                                        >  1:2838548/1‑148 (MQ=255)
                                                                                                             tATTTACGCAGTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGc                                       <  2:2200925/148‑1 (MQ=255)
                                                                                                              aTTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCa                                      >  2:853261/1‑148 (MQ=255)
                                                                                                                ttACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATg                                    >  2:758193/1‑148 (MQ=255)
                                                                                                                ttACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATg                                    >  2:2480555/1‑148 (MQ=255)
                                                                                                                ttACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATg                                    >  1:1593224/1‑148 (MQ=255)
                                                                                                                ttACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATg                                    >  1:663075/1‑148 (MQ=255)
                                                                                                                 tccGCAGTGTGGTAGCGGGTGGTGGGGCGTCAGAAATTTAACCGCTTTAGACACGTTGTGCAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGa                                   <  1:3089517/146‑1 (MQ=255)
                                                                                                                   cGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGAtt                                 >  2:2449154/1‑148 (MQ=255)
                                                                                                                    gCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATtg                                >  1:1391342/1‑148 (MQ=255)
                                                                                                                        tgtgGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGAt                            >  1:3462942/1‑148 (MQ=255)
                                                                                                                        tgtgGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGAt                            >  2:622815/1‑148 (MQ=255)
                                                                                                                        tgtgGTAGCCGGTGGTGCTGCGTCAGAAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGAt                            <  2:1704309/148‑1 (MQ=255)
                                                                                                                         gtgGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGAtt                           >  1:51912/1‑148 (MQ=255)
                                                                                                                           ggTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTaa                         >  1:3456609/1‑148 (MQ=255)
                                                                                                                           ggTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTaa                         >  1:2406871/1‑148 (MQ=255)
                                                                                                                            tttgccGGTGGTGCTGCGTCAGAAAATTACGCGCATGATCCACTATCGGCAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAAc                        <  1:1288742/144‑1 (MQ=255)
                                                                                                                               tCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACccg                     >  2:3451891/1‑148 (MQ=255)
                                                                                                                                  ggtggtGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACccgccg                  <  2:2527124/148‑1 (MQ=255)
                                                                                                                                  ggtgggGCGGCGTCCCCAAATTACGCGAAGGAGCCACTTTGTGCACGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACccgccg                  <  1:106705/148‑1 (MQ=255)
                                                                                                                                   gtggtgGGGCGTCAGCAAATTACGCGCATGATCCACGATCGGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGt                 <  2:1599720/148‑1 (MQ=255)
                                                                                                                                    tggtgCTGCTTCAGCAAATTAAGCGAATGAACCCCGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTg                <  1:812524/148‑1 (MQ=255)
                                                                                                                                    tggtgCTGCTTCAGAAAATTACGCGAATGATCCACGATGTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTg                <  1:2570882/148‑1 (MQ=255)
                                                                                                                                       tgctgcGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGaa             >  1:2821655/1‑148 (MQ=255)
                                                                                                                                       tgctgcGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGACAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGAGTGTGATTAACCCGCCGTGGaa             >  1:1965823/1‑148 (MQ=255)
                                                                                                                                         ctgcGTCAGAAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAAc           <  1:2197864/148‑1 (MQ=255)
                                                                                                                                               cAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGaac     >  2:2662066/1‑148 (MQ=255)
                                                                                                                                                  cAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGaacaac  >  2:1911273/1‑148 (MQ=255)
                                                                                                                                                  cAAATTAAGCGAATGATCCACGATCTGTAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGaacaac  <  2:1015013/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGTTTCCCGCCGTGGTTTAATCCTTATCGACCCGCCGTATGAAATGAAAACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAAC  >  CP009273/3639205‑3639498

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