Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,718,111 G→C 5.9% L256F (TTG→TTC rsmJ → ribosomal RNA small subunit methyltransferase J

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,718,1110GC 5.9% 617.9 / 21.3 185L256F (TTG→TTCrsmJribosomal RNA small subunit methyltransferase J
Reads supporting (aligned to +/- strand):  ref base G (84/90);  new base C (5/6);  total (89/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 192 total aligned reads displayed.

CCGGCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCA  >  NC_002947/1717962‑1718260
                                                                                                                                                     |                                                                                                                                                     
ccGGCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTg                                                                                                                                                       <  1:2019599/150‑1 (MQ=255)
ccGGCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTg                                                                                                                                                       >  1:2051629/1‑150 (MQ=255)
  ggCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGaa                                                                                                                                                     <  2:3964008/150‑1 (MQ=255)
    cgcTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAgg                                                                                                                                                   <  1:349447/150‑1 (MQ=255)
     gctgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGc                                                                                                                                                  <  1:4970978/150‑1 (MQ=255)
     gctgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGc                                                                                                                                                  <  2:2051329/150‑1 (MQ=255)
       tgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCt                                                                                                                                                >  1:4275548/1‑150 (MQ=255)
       tgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCt                                                                                                                                                <  1:5065139/150‑1 (MQ=255)
         ctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttc                                                                                                                                              <  2:3251594/150‑3 (MQ=255)
          tCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGAc                                                                                                                                             >  1:3366411/1‑150 (MQ=255)
          tCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGAc                                                                                                                                             >  1:361168/1‑150 (MQ=255)
          tCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGAc                                                                                                                                             <  2:1083769/150‑1 (MQ=255)
           cGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttctt                                                                                                                                            <  2:1238131/150‑5 (MQ=255)
            gAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCt                                                                                                                                           >  2:3111228/1‑150 (MQ=255)
            gAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCt                                                                                                                                           >  2:787296/1‑150 (MQ=255)
             aaGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTg                                                                                                                                          >  1:2222291/1‑150 (MQ=255)
              aGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGc                                                                                                                                         <  1:3560623/150‑1 (MQ=255)
              aGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcgg                                                                                                                                         <  2:1029438/150‑8 (MQ=255)
              aGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcgg                                                                                                                                         <  1:452947/150‑8 (MQ=255)
                ccGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGccc                                                                                                                                       <  2:601460/150‑1 (MQ=255)
                  gCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCt                                                                                                                                     <  1:5036798/150‑1 (MQ=255)
                            tGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCgg                                                                                                                           <  2:2538893/150‑1 (MQ=255)
                              gccagccaCCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCt                                                                                                                         >  1:3460647/1‑150 (MQ=255)
                                cagccaCCGGGTGGTAGTGAAGCGCCCTCGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAg                                                                                                                       >  1:3909917/1‑150 (MQ=255)
                                  gccaCCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCa                                                                                                                     >  1:5207538/1‑150 (MQ=255)
                                    caCCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACg                                                                                                                   >  2:3544027/1‑150 (MQ=255)
                                         ggtggtAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGt                                                                                                              <  1:4828210/150‑1 (MQ=255)
                                         ggtggtAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGt                                                                                                              >  1:3100680/1‑150 (MQ=255)
                                          gtggtAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTa                                                                                                             <  2:3621732/150‑1 (MQ=255)
                                          gtggtAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTa                                                                                                             <  2:2407644/150‑1 (MQ=255)
                                          gtggtAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTa                                                                                                             >  1:4885713/1‑150 (MQ=255)
                                            ggtAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAgg                                                                                                           <  1:4424260/150‑1 (MQ=255)
                                              tAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAg                                                                                                         <  1:5138706/150‑1 (MQ=255)
                                                   aaGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCt                                                                                                    >  1:3368314/1‑150 (MQ=255)
                                                   aaGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCATATCTGGCCGGCTAGCACGATTGTAGGAGCGGCtt                                                                                                   <  2:968358/150‑1 (MQ=255)
                                                     gcgcCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCt                                                                                                                                                <  2:5041858/104‑1 (MQ=255)
                                                       gcCCACGCAAGGCACCGATCATCACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGc                                                                                               >  1:290337/1‑150 (MQ=255)
                                                          cACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCcg                                                                                             >  2:4425470/1‑150 (MQ=255)
                                                           aCGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCcgc                                                                                            <  2:2254700/150‑1 (MQ=255)
                                                           aCGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCcgc                                                                                            <  1:1074292/150‑1 (MQ=255)
                                                           aCGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCcgc                                                                                            <  1:2294424/150‑1 (MQ=255)
                                                              cAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGaa                                                                                         <  2:4029953/150‑1 (MQ=255)
                                                              cAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGaa                                                                                         <  1:2677974/150‑1 (MQ=255)
                                                                aGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCcgc                                                                                            <  2:2078146/145‑1 (MQ=255)
                                                                 ggCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAacac                                                                                      <  2:1624765/150‑1 (MQ=255)
                                                                    aCCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGa                                                                                   <  1:2184894/150‑1 (MQ=255)
                                                                     ccGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGAc                                                                                  <  1:1618531/150‑1 (MQ=255)
                                                                          catcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGc                                                                             <  2:1047974/150‑1 (MQ=255)
                                                                           atcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGcc                                                                            >  2:3014514/1‑150 (MQ=255)
                                                                            tcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCg                                                                           >  2:492491/1‑150 (MQ=255)
                                                                            tcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCg                                                                           >  2:1845170/1‑150 (MQ=255)
                                                                            tcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCg                                                                           <  1:3164644/150‑1 (MQ=255)
                                                                            tcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCg                                                                           >  1:5426249/1‑150 (MQ=255)
                                                                                cGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGc                                                                       <  1:3757454/150‑1 (MQ=255)
                                                                                 gggCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGcc                                                                      >  1:1484600/1‑150 (MQ=255)
                                                                                  ggCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCa                                                                     <  2:1697737/150‑1 (MQ=255)
                                                                                    ccGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATc                                                                   >  2:4796722/1‑150 (MQ=255)
                                                                                      gAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGc                                                                                                                                                  >  2:4113571/1‑69 (MQ=255)
                                                                                      gAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGc                                                                                                                                                  <  1:4113571/69‑1 (MQ=255)
                                                                                          ccGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGc                                                             <  2:3397592/150‑1 (MQ=255)
                                                                                            gAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGAtt                                                                                                                <  2:3214698/97‑1 (MQ=255)
                                                                                            gAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGAtt                                                                                                                >  1:3214698/1‑97 (MQ=255)
                                                                                            gAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGcgt                                                           >  2:3280014/1‑150 (MQ=255)
                                                                                               ccATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGcgtcgt                                                        <  1:2441331/150‑1 (MQ=255)
                                                                                                cATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCgtcgtc                                                       <  1:4588821/150‑1 (MQ=255)
                                                                                                cATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCgtcgtc                                                       <  1:3699033/150‑1 (MQ=255)
                                                                                                cATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCgtcgtc                                                       <  1:1847892/150‑1 (MQ=255)
                                                                                                 aTAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGtct                                                      >  2:3525959/1‑150 (MQ=255)
                                                                                                  tAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGtctc                                                     >  2:1224433/1‑150 (MQ=255)
                                                                                                    gCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCt                                                   <  1:702705/150‑1 (MQ=255)
                                                                                                    gCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCt                                                   <  2:3496322/150‑1 (MQ=255)
                                                                                                       ttGAGGGTAAGTCGAGCCGGCATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTgcg                                                >  2:1368357/1‑150 (MQ=255)
                                                                                                          aGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGt                                             <  2:2364771/150‑1 (MQ=255)
                                                                                                             gTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTa                                                                                                                        >  2:2798176/1‑72 (MQ=255)
                                                                                                             gTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTa                                                                                                                        <  1:2798176/72‑1 (MQ=255)
                                                                                                                aGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAAcc                                       >  2:5201334/1‑150 (MQ=255)
                                                                                                                  tCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTg                                     >  2:2316584/1‑150 (MQ=255)
                                                                                                                          gTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTa                               >  1:1241199/1‑148 (MQ=255)
                                                                                                                              gACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGCAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGa                         <  1:1891107/150‑1 (MQ=255)
                                                                                                                                 aTCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCt                      >  2:1471065/1‑150 (MQ=255)
                                                                                                                                  tCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTg                     >  2:3796653/1‑150 (MQ=255)
                                                                                                                                    tACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGtt                   >  1:3309534/1‑150 (MQ=255)
                                                                                                                                     aCCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTa                  >  1:1166699/1‑150 (MQ=255)
                                                                                                                                      cccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTAt                 >  2:2159388/1‑150 (MQ=255)
                                                                                                                                       ccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACACTGACCTGTTAtt                >  2:4786660/1‑150 (MQ=255)
                                                                                                                                         gaagaaGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTc              >  1:128962/1‑150 (MQ=255)
                                                                                                                                          ggtcagGCCTTCAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCt             >  1:1989772/5‑150 (MQ=255)
                                                                                                                                          ggtcagGCCTTCAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCt             >  1:886398/5‑150 (MQ=255)
                                                                                                                                           gtcagGCCTTCAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTc              <  2:4506635/145‑1 (MQ=255)
                                                                                                                                              agGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCa         >  2:2386289/1‑150 (MQ=255)
                                                                                                                                              agGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCa         >  2:1541189/1‑150 (MQ=255)
                                                                                                                                                 ccTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCg                                                                           >  1:4698542/1‑81 (MQ=255)
                                                                                                                                                 ccTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagta      <  2:1988660/150‑1 (MQ=255)
                                                                                                                                                 ccTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagta      <  1:4707237/150‑1 (MQ=255)
                                                                                                                                                 ccTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagta      >  2:229294/1‑150 (MQ=255)
                                                                                                                                                 ccTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGAGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACATGTAATTTCTGCCagta      >  2:3096610/1‑150 (MQ=255)
                                                                                                                                                  cTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagtag     <  2:4678470/150‑1 (MQ=255)
                                                                                                                                                  cTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagtag     >  1:2967012/1‑150 (MQ=255)
                                                                                                                                                   ttGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagtagt    >  1:2829728/1‑150 (MQ=255)
                                                                                                                                                     gAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCa  >  1:2049151/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCA  >  NC_002947/1717962‑1718260

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: