Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,726,152 T→G S44A (TCC→GCC)  nemR → DNA‑binding transcriptional repressor NemR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,1520TG98.3% 225.1 / ‑3.3 60coding (130/600 nt)nemRDNA‑binding transcriptional repressor NemR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (29/30);  minor base . (0/1);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

TCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGC  >  NC_000913/1726015‑1726285
                                                                                                                                         |                                                                                                                                     
tCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGcc                                                                                                                                     >  2:984157/1‑140 (MQ=255)
tCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGcc                                                                                                                                     >  1:1137616/1‑140 (MQ=255)
tCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGcc                                                                                                                                     >  1:131080/1‑140 (MQ=255)
tCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGcc                                                                                                                                     >  1:247992/1‑140 (MQ=255)
tCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGcc                                                                                                                                     >  1:337495/1‑140 (MQ=255)
 ccAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCt                                                                                                                                    <  2:247992/140‑1 (MQ=255)
    aCGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCt                                                                                                                                 >  2:939113/1‑140 (MQ=255)
    aCGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCt                                                                                                                                 >  2:365385/1‑140 (MQ=255)
    aCGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCt                                                                                                                                 <  2:1137616/140‑1 (MQ=255)
         tGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCac                                                                                                                            >  1:354694/1‑140 (MQ=255)
          gAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCact                                                                                                                           >  2:720199/1‑140 (MQ=255)
          gAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCact                                                                                                                           >  1:1146126/1‑140 (MQ=255)
          gAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCact                                                                                                                           >  1:597069/1‑140 (MQ=255)
            acaaacaCACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCactac                                                                                                                         >  1:929831/1‑140 (MQ=255)
                 cacacCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCg                                                                                                                    <  2:653292/140‑1 (MQ=255)
                 cacacCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCg                                                                                                                    <  1:365385/140‑1 (MQ=255)
                 cacacCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCg                                                                                                                    <  1:720199/140‑1 (MQ=255)
                  acacCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGCCTTCTATCACTACTTTCGct                                                                                                                  <  2:743864/140‑1 (MQ=255)
                      cGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTa                                                                                                               >  1:201265/1‑140 (MQ=255)
                      cGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTa                                                                                                               >  1:724415/1‑140 (MQ=255)
                          cATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAaaga                                                                                                           <  1:64554/140‑1 (MQ=255)
                          cATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAaaga                                                                                                           <  1:646790/140‑1 (MQ=255)
                               tACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGt                                                                                                      <  2:337495/140‑1 (MQ=255)
                                aCTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGtt                                                                                                     <  2:354694/140‑1 (MQ=255)
                                  tCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTg                                                                                                   >  1:1081537/1‑140 (MQ=255)
                                    gcgAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGc                                                                                                 >  2:957008/1‑140 (MQ=255)
                                     cgAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCg                                                                                                >  2:248434/1‑140 (MQ=255)
                                      gAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGt                                                                                               <  1:984157/140‑1 (MQ=255)
                                          aTCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGcc                                                                                           <  2:977533/140‑1 (MQ=255)
                                          aTCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGcc                                                                                           <  1:492538/140‑1 (MQ=255)
                                          aTCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGcc                                                                                           <  2:686098/140‑1 (MQ=255)
                                             tcCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATg                                                                                        <  2:597069/140‑1 (MQ=255)
                                                 ggCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTg                                                                                    <  1:248434/140‑1 (MQ=255)
                                                   cGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAg                                                                                  <  2:826785/140‑1 (MQ=255)
                                                    gACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGc                                                                                 >  1:260507/1‑140 (MQ=255)
                                                      cGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGt                                                                               >  2:1076357/1‑140 (MQ=255)
                                                              cAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTAcgc                                                                       <  2:1181508/140‑1 (MQ=255)
                                                                    tGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCata                                                                 <  1:696744/140‑1 (MQ=255)
                                                                            aCGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGc                                                         <  1:793801/140‑1 (MQ=255)
                                                                                   ttCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCgactgac                                                  <  1:245592/140‑1 (MQ=255)
                                                                                    tCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCgactgact                                                 >  2:817861/1‑140 (MQ=255)
                                                                                       ccGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAg                                              <  1:939113/140‑1 (MQ=255)
                                                                                        cGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGt                                             <  1:1182677/140‑1 (MQ=255)
                                                                                         gggATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGCGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGtt                                            >  1:1102900/1‑140 (MQ=255)
                                                                                           gATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTtgc                                          >  1:543005/1‑140 (MQ=255)
                                                                                               gggCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCa                                      >  2:783441/1‑140 (MQ=255)
                                                                                                ggCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCaa                                     >  1:297162/1‑140 (MQ=255)
                                                                                                       cGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCg                              <  2:1081537/140‑1 (MQ=255)
                                                                                                       cGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCg                              <  2:1002189/140‑1 (MQ=255)
                                                                                                       cGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCg                              <  1:1076357/140‑1 (MQ=255)
                                                                                                       cGAATTACTAAAAACCGCTGAAGTGCCGAAAGCGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCg                              <  2:1102900/140‑1 (MQ=255)
                                                                                                         aaTTACTAAAAACCGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTAATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGaa                            <  1:957008/140‑1 (MQ=255)
                                                                                                                      cGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACcgcg               <  2:276780/140‑1 (MQ=255)
                                                                                                                      cGCTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACcgcg               >  1:276780/1‑140 (MQ=255)
                                                                                                                        cTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGa              <  2:733561/139‑1 (MQ=255)
                                                                                                                        cTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGa              >  1:733561/1‑139 (MQ=255)
                                                                                                                        cTGAAGTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGAc             >  2:26492/1‑140 (MQ=255)
                                                                                                                             gTGCCGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCAt        >  1:1230197/1‑140 (MQ=255)
                                                                                                                               gCCGAAAGGGGCCCTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATAc      <  2:831362/140‑1 (MQ=255)
                                                                                                                                ccGAAAGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACt     >  1:110386/1‑140 (MQ=255)
                                                                                                                                   aaaGGGGCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGc  <  1:872834/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGC  >  NC_000913/1726015‑1726285

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: