Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 77,500 A→T 10.8% L15H (CTC→CAC)  rsmB ← ribosomal RNA small subunit methyltransferase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294777,5000AT10.8% 226.5 / 15.7 85L15H (CTC→CAC) rsmBribosomal RNA small subunit methyltransferase B
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base T (6/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

GCTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGG  >  NC_002947/77351‑77648
                                                                                                                                                      |                                                                                                                                                     
gCTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGa                                                                                                                                                       <  2:751298/150‑1 (MQ=255)
        agcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                              <  2:1636650/150‑1 (MQ=255)
            gcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCa                                                                                                                                          <  1:102333/150‑1 (MQ=255)
            gcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCa                                                                                                                                          >  1:158795/1‑150 (MQ=255)
                tgCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAgcgc                                                                                                                                      >  2:1750679/1‑150 (MQ=255)
                   ggCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACg                                                                                                                                   <  2:158795/150‑1 (MQ=255)
                    gCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACgg                                                                                                                                  <  1:1240974/150‑1 (MQ=255)
                     ccAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGGCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACggg                                                                                                                                 >  2:742256/1‑150 (MQ=255)
                      cagcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggc                                                                                                                                >  1:279410/1‑150 (MQ=255)
                      cagcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGGCAGGCCGCGGTCGCGCTCATCTACCTTGTCCAGTTGCGCCGGCACCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggc                                                                                                                                >  2:1267893/1‑150 (MQ=255)
                            gTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                              >  2:589362/1‑130 (MQ=255)
                            gTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                              <  1:589362/130‑1 (MQ=255)
                               gAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAg                                                                                                                       <  1:984014/150‑1 (MQ=255)
                               gAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTCTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAg                                                                                                                       <  1:1877650/150‑1 (MQ=255)
                                aGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAgg                                                                                                                      <  2:1460717/150‑1 (MQ=255)
                                   cGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGTGGCCAGCGCACGGGCGGCGGCGAGGCGt                                                                                                                   <  1:1632004/150‑1 (MQ=255)
                                    gTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTg                                                                                                                  >  2:1800650/1‑150 (MQ=255)
                                         tGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACCGCGGCCCGCGCACGGGCGGCGGCGGGGCGTGGGTTc                                                                                                             >  1:950057/1‑150 (MQ=255)
                                         tGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTc                                                                                                             >  1:1070288/1‑150 (MQ=255)
                                             aGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTg                                                                                                         <  1:1546293/150‑1 (MQ=255)
                                             aGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCACCTCCTTGCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTg                                                                                                         <  1:19273/150‑1 (MQ=255)
                                              gCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTgg                                                                                                        >  1:1740172/1‑150 (MQ=255)
                                                gttcGGTACCGCACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGcc                                                                                                      >  1:666996/4‑150 (MQ=255)
                                                    ggTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGc                                                                                                  >  2:1327579/1‑150 (MQ=255)
                                                      tACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCAc                                                                                                >  1:1549536/1‑150 (MQ=255)
                                                          gAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCAGCTTg                                                                                            >  1:755308/1‑150 (MQ=255)
                                                            aCGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGcc                                                                                          <  1:1327579/150‑1 (MQ=255)
                                                               ccAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCcggc                                                                                       <  2:1549536/150‑1 (MQ=255)
                                                                      cTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                >  2:1692472/1‑150 (MQ=255)
                                                                      cTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                >  2:1692328/1‑150 (MQ=255)
                                                                          gTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                              >  1:1351186/1‑84 (MQ=255)
                                                                          gTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                              <  2:1351186/84‑1 (MQ=255)
                                                                           tCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGCGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTtctc                                                                           <  1:808461/150‑1 (MQ=255)
                                                                               gCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGTCATACATCCATCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCgcgg                                                                       >  2:543444/1‑150 (MQ=255)
                                                                                ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCgcggc                                                                      <  1:156136/150‑1 (MQ=255)
                                                                                ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCgcggc                                                                      <  1:755933/150‑1 (MQ=255)
                                                                                     gTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTg                                                                 <  2:484080/150‑1 (MQ=255)
                                                                                      tCGCGTTCATCTACCTTGTCAAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGt                                                                >  2:691321/1‑150 (MQ=255)
                                                                                       cgcgTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGtt                                                               <  2:279410/150‑1 (MQ=255)
                                                                                          gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGaa                                                            >  1:1511045/1‑150 (MQ=255)
                                                                                          gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGaa                                                            <  1:194047/150‑1 (MQ=255)
                                                                                             cATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAg                                                         <  2:1740172/150‑1 (MQ=255)
                                                                                              aTCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAgg                                                        >  2:109301/1‑150 (MQ=255)
                                                                                                   ccTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACt                                                   >  1:545460/1‑150 (MQ=255)
                                                                                                   ccTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGGGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACt                                                   >  1:76284/1‑150 (MQ=255)
                                                                                                     ttGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGa                                                 >  2:703080/1‑150 (MQ=255)
                                                                                                     ttGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTAACtaa                                                 >  1:789657/1‑148 (MQ=255)
                                                                                                     ttGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTACCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGa                                                 >  1:903909/1‑150 (MQ=255)
                                                                                                       tGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                                <  2:1070288/150‑1 (MQ=255)
                                                                                                       tGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                                <  1:1762181/150‑1 (MQ=255)
                                                                                                            aGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGcaa                                                                                                                                          >  1:1046581/1‑53 (MQ=37)
                                                                                                            aGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGcaa                                                                                                                                          <  2:1046581/55‑3 (MQ=37)
                                                                                                            aGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGt                                                                >  2:201665/1‑129 (MQ=255)
                                                                                                            aGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGt                                                                <  1:201665/129‑1 (MQ=255)
                                                                                                            aGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCa                                           >  2:1573236/1‑150 (MQ=255)
                                                                                                                  gcCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCt                                     <  1:1579063/150‑1 (MQ=255)
                                                                                                                   ccGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGTCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCtt                                    <  2:580714/150‑1 (MQ=255)
                                                                                                                    cGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTg                                   <  2:1684221/150‑1 (MQ=255)
                                                                                                                    cGGCAGCGAGCTGTTCAGCGAAGCCTTGACGATGAGGGCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTg                                   <  2:152552/150‑1 (MQ=255)
                                                                                                                            aGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                >  2:1688457/1‑97 (MQ=255)
                                                                                                                            aGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                <  1:1688457/97‑1 (MQ=255)
                                                                                                                              cTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCa                         <  2:47071/150‑1 (MQ=255)
                                                                                                                               tGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAg                        >  1:1385648/1‑150 (MQ=255)
                                                                                                                                gTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGt                       <  2:625012/150‑1 (MQ=255)
                                                                                                                                 ttCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGtt                      <  1:278601/150‑1 (MQ=255)
                                                                                                                                    aGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCa                   >  1:1272018/1‑150 (MQ=255)
                                                                                                                                    aGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCa                   <  2:250153/150‑1 (MQ=255)
                                                                                                                                     gCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAg                  >  1:1315942/1‑150 (MQ=255)
                                                                                                                                           gctggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGAGCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGt            >  2:150833/5‑150 (MQ=255)
                                                                                                                                           gctggCCGCTGTGCTCAGCGGCAAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAg                            >  1:113831/5‑134 (MQ=255)
                                                                                                                                           gctggCCGCTGTGCTCAGCGGCAAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAg                            <  2:113831/130‑1 (MQ=255)
                                                                                                                                           ccTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGtt                      <  1:1237607/140‑1 (MQ=255)
                                                                                                                                           ccTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGtt                      >  2:1237607/1‑140 (MQ=255)
                                                                                                                                           ccTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGt            >  1:963894/1‑150 (MQ=255)
                                                                                                                                            cTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTc           >  1:63469/1‑150 (MQ=255)
                                                                                                                                             ttGCCGCTGGTCACAGCGGCGAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCGCGCGGCGGCTGTTGAACAGGTCACTGGAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa          <  2:950057/150‑1 (MQ=255)
                                                                                                                                             ttGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa          >  1:1414368/1‑150 (MQ=255)
                                                                                                                                             ttGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa          >  1:25861/1‑150 (MQ=255)
                                                                                                                                             tggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGCAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa          >  2:937879/3‑150 (MQ=255)
                                                                                                                                             tggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGtac          >  1:27783/3‑148 (MQ=255)
                                                                                                                                             tggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATAGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGGTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa          >  2:853862/3‑150 (MQ=255)
                                                                                                                                              tGCCGCTGAGCACAGCGGCCAGCGAACGGGCGGCGGCGAGGCGTGCGTTCATTGGCCCAGCACCTTGCCGGAACCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAg         <  1:125064/150‑1 (MQ=255)
                                                                                                                                               gCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAgg        <  2:63469/150‑1 (MQ=255)
                                                                                                                                                  gCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGCCAATCTCTCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTc     >  2:65043/1‑150 (MQ=255)
                                                                                                                                                   cTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCa    <  2:27783/150‑1 (MQ=255)
                                                                                                                                                     gAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCGATGCCTTGCCGCCAGGCAGTGGCAGGCGGGTCAGGCTCAgg  <  2:545460/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCTTTTGCAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGG  >  NC_002947/77351‑77648

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