Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,498 1 bp→CG intergenic (+75/‑89) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,4971.C89.8% 117.3 / ‑0.9 52intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  major base C (18/26);  minor base A (5/0);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4980TG98.1% 146.8 / NA 52intergenic (+75/‑89)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (22/29);  minor base C (1/0);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG  >  NC_002947/1499359‑1499644
                                                                                                                                                ||                                                                                                                                                   
aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAAGCGCCCCCGTGCAATGCACCCGCTTTTCCGGTGTTAGCGGGCAc                                                                                                                                              >  8:43506/1‑148 (MQ=255)
aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCATGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGCCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAc                                                                                                                                              <  1:123558/148‑1 (MQ=37)
 cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAc                                                                                                                                              <  3:416373/147‑1 (MQ=37)
  gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcg                                                                                                                                             <  4:362553/147‑1 (MQ=37)
  gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGAAACAGGCGCCGCGGGGCAAGGCACCGGCTTTGCCGGTGTTCGCGGGCAcg                                                                                                                                             >  3:288984/1‑147 (MQ=37)
         cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcccgctc                                                                                                                                     <  4:69597/147‑6 (MQ=25)
         cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGTTCGCGGGCAcgcccgctc                                                                                                                                     <  3:31443/147‑6 (MQ=25)
          ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTcc                                                                                                                                     <  3:81482/147‑1 (MQ=25)
           cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCca                                                                                                                                   >  2:290464/1‑148 (MQ=25)
            gttaTGCGGCCATTCAGAGCCTCAACGACGCCTTGGCACCACAGCAGCACTGCGGCTGCCAACCGGGGGCCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCca                                                                                                                                   <  4:139383/144‑1 (MQ=18)
               aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                >  2:106935/1‑147 (MQ=21)
                    gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACACGCGCCCCGGTGCATTGCACCCGCTTTTCCGGTGTTCGCCGGGACCCCCGCGCCCACAGACTGAt                                                                                                                          >  7:20125/1‑148 (MQ=255)
                      cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTg                                                                                                                        >  2:289240/1‑148 (MQ=21)
                             agCCTGTACCAAGCGTTGGCCCAACCGCAGGATTGAGGCGGCCAACGGGGGGTCGGGGGCGTGCCCGCTAAACCGGCGCCGGGGGGCTTGGCACCGGCTTTGCCGGTGTTCGCGGCCACGCCCGCTCCCACAGACTGATTGCCATTaa                                                                                                                 <  8:189173/148‑1 (MQ=12)
                             agCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGGGGGGGATGGCACCGGCTTTGGCGGTGTTAGCGGGCACGCCCGCGGCCCCAGGCGGATTGCCAt                                                                                                                    >  3:67096/1‑145 (MQ=21)
                              gCCTGAACGAAGCCTTGGCACCACATCAGGACTGCGGCTGCAAAAGGGTGGGTGCGGGCGTGCCCGCGAAACTGGGGCCGCGGTGCATGGCACGGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACATACTGATTGCCAt                                                                                                                    <  4:288984/144‑1 (MQ=18)
                                ctcaaCGAAGCCTTGGCACCACAGCTGTACGGAGGCTGCCAAAGGGTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGC‑‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                                                              <  1:163926/145‑1 (MQ=18)
                                cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                                                               <  3:168591/147‑1 (MQ=21)
                                    aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCACGCCCCCTCCCACAGCCTGATTGCCATTAAGCGGGcc                                                                                                          >  2:55383/1‑148 (MQ=21)
                                        ccccTTGGCACCACCCCCGCACGCCGGCGCCCAACCGGTGGCCGCGGGCGGCCCCGCTAAACCGCCGCCGGGGGCCTTGCCACCGGCTTTCCGGGTGTTCGCGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                      <  8:13134/146‑1 (MQ=255)
                                                                aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                               <  7:118651/147‑1 (MQ=21)
                                                                aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                               <  1:125642/147‑1 (MQ=21)
                                                                aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                               <  1:290464/147‑1 (MQ=21)
                                                                aGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                                                <  7:228410/146‑1 (MQ=21)
                                                                 ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCCGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                             >  1:288673/1‑148 (MQ=21)
                                                                  gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                            <  1:55383/148‑1 (MQ=21)
                                                                     gCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCGCGGTGCATTGCACCCGCGTTTCCGGTGTTAGCGGGCACGCCCCCTCCCCCACACTGATTGCCAATTAGCGGGCCGATCGCTACAGTGACGGCCCGGtttttttttc                                                                          >  6:385937/1‑145 (MQ=255)
                                                                           cGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTGTTAGCCGGCACCCCCGCTCCCCCAGACAGATTTCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGccagggt                                                                    >  2:324319/1‑142 (MQ=21)
                                                                           cGGTGGGAGCGGGCGGGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGGGTGC‑‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTa                                                                   <  7:217203/148‑1 (MQ=21)
                                                                            ggggggAGGGGGCGGGCCCGCCAAACCGGCGCCCCGGCGCATGGCCCCGGCTTTGCCGGTGTTCGCGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                                                  <  7:148340/145‑1 (MQ=25)
                                                                                    cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                           >  3:358810/1‑147 (MQ=25)
                                                                                    cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                           <  5:385937/147‑1 (MQ=25)
                                                                                      tgCGTGCCCTCGAAACAGGAGCCGCGGTGCATGGCAACGGCTTTGCCGTTGTTCGCGGGGACGCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACATTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATg                                                         <  4:330144/146‑1 (MQ=255)
                                                                                             ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                   >  4:247275/1‑146 (MQ=37)
                                                                                              cgctAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCCCGCCCGCGCCCACCGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTGTTTGTTGCCTCGCTAACGAAGCATGACTccc                                                   >  5:270716/1‑145 (MQ=25)
                                                                                              cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                   >  3:5124/1‑145 (MQ=37)
                                                                                               gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                         >  7:326290/1‑122 (MQ=25)
                                                                                               gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                         <  8:326290/122‑1 (MQ=25)
                                                                                               gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                   >  5:406656/1‑144 (MQ=37)
                                                                                                cgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                >  4:284930/1‑146 (MQ=37)
                                                                                                cgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                              <  2:288673/148‑1 (MQ=37)
                                                                                                  aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc                                            <  7:43506/148‑1 (MQ=37)
                                                                                                         cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                      >  8:321553/1‑147 (MQ=37)
                                                                                                           ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                                      >  3:282748/1‑129 (MQ=25)
                                                                                                           ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                                      <  4:282748/129‑1 (MQ=25)
                                                                                                                gTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGa                               <  3:247275/147‑1 (MQ=38)
                                                                                                                   cATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTACCTTAGAACGCTATCAAGCAGATc                             >  1:409766/1‑146 (MQ=37)
                                                                                                                   cATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTACCTTAGAACGCTATCAAGCAGATc                             <  2:409766/146‑1 (MQ=37)
                                                                                                                    aTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGa                          <  1:289240/148‑1 (MQ=38)
                                                                                                                             ggCTTTGCCGGTGGTC‑‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                   <  4:67096/147‑1 (MQ=39)
                                                                                                                               cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCGCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa                >  2:38228/1‑148 (MQ=39)
                                                                                                                                    ccGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtc            >  6:30473/1‑147 (MQ=255)
                                                                                                                                              tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAt      <  3:284930/143‑1 (MQ=255)
                                                                                                                                                gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGa                               <  6:109137/118‑1 (MQ=255)
                                                                                                                                                gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGa                               >  5:109137/2‑119 (MQ=255)
                                                                                                                                                 cGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcg  <  6:406656/147‑1 (MQ=255)
                                                                                                                                                ||                                                                                                                                                   
ACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG  >  NC_002947/1499359‑1499644

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