Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 77,503 A→T 10.8% V14E (GTG→GAG)  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,5030AT10.8% 238.5 / 15.1 85V14E (GTG→GAG) 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 = 9.75e-01

AGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCTT  >  NC_002947/77359‑77652
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                                  <  1:589362/130‑1 (MQ=255)
                    gTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                                  >  2:589362/1‑130 (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:1692328/1‑150 (MQ=255)
                                                              cTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                    >  2:1692472/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:755933/150‑1 (MQ=255)
                                                                        ccGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCgcggc                                                                          <  1:156136/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                                                                                                                                              <  2:1046581/55‑3 (MQ=37)
                                                                                                    aGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGcaa                                                                                                                                              >  1:1046581/1‑53 (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                                                                                    <  1:1688457/97‑1 (MQ=255)
                                                                                                                    aGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                    >  2:1688457/1‑97 (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:25861/1‑150 (MQ=255)
                                                                                                                                     ttGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa              >  1:1414368/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)
                                                                                                                                                 cacaGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCt   >  2:637720/1‑150 (MQ=255)
                                                                                                                                                  acaGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                             >  2:1056461/1‑107 (MQ=255)
                                                                                                                                                  acaGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                             <  1:1056461/107‑1 (MQ=255)
                                                                                                                                                  acaGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTGCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCtt  <  1:1875069/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCTT  >  NC_002947/77359‑77652

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: