Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,639,551 C→A 26.8% T269T (ACC→ACA 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,5510CA26.8% 725.6 / 237.6 323T269T (ACC→ACArlmJ23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (115/121);  new base A (53/34);  total (168/155)
Fisher's exact test for biased strand distribution p-value = 5.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

CTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAGC  >  CP009273/3639402‑3639700
                                                                                                                                                     |                                                                                                                                                     
   cAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCggg                                                                                                                                                    <  2:4421349/150‑1 (MQ=255)
   cAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCggg                                                                                                                                                    <  1:4784115/150‑1 (MQ=255)
    aaaTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGc                                                                                                                                                   >  2:2294897/1‑150 (MQ=255)
     aaTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCa                                                                                                                                                  >  2:4486840/1‑150 (MQ=255)
      aTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCAc                                                                                                                                                 <  1:142706/150‑1 (MQ=255)
         gAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGcc                                                                                                                                              >  2:3784498/1‑150 (MQ=255)
         gAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGcc                                                                                                                                              >  1:2974162/1‑150 (MQ=255)
           aCTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCAc                                                                                                                                            >  1:2154762/1‑150 (MQ=255)
              ggCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGt                                                                                                                                         >  2:4652502/1‑150 (MQ=255)
              ggCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGt                                                                                                                                         >  1:1199961/1‑150 (MQ=255)
                cGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTaa                                                                                                                                       >  2:1351025/1‑150 (MQ=255)
                cGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTaa                                                                                                                                       <  1:2193013/150‑1 (MQ=255)
                cGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTaa                                                                                                                                       >  2:1717431/1‑150 (MQ=255)
                cGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTaa                                                                                                                                       >  2:1718746/1‑150 (MQ=255)
                   tACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCt                                                                                                                                    <  2:1627047/150‑1 (MQ=255)
                   tACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCt                                                                                                                                    >  1:3495753/1‑150 (MQ=255)
                       gCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGAt                                                                                                                                >  2:1270658/1‑150 (MQ=255)
                          aGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGt                                                                                                                             <  2:534668/150‑1 (MQ=255)
                          aGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGt                                                                                                                             >  1:2502402/1‑150 (MQ=255)
                             cAGCGATCGCCGTGGCATGACCGCTTCCTGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGcc                                                                                                                          <  1:1032757/150‑1 (MQ=255)
                             cAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGcc                                                                                                                          <  2:1227955/150‑1 (MQ=255)
                              aGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCg                                                                                                                         >  2:1250430/1‑150 (MQ=255)
                                cGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGa                                                                                                                       <  2:3353999/150‑1 (MQ=255)
                                  aTCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGt                                                                                                                     >  1:3079830/1‑150 (MQ=255)
                                     gCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAAt                                                                                                                  >  2:3881011/1‑150 (MQ=255)
                                         tGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCa                                                                                                              >  2:4685544/1‑150 (MQ=255)
                                          ggCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAg                                                                                                             >  1:4503928/1‑150 (MQ=255)
                                                 ccGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGc                                                                                                      >  2:2473690/1‑150 (MQ=255)
                                                  cGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCt                                                                                                     >  1:2880651/1‑150 (MQ=255)
                                                   gCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTg                                                                                                    <  1:2032193/150‑1 (MQ=255)
                                                        cGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCAcc                                                                                               <  2:1880629/150‑1 (MQ=255)
                                                        cGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCAcc                                                                                               >  2:3967497/1‑150 (MQ=255)
                                                        cGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCAcc                                                                                               <  2:4583560/150‑1 (MQ=255)
                                                         ggCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCt                                                                                              >  1:737008/1‑150 (MQ=255)
                                                             tGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTa                                                                                          <  1:2815823/150‑1 (MQ=255)
                                                              gATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTAc                                                                                         >  1:3604852/1‑150 (MQ=255)
                                                               aTTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACg                                                                                        >  2:1436430/1‑150 (MQ=255)
                                                               aTTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACg                                                                                        >  2:2450251/1‑150 (MQ=255)
                                                               aTTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACg                                                                                        >  1:2614013/1‑150 (MQ=255)
                                                                   tgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGtctc                                                                                    <  2:92366/150‑1 (MQ=255)
                                                                   tgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGtctc                                                                                    >  1:4536956/1‑150 (MQ=255)
                                                                         aCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTa                                                                              >  2:2282922/1‑150 (MQ=255)
                                                                          cccGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTAc                                                                             <  1:848063/150‑1 (MQ=255)
                                                                           ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACa                                                                            <  1:2812094/150‑1 (MQ=255)
                                                                             gccgTGGAAACTGGAACTACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAAt                                                                          >  1:31709/1‑150 (MQ=255)
                                                                              ccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATc                                                                         <  1:1333053/150‑1 (MQ=255)
                                                                                  ggAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCgg                                                                     >  1:770963/1‑150 (MQ=255)
                                                                                       cTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATc                                                                <  1:4209537/150‑1 (MQ=255)
                                                                                       cTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATc                                                                >  2:1396554/1‑150 (MQ=255)
                                                                                        tGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCa                                                               <  1:2457817/150‑1 (MQ=255)
                                                                                         ggAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCaa                                                              >  2:1522603/1‑150 (MQ=255)
                                                                                           aacaacAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACg                                                            <  2:4180496/150‑1 (MQ=255)
                                                                                           aacaacAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACg                                                            <  1:4121981/150‑1 (MQ=255)
                                                                                              aacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATa                                                         >  1:4460781/1‑150 (MQ=255)
                                                                                                 aGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                                                      <  2:960850/150‑1 (MQ=255)
                                                                                                   aTGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGac                                                    <  2:2662025/150‑1 (MQ=255)
                                                                                                   aTGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGac                                                    >  1:923141/1‑150 (MQ=255)
                                                                                                     gAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGacac                                                  >  2:588931/1‑150 (MQ=255)
                                                                                                      aaTAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACt                                                 >  1:2240668/1‑150 (MQ=255)
                                                                                                        tAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACttt                                               >  1:4350471/1‑150 (MQ=255)
                                                                                                        tAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACttt                                               <  1:1899454/150‑1 (MQ=255)
                                                                                                        tAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACttt                                               >  1:4399023/1‑150 (MQ=255)
                                                                                                         aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTg                                              <  2:3844261/150‑1 (MQ=255)
                                                                                                         aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTg                                              <  2:4006220/150‑1 (MQ=255)
                                                                                                          aCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGt                                             >  1:3943134/1‑150 (MQ=255)
                                                                                                            gTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCa                                           >  1:169321/1‑150 (MQ=255)
                                                                                                                tgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGAc                                       <  2:450081/150‑1 (MQ=255)
                                                                                                                tgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGAc                                       <  2:46407/150‑1 (MQ=255)
                                                                                                                tgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGAc                                       >  2:2761545/1‑150 (MQ=255)
                                                                                                                tgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGAc                                       >  2:680128/1‑150 (MQ=255)
                                                                                                                 gcCGTGGCTGCACAGCAAACTGGGTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACt                                      <  2:2793768/150‑1 (MQ=255)
                                                                                                                  ccGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTa                                     >  1:1096684/1‑150 (MQ=255)
                                                                                                                  ccGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTa                                     >  2:1061524/1‑150 (MQ=255)
                                                                                                                  ccGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTa                                     >  2:1061465/1‑150 (MQ=255)
                                                                                                                   cGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaa                                    >  2:1095676/1‑150 (MQ=255)
                                                                                                                    gTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaaa                                   >  1:1550856/1‑150 (MQ=255)
                                                                                                                     tGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAac                                  >  1:658934/1‑150 (MQ=255)
                                                                                                                       gCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAacac                                <  1:477104/150‑1 (MQ=255)
                                                                                                                         tGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTa                              <  2:3648877/150‑1 (MQ=255)
                                                                                                                         tGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTa                              <  2:428856/150‑1 (MQ=255)
                                                                                                                           cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                            <  1:2640797/150‑1 (MQ=255)
                                                                                                                           cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                            >  1:1794427/1‑150 (MQ=255)
                                                                                                                             caGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAt                          >  1:3148532/1‑150 (MQ=255)
                                                                                                                               gCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTa                        <  2:3459492/150‑1 (MQ=255)
                                                                                                                                 aaaCTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACa                      <  2:2115274/150‑1 (MQ=255)
                                                                                                                                   aCTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATc                    <  1:3825544/150‑1 (MQ=255)
                                                                                                                                      ggTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCCTTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGcc                 >  2:831146/1‑150 (MQ=255)
                                                                                                                                       gTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCa                >  1:458109/1‑150 (MQ=255)
                                                                                                                                          ccGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcg             <  2:298738/150‑1 (MQ=255)
                                                                                                                                           cGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcgg            <  1:4652487/150‑1 (MQ=255)
                                                                                                                                           cGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcgg            >  2:746688/1‑150 (MQ=255)
                                                                                                                                           cGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcgg            >  2:3068578/1‑150 (MQ=255)
                                                                                                                                             gcaggcaCAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcg          >  2:294680/1‑150 (MQ=255)
                                                                                                                                              caggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcgg         <  1:3721169/150‑1 (MQ=255)
                                                                                                                                              caggcaCAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcgg         >  1:3991204/1‑150 (MQ=255)
                                                                                                                                               aggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcggc        >  1:1646017/1‑150 (MQ=255)
                                                                                                                                                ggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcggcg       >  2:3183376/1‑150 (MQ=255)
                                                                                                                                                ggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcggcg       >  2:4082505/1‑150 (MQ=255)
                                                                                                                                                   aCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCa    <  2:2154762/150‑1 (MQ=255)
                                                                                                                                                     cGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAgc  <  2:369057/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAGC  >  CP009273/3639402‑3639700

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: