Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,144 G→A 7.5% R2H (CGT→CAT)  rlmKL → 23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,389,1440GA 7.5% 151.3 / 2.3 54R2H (CGT→CAT) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base A (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

ATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCG  >  NC_002947/2389005‑2389292
                                                                                                                                            |                                                                                                                                                    
atgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                             <  4:307915/149‑1 (MQ=255)
   atCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCATTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                             >  1:84928/1‑147 (MQ=255)
   atCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCACGGGGCGACGCATCGCCCCgcga                                                                                                                                           <  3:311491/148‑4 (MQ=255)
    tCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagt                                                                                                                                         >  3:138386/1‑149 (MQ=255)
     cTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagt                                                                                                                                         >  6:102243/1‑148 (MQ=255)
       aaGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagtca                                                                                                                                       >  8:242494/1‑148 (MQ=255)
        aGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaaa                                                                                                                                     <  5:149193/149‑1 (MQ=255)
               cTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGccc                                                                                                                               <  8:111985/148‑1 (MQ=255)
                         ggCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcag                                                                                                                                          <  5:42692/127‑1 (MQ=255)
                         ggCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcag                                                                                                                                          >  6:42692/1‑127 (MQ=255)
                         ggCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTg                                                                                                                      <  4:138386/147‑1 (MQ=255)
                              ggCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCTAATAGCCCAAGGCGCTGGAATACt                                                                                                               <  2:84928/149‑1 (MQ=255)
                              ggCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                               <  4:96941/149‑1 (MQ=255)
                              ggCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                               <  7:242494/149‑1 (MQ=255)
                               gCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                               >  3:79409/1‑148 (MQ=255)
                                        ggTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGAc                                                                                                     <  5:102243/149‑1 (MQ=255)
                                         gTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGGCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGAcc                                                                                                    >  1:98644/1‑149 (MQ=255)
                                          tCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGc                                                                                                                                 <  4:18549/119‑1 (MQ=255)
                                            aGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGtt                                                                                                 <  6:173085/149‑1 (MQ=255)
                                                  cccGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAAc                                                                                           >  4:41232/1‑149 (MQ=255)
                                                   ccGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACt                                                                                          >  5:67045/1‑149 (MQ=255)
                                                   ccGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACt                                                                                          >  7:209906/1‑149 (MQ=255)
                                                    cGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                         >  6:289135/1‑149 (MQ=255)
                                                     ggC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACttt                                                                                        <  3:41232/149‑1 (MQ=255)
                                                     ggC‑TTTTCTGCAGCAAGTTGTGTGTCGCAGAAGGCCACCAGCACGACTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                         >  4:273411/1‑148 (MQ=255)
                                                                cagcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTg                                                                              >  2:10353/1‑149 (MQ=255)
                                                                   caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                                       <  2:83615/105‑1 (MQ=255)
                                                                   caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTcc                                                                           <  8:209906/149‑1 (MQ=255)
                                                                   caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGAAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                                       >  1:83615/1‑105 (MQ=255)
                                                                      gTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCa                                                                         >  5:80731/1‑148 (MQ=255)
                                                                       ttGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAt                                                                                                            <  1:284045/112‑1 (MQ=255)
                                                                       ttGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAt                                                                                                            >  2:284045/1‑112 (MQ=255)
                                                                              tCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                 >  6:36529/1‑148 (MQ=255)
                                                                                  aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCCCCTGTCCCAAAGGCCTTGAAAg                                                            >  6:14011/1‑149 (MQ=255)
                                                                                  aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCCTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACCCATGTCGGACCGTTTAGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAg                                                            >  7:230810/1‑149 (MQ=255)
                                                                                       gCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGc                                                       >  1:18115/1‑149 (MQ=255)
                                                                                            ccGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGcc                                                  <  5:289135/149‑1 (MQ=255)
                                                                                                aCGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCg                                                 >  4:288589/1‑146 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                     <  2:98644/148‑1 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                     <  6:67045/148‑1 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                     <  5:14011/148‑1 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAg                                      <  8:21838/147‑1 (MQ=255)
                                                                                                               gtgtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCg                               <  6:80731/149‑1 (MQ=255)
                                                                                                                  tACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                     <  3:134864/140‑1 (MQ=255)
                                                                                                                  tACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                     >  4:134864/1‑140 (MQ=255)
                                                                                                                      ggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGa                        <  3:273411/149‑1 (MQ=255)
                                                                                                                            cATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGt                  <  3:288589/149‑1 (MQ=255)
                                                                                                                              tgacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAagag                                            >  6:208320/4‑121 (MQ=255)
                                                                                                                              tgacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAagag                                            <  5:208320/118‑1 (MQ=255)
                                                                                                                               gacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcg               >  2:181998/3‑149 (MQ=255)
                                                                                                                                gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc              >  3:285611/1‑149 (MQ=255)
                                                                                                                                gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc              >  6:99058/1‑149 (MQ=255)
                                                                                                                                gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc              <  2:18115/149‑1 (MQ=255)
                                                                                                                                 cgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGcgcg             <  4:79409/149‑1 (MQ=255)
                                                                                                                                            gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTc   >  7:104964/1‑148 (MQ=255)
                                                                                                                                            gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTc   <  8:104964/148‑1 (MQ=255)
                                                                                                                                            gTCGCCCCGACAGTCAAATAGCCCAAGGCGCTGGAAAACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCCAAGAGGCCAAGCGCCACGGCCTTGACGAGGTGCGCGACCACACCTCg  >  1:69672/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
ATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGC‑TTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCG  >  NC_002947/2389005‑2389292

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