Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 76,886 A→G 7.1% F220L (TTC→CTC)  rsmB ← ribosomal RNA small subunit methyltransferase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294776,8860AG 7.1% 160.8 / 2.2 56F220L (TTC→CTC) rsmBribosomal RNA small subunit methyltransferase B
Reads supporting (aligned to +/- strand):  ref base A (30/22);  new base G (2/2);  total (32/24)
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.97e-01

GCCAGGCCAGCCTCGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAATCA  >  NC_002947/76738‑77032
                                                                                                                                                    |                                                                                                                                                  
gccaggccagCCTCGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGaa                                                                                                                                                    >  4:13529/1‑149 (MQ=255)
 ccaggccagCCTCGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGaa                                                                                                                                                    <  8:58866/148‑1 (MQ=255)
  caggccagCCTCGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGc                                                                                                                                                  >  5:297107/1‑149 (MQ=255)
       cagCCTCGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGc                                                                                                                                             >  8:232886/1‑149 (MQ=255)
             cGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCAGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCcg                                                                                                                                       >  8:183012/1‑149 (MQ=255)
              ggCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCcg                                                                                                                                       <  3:13529/148‑1 (MQ=255)
                ccTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCgcgc                                                                                                                                    <  4:158047/149‑1 (MQ=255)
                                 gTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGc                                                                                                                   >  2:51976/1‑149 (MQ=255)
                                      gccgccCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCa                                                                                                                 <  4:101333/146‑1 (MQ=255)
                                       ccgccCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCAc                                                                                                             >  3:88237/1‑149 (MQ=255)
                                        cgccCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCAc                                                                                                             >  2:19211/1‑148 (MQ=255)
                                         gccCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGa                                                                                                           >  3:202158/1‑149 (MQ=255)
                                          cccGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGAt                                                                                                          <  1:113668/149‑1 (MQ=255)
                                           ccGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATg                                                                                                         <  6:297107/149‑1 (MQ=255)
                                             ggTGAGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTGCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCCCGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGc                                                                                                        >  5:174752/1‑148 (MQ=255)
                                                    atcaagcaTAGAGAACTCGTTTGCCGGGGGCCAGTTCCAGCATGTCTGCAGATATCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCTGCAAGCACGATGCCATCCCgg                                                                                                <  4:167819/147‑1 (MQ=255)
                                                     gcaagcaTCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGc                                                                                               <  4:88237/149‑1 (MQ=255)
                                                      caagcaTCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCCCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCATGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCt                                                                                              >  7:144059/1‑149 (MQ=255)
                                                       aagcaTCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCt                                                                                              <  7:183012/148‑1 (MQ=255)
                                                       aagcaTCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTg                                                                                             >  5:255091/1‑149 (MQ=255)
                                                            tCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAGGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTg                                                                                        >  5:298606/1‑149 (MQ=255)
                                                            tCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAGGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTg                                                                                        >  5:294388/1‑149 (MQ=255)
                                                            tCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTg                                                                                        <  1:51976/149‑1 (MQ=255)
                                                              gAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGc                                                                                      <  3:235899/149‑1 (MQ=255)
                                                                  aCGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATgc                                                                                  >  5:212219/1‑149 (MQ=255)
                                                                    gcgTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCTGCCAGCACGATGCCATCCCGGCTGTACTGGCATgcgc                                                                                >  8:260602/1‑149 (MQ=255)
                                                                      gTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCTTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTg                                                                              <  8:144059/149‑1 (MQ=255)
                                                                                 gCCAGTTACAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGcc                                                                   >  5:90037/1‑149 (MQ=255)
                                                                                      ttCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCtt                                                              >  3:133906/1‑149 (MQ=255)
                                                                                        ccATCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGACGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTc                                                             <  7:279535/148‑1 (MQ=255)
                                                                                         cagcagGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCgg                                                           >  3:35434/1‑149 (MQ=255)
                                                                                           gcagGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGGTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGcc                                                         <  6:90037/149‑1 (MQ=255)
                                                                                               gTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCcagca                                                     <  4:133906/149‑1 (MQ=255)
                                                                                                     aGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTg                                                  >  2:292166/1‑146 (MQ=255)
                                                                                                     aGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTg                                                  <  1:292166/146‑1 (MQ=255)
                                                                                                          gCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAg                                          >  3:226512/1‑149 (MQ=255)
                                                                                                                        tCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATg                            <  1:19211/149‑1 (MQ=255)
                                                                                                                        tCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATg                            >  2:287521/1‑149 (MQ=255)
                                                                                                                          cTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGCTGCCCTCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGt                          >  5:218386/1‑149 (MQ=255)
                                                                                                                          cTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGt                          <  5:19730/149‑1 (MQ=255)
                                                                                                                                gCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCa                                               <  6:168907/122‑1 (MQ=255)
                                                                                                                                gCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCa                                               >  5:168907/1‑122 (MQ=255)
                                                                                                                                gCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGAcggc                    <  6:255091/149‑1 (MQ=255)
                                                                                                                                  tCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGAcggcgg                  >  1:267367/1‑149 (MQ=255)
                                                                                                                                    aCCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGcc                                                                   <  6:96150/98‑1 (MQ=255)
                                                                                                                                    aCCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGcc                                                                   >  5:96150/1‑98 (MQ=255)
                                                                                                                                    aCCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGAcggcgg                  <  7:260602/147‑1 (MQ=255)
                                                                                                                                    aCCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGtt                >  4:134416/1‑149 (MQ=255)
                                                                                                                                      ccAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGa              >  3:5832/1‑149 (MQ=255)
                                                                                                                                       cAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGAc             >  2:231107/1‑149 (MQ=255)
                                                                                                                                              cGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCa        <  2:267367/147‑1 (MQ=255)
                                                                                                                                               ggCGAGGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGaa     <  6:298606/149‑1 (MQ=255)
                                                                                                                                               ggCGAGGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGaa     <  6:294388/149‑1 (MQ=255)
                                                                                                                                                gCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAat    >  2:134175/1‑149 (MQ=255)
                                                                                                                                                gCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAat    >  6:286009/1‑149 (MQ=255)
                                                                                                                                                gCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAat    >  1:299204/1‑149 (MQ=255)
                                                                                                                                                  gAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAatca  >  1:247489/1‑149 (MQ=255)
                                                                                                                                                  gAAGCCTGGCAGGCCGCGCACGTCGCAGGCNTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAatca  >  1:27221/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCAGGCCAGCCTCGGCCTCCAGCAGGTGGCAGGTCTTGCCGCCCGGTGCGCAGCAAGCATCGAGCACGCGTTGGCCGGGGGCCAGTTCCAGCAGGTCGGCAGATAGCTGCGCAGCTTCGTCCTGCACGCTCACCCAGCCTTCGGCGAAGCCTGGCAGGCCGCGCACGTCGCAGGCTTCGGCCAGCACGATGCCATCCCGGCTGTACTGGCATGCGCTGGCGCCGATGCCGGCTTCGGCCAGCAGTGCCAGGTAGGCATCGCGGCTATGGTGACGGCGGTTGACCCGCAGAATCA  >  NC_002947/76738‑77032

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