Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,646,1960AG48.8% 18.5 / 123.7 84H271R (CAC→CGC) rlmJ23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base G (20/21);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

AACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAGCGGC  >  NC_000913/3646052‑3646343
                                                                                                                                                |                                                                                                                                                   
aaCTGGGGGTACGGCCAGACAGCGATCGCCGTGGCATGCCCGCTTCCGGGATTATTGTTATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCAc                                                                                                                                               <  2:555398/151‑1 (MQ=255)
aaCTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCAc                                                                                                                                               >  2:202769/1‑151 (MQ=255)
       ggTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGc                                                                                                                                        >  2:243632/1‑151 (MQ=255)
       ggTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGc                                                                                                                                        >  2:247571/1‑151 (MQ=255)
           cTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTg                                                                                                                                      <  1:117338/149‑1 (MQ=255)
             gCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGAt                                                                                                                                   <  2:280313/150‑1 (MQ=255)
             gCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATc                                                                                                                                  >  2:34387/1‑151 (MQ=255)
              ccAGACGGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATc                                                                                                                                  <  2:21950/150‑1 (MQ=255)
               cAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGt                                                                                                                                >  1:19644/1‑151 (MQ=255)
                aGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGt                                                                                                                                <  1:501472/150‑1 (MQ=255)
                aGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTg                                                                                                                               >  1:454943/1‑151 (MQ=255)
                   cAGCGATCGCCGTGGCATGCCCGCTCCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCg                                                                                                                            <  2:608560/151‑1 (MQ=255)
                      cGATCGCCGTGGCATGACCGCTCCCGGAATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGa                                                                                                                          <  2:3006/150‑1 (MQ=255)
                       gATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAg                                                                                                                         >  1:245376/1‑150 (MQ=255)
                        aTCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGt                                                                                                                        <  1:490753/150‑1 (MQ=255)
                            ccGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTg                                                                                                                   <  1:36623/151‑1 (MQ=255)
                            ccGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTg                                                                                                                   <  1:592036/151‑1 (MQ=255)
                            ccGGGGCATGACCGCTTCCGGCAGGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTaa                                                                                                                      <  1:589029/148‑1 (MQ=255)
                              gTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGc                                                                                                                  >  2:283998/1‑150 (MQ=255)
                                ggCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCa                                                                                                                 <  2:223886/149‑1 (MQ=255)
                                  cATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCa                                                                                                             <  1:50702/151‑1 (MQ=255)
                                     gcccGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATATCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCAtt                                                                                                           <  2:473879/148‑1 (MQ=255)
                                     gACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTg                                                                                                          >  2:595544/1‑151 (MQ=255)
                                     gACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATAAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTg                                                                                                          >  2:462714/1‑151 (MQ=255)
                                      atcGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGc                                                                                                         >  1:560406/3‑151 (MQ=255)
                                          cTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGc                                                                                                     >  1:230935/1‑151 (MQ=255)
                                           ttCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGc                                                                                                     >  1:547488/1‑150 (MQ=255)
                                           ttCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGc                                                                                                     >  1:260761/1‑150 (MQ=255)
                                            tCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCAGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCAc                                                                                                   >  2:186608/1‑151 (MQ=255)
                                             ccGGCATGTTTGTGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCa                                                                                                    <  2:633005/149‑1 (MQ=255)
                                             ccGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCAc                                                                                                   <  2:22096/150‑1 (MQ=255)
                                               ggCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCt                                                                                                 <  2:133428/150‑1 (MQ=255)
                                                  aTGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTAtt                                                                                              >  1:431624/1‑150 (MQ=255)
                                                   tGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTAc                                                                                            >  1:233014/1‑151 (MQ=255)
                                                      ttggtttaaaCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTa                                                                                                <  2:214636/136‑1 (MQ=255)
                                                        gtgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGtctc                                                                                       <  1:689/151‑1 (MQ=255)
                                                        gtgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGtctc                                                                                       >  1:460083/1‑151 (MQ=255)
                                                           aTTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTcgcg                                                                                    >  1:615080/1‑151 (MQ=255)
                                                              accccGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTa                                                                                 <  2:200194/149‑1 (MQ=255)
                                                              aaCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCt                                                                                  <  2:608118/150‑1 (MQ=255)
                                                                cccGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTAc                                                                                >  1:148897/1‑150 (MQ=255)
                                                                 ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACaa                                                                              <  2:143503/151‑1 (MQ=255)
                                                                 ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACa                                                                               >  1:450111/1‑150 (MQ=255)
                                                                   gccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATc                                                                            >  1:322488/1‑151 (MQ=255)
                                                                   gccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATc                                                                            >  1:7651/1‑151 (MQ=255)
                                                                    ccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATcg                                                                           >  1:484307/1‑151 (MQ=255)
                                                                    ccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATcg                                                                           <  2:528240/151‑1 (MQ=255)
                                                                       tGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCgg                                                                        >  1:513227/1‑151 (MQ=255)
                                                                         gAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGt                                                                       >  2:456658/1‑150 (MQ=255)
                                                                          aaaCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTa                                                                      >  2:145187/1‑150 (MQ=255)
                                                                           aaCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAAt                                                                    >  1:199773/1‑151 (MQ=255)
                                                                            aCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATc                                                                   >  2:621548/1‑151 (MQ=255)
                                                                            aCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATc                                                                   <  1:214838/151‑1 (MQ=255)
                                                                                gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACg                                                               <  2:230930/151‑1 (MQ=255)
                                                                                          tGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGacac                                                     >  2:501480/1‑151 (MQ=255)
                                                                                            aaTAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACtt                                                   <  1:600737/151‑1 (MQ=255)
                                                                                                aCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTc                                               >  2:455154/1‑151 (MQ=255)
                                                                                                 cGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTc                                               >  2:262971/1‑150 (MQ=255)
                                                                                                  gTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCAt                                             <  2:317653/151‑1 (MQ=255)
                                                                                                  gTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCAt                                             >  1:151011/1‑151 (MQ=255)
                                                                                                   tgctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                            >  1:556132/1‑151 (MQ=255)
                                                                                                   tgctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                            >  1:470466/1‑151 (MQ=255)
                                                                                                    gctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                            <  1:89590/150‑1 (MQ=255)
                                                                                                    gctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                            <  2:452673/150‑1 (MQ=255)
                                                                                                         cGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaaa                                      <  2:142039/151‑1 (MQ=255)
                                                                                                          gTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAac                                     >  2:233430/1‑151 (MQ=255)
                                                                                                            ggCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAaca                                    <  1:97748/150‑1 (MQ=255)
                                                                                                            ggCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAacac                                   <  2:366155/151‑1 (MQ=255)
                                                                                                             gCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACt                                  <  2:222652/151‑1 (MQ=255)
                                                                                                               tGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTAt                                >  1:5689/1‑151 (MQ=255)
                                                                                                                gCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTAt                                <  2:119666/150‑1 (MQ=255)
                                                                                                                gCACAGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                               >  2:585663/1‑151 (MQ=255)
                                                                                                                gCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                               >  2:599117/1‑151 (MQ=255)
                                                                                                                 cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                               >  2:608549/1‑150 (MQ=255)
                                                                                                                  acaGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGa                              <  1:42713/150‑1 (MQ=255)
                                                                                                                    aGCAAACTGGTTCCGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAtt                            >  2:430808/1‑150 (MQ=255)
                                                                                                                    aGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTa                           <  1:591486/151‑1 (MQ=255)
                                                                                                                         aCTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCg                      >  2:493806/1‑151 (MQ=255)
                                                                                                                                ccGGCAGGCACCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcgg               <  2:547477/151‑1 (MQ=255)
                                                                                                                                  ggcaggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcg             >  1:531291/1‑151 (MQ=255)
                                                                                                                                     aggcaCCGGGCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcggcg          <  2:257351/151‑1 (MQ=255)
                                                                                                                                         aCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAg      <  1:537691/151‑1 (MQ=255)
                                                                                                                                            gggCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAgcgg   >  2:597514/1‑151 (MQ=255)
                                                                                                                                             ggCGCGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAgcggc  <  2:376193/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAGCGGC  >  NC_000913/3646052‑3646343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: