Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,940,965 T→C L93L (TTA→TTG rlmM ← 23S rRNA 2'‑O‑ribose C2498 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,940,9650TC100.0% 220.7 / NA 63L93L (TTA→TTGrlmM23S rRNA 2'‑O‑ribose C2498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/35);  total (28/35)

TTCGCCAGCACCCCGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGTAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCGCCATCATCAGG  >  NC_000913/2940828‑2941099
                                                                                                                                         |                                                                                                                                      
ttCGCCAGCACCCCGGCATCGCGCAAGGCAGCGCGTAGCGTAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCaa                                                                                                                                      >  1:1026780/1‑140 (MQ=255)
ttCGCCAGCACCCCGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACTGCCCTTCTCTACTACGCCCTGCaa                                                                                                                                      >  1:171046/1‑140 (MQ=255)
ttCGCCAGCACCCCGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCaa                                                                                                                                      >  1:685265/1‑140 (MQ=255)
       gCACCCCGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCg                                                                                                                               <  1:1112663/140‑1 (MQ=255)
           cccGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACaa                                                                                                                           <  1:369591/140‑1 (MQ=255)
            ccGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAAt                                                                                                                          <  1:968018/140‑1 (MQ=255)
              ggCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATgg                                                                                                                        <  1:140497/140‑1 (MQ=255)
               gCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATggg                                                                                                                       <  1:357634/140‑1 (MQ=255)
                        aaGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACg                                                                                                              <  2:391241/140‑1 (MQ=255)
                        aaGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACg                                                                                                              <  1:343794/140‑1 (MQ=255)
                                gcgTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTg                                                                                                      >  1:302386/1‑140 (MQ=255)
                                        gAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCaaa                                                                                              <  2:1131867/140‑1 (MQ=255)
                                           cAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAAtgc                                                                                           <  2:903207/140‑1 (MQ=255)
                                              tAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAAtgctgc                                                                                        >  2:1141849/1‑140 (MQ=255)
                                              tAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAAtgctgc                                                                                        >  1:364058/1‑140 (MQ=255)
                                                        gCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACGCGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTccc                                                                              <  2:797556/140‑1 (MQ=255)
                                                        gCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAAtgc                                                                                           >  2:37821/1‑127 (MQ=255)
                                                        gCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAAtgc                                                                                           <  1:37821/127‑1 (MQ=255)
                                                        gCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTccc                                                                              <  2:1074982/140‑1 (MQ=255)
                                                        gCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTccc                                                                              >  1:407422/1‑140 (MQ=255)
                                                         cAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTcccc                                                                             <  2:685265/140‑1 (MQ=255)
                                                         cAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTcccc                                                                             <  2:364058/140‑1 (MQ=255)
                                                         cAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTcccc                                                                             <  2:407422/140‑1 (MQ=255)
                                                         cAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTcccc                                                                             <  2:80835/140‑1 (MQ=255)
                                                              tttCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACTGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCccacca                                                                        <  2:171046/140‑1 (MQ=255)
                                                                      aCTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCAc                                                                >  1:1056526/1‑140 (MQ=255)
                                                                      aCTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCAc                                                                >  1:936856/1‑140 (MQ=255)
                                                                        tctTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTg                                                              >  1:28335/1‑140 (MQ=255)
                                                                        tctTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTg                                                              >  1:785677/1‑140 (MQ=255)
                                                                         ctTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTgg                                                             <  2:936856/140‑1 (MQ=255)
                                                                            tGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCgg                                                          >  2:1230469/1‑140 (MQ=255)
                                                                                  cGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCaaaa                                                    >  1:1008355/1‑140 (MQ=255)
                                                                                  cGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCaaaa                                                    >  1:1140978/1‑140 (MQ=255)
                                                                                    ttGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAAt                                                  >  2:567338/1‑140 (MQ=255)
                                                                                      ggTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTa                                                <  2:1008355/140‑1 (MQ=255)
                                                                                           cGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAAtgaa                                           >  1:195898/1‑140 (MQ=255)
                                                                                                 cTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAAtgaactgaac                                     >  1:433003/1‑140 (MQ=255)
                                                                                                  ttCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACg                                    <  2:443296/140‑1 (MQ=255)
                                                                                                  ttCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACg                                    >  1:245854/1‑140 (MQ=255)
                                                                                                   tCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACgg                                   <  1:567338/140‑1 (MQ=255)
                                                                                                     aaCACGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCa                                 >  1:367461/1‑140 (MQ=255)
                                                                                                      acacGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAg                                <  2:470862/140‑1 (MQ=255)
                                                                                                        acGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCt                              <  2:367461/140‑1 (MQ=255)
                                                                                                        acGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCt                              <  1:1141849/140‑1 (MQ=255)
                                                                                                        acGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCt                              <  2:195898/140‑1 (MQ=255)
                                                                                                         cGCAGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTc                             <  2:328528/140‑1 (MQ=255)
                                                                                                            aGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACg                          >  1:247498/1‑140 (MQ=255)
                                                                                                            aGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACg                          >  1:1057535/1‑140 (MQ=255)
                                                                                                            aGTTCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACg                          >  2:872814/1‑140 (MQ=255)
                                                                                                              ttCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGa                        <  2:245854/140‑1 (MQ=255)
                                                                                                              ttCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGa                        >  2:404335/1‑140 (MQ=255)
                                                                                                              ttCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGa                        <  2:166698/140‑1 (MQ=255)
                                                                                                               tCACCGCCCTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGAt                       <  2:28335/140‑1 (MQ=255)
                                                                                                                      ccTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTa                <  2:785677/140‑1 (MQ=255)
                                                                                                                      ccTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTa                <  1:404335/140‑1 (MQ=255)
                                                                                                                      ccTTCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTa                <  2:1026780/140‑1 (MQ=255)
                                                                                                                        ttCTCTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGTGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATc              >  2:170144/1‑140 (MQ=255)
                                                                                                                         tctcTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCg             <  2:799188/140‑1 (MQ=255)
                                                                                                                         tctcTACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCg             <  2:1056526/140‑1 (MQ=255)
                                                                                                                          ctctACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCGc            >  2:722455/1‑140 (MQ=255)
                                                                                                                           tctACTACGCCCTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCGcc           <  2:206838/140‑1 (MQ=255)
                                                                                                                                    cccTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCGCCATCATCAgg  >  1:1153744/1‑140 (MQ=255)
                                                                                                                                    cccTGCAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCGCCATCATCAgg  >  1:928991/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TTCGCCAGCACCCCGGCATCGCGCAAGGCAGCGCGTAGCGGAACAGTAAATTTACGGCAGAATTTCAGTAACTCTTTGCTTTCGTTGGTATCGGCAACTTCAACACGCAGTTCACCGCCCTTCTCTACTACGCCCTGTAACATGCCGACAATGGGGGTAATACGATCTTCTGGCGGCAAATGCTGCAGGAGTTCCCCCACCACAAACCACTGGCGGGCAAAAATTAATGAACTGAACGGCAGCTCACGGATTAACTTATCGCCATCATCAGG  >  NC_000913/2940828‑2941099

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: