Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,179,082 G→A 6.8% A164T (GCG→ACG)  nagK → N‑acetyl‑D‑glucosamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,179,0820GA 6.8% 186.5 / 2.4 59A164T (GCG→ACG) nagKN‑acetyl‑D‑glucosamine kinase
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (2/2);  total (29/30)
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.88e-01

CTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTA  >  NC_000913/1178934‑1179226
                                                                                                                                                    |                                                                                                                                                
cTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATg                                                                                                                                                  >  3:307665/1‑149 (MQ=255)
  gggATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCg                                                                                                                                                >  2:170094/1‑149 (MQ=255)
     atgatgAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTa                                                                                                                                             <  7:93104/149‑1 (MQ=255)
       gatgaATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCCGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAAc                                                                                                                                           >  1:284922/1‑149 (MQ=255)
        atgaATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAAc                                                                                                                                           >  5:145858/1‑148 (MQ=255)
           aaTGTACTCAATCTCCACTGGTGATGGGGGTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCatg                                                                                                                                       <  6:73199/149‑1 (MQ=255)
                          cACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTc                                                                                                                        <  5:161995/149‑1 (MQ=255)
                           aCTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTcc                                                                                                                       >  8:379531/1‑149 (MQ=255)
                           aCTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGAGGGGGTTGGATTTcc                                                                                                                       >  4:5866/1‑149 (MQ=255)
                                tGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTa                                                                                                                  <  1:170094/149‑1 (MQ=255)
                                 gATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTg                                                                                                                              <  1:281832/136‑1 (MQ=255)
                                 gATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTg                                                                                                                              <  7:146354/136‑1 (MQ=255)
                                 gATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTg                                                                                                                              >  2:281832/1‑136 (MQ=255)
                                 gATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTg                                                                                                                              >  8:146354/1‑136 (MQ=255)
                                       gTTGATTCTCGGCACCGGTGTTGGCGGCTGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATgggg                                                                                                                                 <  7:262802/127‑1 (MQ=255)
                                       gTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATgggg                                                                                                                                 >  8:262802/1‑127 (MQ=255)
                                                cGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTaa                                                                                                                                            >  3:258962/1‑107 (MQ=255)
                                                cGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTaa                                                                                                                                            <  4:258962/107‑1 (MQ=255)
                                                 ggCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCtgtg                                                                                                  <  6:245176/148‑1 (MQ=255)
                                                 ggCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCtg                                                                                                    <  3:5866/146‑1 (MQ=255)
                                                    aCCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGt                                                                                                >  6:283158/1‑147 (MQ=255)
                                                    aCCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGt                                                                                                >  1:349092/1‑147 (MQ=255)
                                                             ggcggcGGGCTGATTTTCAATGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTg                                                                                     <  2:335877/149‑1 (MQ=255)
                                                              gcggcgGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCAtatggc                                                                                     >  6:221709/1‑143 (MQ=255)
                                                                ggcggGCCGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCAt                                                                                  <  2:284922/149‑1 (MQ=255)
                                                                    ggCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGaa                                                                              <  7:283555/149‑1 (MQ=255)
                                                                     gCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGaa                                                                              >  4:183099/1‑148 (MQ=255)
                                                                         aTTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTa                                                                         <  4:307665/149‑1 (MQ=255)
                                                                          ttttCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTAt                                                                        <  6:145858/149‑1 (MQ=255)
                                                                               aaCGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTAtctgt                                                                    <  4:144453/148‑1 (MQ=255)
                                                                               aaCGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGAGGCTGTGGTCAGCATGGCTGCATTGAAAATTAtctgtc                                                                   >  4:60181/1‑149 (MQ=255)
                                                                                      aaCCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGt                                                                                                                                               <  6:214734/66‑1 (MQ=255)
                                                                                      aaCCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGt                                                                                                                                               >  5:214734/1‑66 (MQ=255)
                                                                                              aCCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGc                                                    >  7:74793/1‑149 (MQ=255)
                                                                                                cGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGCTGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTGTTGCGt                                                  >  8:73940/1‑149 (MQ=255)
                                                                                                cGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGc                                                    >  8:327096/1‑147 (MQ=255)
                                                                                                cGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGc                                                    <  7:327096/147‑1 (MQ=255)
                                                                                                cGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGt                                                  >  2:324287/1‑149 (MQ=255)
                                                                                                  gTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTgg                                                >  7:74195/1‑149 (MQ=255)
                                                                                                   tAAACGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGGTGAGGCGTTAACCATGATGGGGTTGGATGTCCCGTTACGCCGCTGCGGCTGTGGTCAGCTTGGCTGCATGGAAAATTATCTGTCGGGTCGCGGTTTTGCGTGGc                                               >  5:300801/1‑149 (MQ=255)
                                                                                                    aaaaGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCt                                              >  4:133460/1‑149 (MQ=255)
                                                                                                       aGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTa                                           >  1:347621/1‑149 (MQ=255)
                                                                                                          tACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATc                                         >  6:354645/1‑148 (MQ=255)
                                                                                                             aTTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAaca                                     <  5:283158/149‑1 (MQ=255)
                                                                                                              ttACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGAGTTACCCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAacac                                    >  6:95522/1‑149 (MQ=255)
                                                                                                                aCCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACt                                   <  3:183099/148‑1 (MQ=255)
                                                                                                                 ccGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGGGGCTGTATCAACACta                                  <  7:379531/148‑1 (MQ=255)
                                                                                                                  cGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATACGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAacac                                    <  2:46155/145‑1 (MQ=255)
                                                                                                                  cGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATACGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAacac                                    >  1:46155/1‑145 (MQ=255)
                                                                                                                   ggCGAGTTTGGCCATATGCGTCTGCCGGTTGATACGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGc                                                    >  8:348310/1‑128 (MQ=255)
                                                                                                                   ggCGAGGTTGGCCATATGCGTCTGCCGGTTGATACGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGc                                                    <  7:348310/128‑1 (MQ=255)
                                                                                                                        gTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGt                                                         >  1:123958/1‑118 (MQ=255)
                                                                                                                        gTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGt                                                         <  2:123958/118‑1 (MQ=255)
                                                                                                                        gTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTatcatc                          <  2:347621/149‑1 (MQ=255)
                                                                                                                            ggCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAAc                       <  1:101896/148‑1 (MQ=255)
                                                                                                                                    gCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAgg              >  7:386465/1‑149 (MQ=255)
                                                                                                                                             ggTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGaaa     <  1:324287/149‑1 (MQ=255)
                                                                                                                                              gTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACtatta                               <  4:256335/122‑1 (MQ=255)
                                                                                                                                                 gATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAatta  <  3:133460/148‑1 (MQ=255)
                                                                                                                                                 gATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAatta  <  8:74195/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTGGGATGATGAATTTACTCAATATCCACTGGTGATGGGGTTGATTCTCGGCACCGGTGTTGGCGGCGGGCTGATTTTCAACGGCAAACCGATTACCGGTAAAAGCTACATTACCGGCGAGTTTGGCCATATGCGTCTGCCGGTTGATGCGTTAACCATGATGGGGTTGGATTTCCCGTTACGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTA  >  NC_000913/1178934‑1179226

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: