Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 282,518 G→A 7.3% G81S (GGC→AGC)  yagE → CP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,5180GA 7.3% 143.0 / 3.4 57G81S (GGC→AGC) yagECP4‑6 prophage; putative 2‑keto‑3‑deoxygluconate aldolase
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-02

CGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAG  >  NC_000913/282376‑282655
                                                                                                                                              |                                                                                                                                         
cgacgcTCTGGTCAAGGCAGGCGTTGACGGCCTGTTCTGCCTGGGGAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGCTCCTGTCGATCGTCGCGTGCCGGTGCTGATCg                                                                                                                                           <  2:712306/143‑1 (MQ=255)
    gaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGGTCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAc                                                                                                                                       <  2:178934/143‑1 (MQ=255)
     aTCTGATCTAAGCTTGCGTTTACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAcc                                                                                                                                      <  1:345041/143‑1 (MQ=255)
       cTGATCAATGCAGGCGTTGACGCCGTGTTTTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCACcgg                                                                                                                                    <  2:972859/143‑1 (MQ=255)
        tGATCAAAGCAGGCGTTGACGGCCTGTCCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCACcggc                                                                                                                                   <  1:708001/143‑1 (MQ=255)
         gATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAACATGTCGATCGTCGCGTGCCGCTGCTGCTCGGCCCcggcg                                                                                                                                  >  2:89136/1‑143 (MQ=255)
          aTCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACcggcgg                                                                                                                                 >  1:843046/1‑143 (MQ=255)
             aaaGCAGGCGTTGACGGCCTGTCCGTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCAc                                                                                                                              <  2:721009/143‑1 (MQ=255)
               ggcgggCGCTGACGGCCTGTTCTGCCTGGGCAGCGGGGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGGCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCa                                                                                                                            <  1:381968/139‑1 (MQ=255)
               cgCAGGCGTTGACGGCCCGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGGTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCa                                                                                                                            <  1:1023449/142‑1 (MQ=255)
                 cAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAc                                                                                                                          <  1:838576/143‑1 (MQ=255)
                 cAGGCGTTGACGGCCTGTTCTGCCTGGGCAGGGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAc                                                                                                                          <  2:530433/143‑1 (MQ=255)
                    gCGTTGGCGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGcc                                                                                                                       <  1:383457/143‑1 (MQ=255)
                      gTTGACGGCCTGTTCTTCCTGGGCAGAGATGGCGAGTTCGCCCGGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTGTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCCCCGGCGGCACCAACACCCg                                                                                                                     >  1:730944/1‑143 (MQ=255)
                       ttGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGTCCgg                                                                                                                    <  2:538251/143‑1 (MQ=255)
                          agggCCTGGTCTTCCTGGGCAGCGGTGGCGAGTGCTCCCGGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTAGCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGaa                                                                                                                 <  2:825708/141‑1 (MQ=255)
                            ggCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGCGCTGATCGGCACCGGCGGCACCAACGCCCGGCAAAc                                                                                                               >  2:37642/1‑143 (MQ=255)
                            ggCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTACCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAc                                                                                                               >  1:611271/1‑143 (MQ=255)
                            ggCCTGTTCGTCCTGGGCAGCGGTGGGGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAc                                                                                                               <  2:139339/143‑1 (MQ=255)
                               cTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAccct                                                                                                            >  1:84701/1‑141 (MQ=255)
                               cTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCAt                                                                                                            <  1:866084/143‑1 (MQ=255)
                                    cttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTTGCGATCATGTCGCTCGTCTCGTGCCGGTGATGATCGGCAAAGGCGGGAACAACGATCGGGAAAACATCggaa                                                                                                       >  1:342886/1‑140 (MQ=255)
                                    cttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAGAGGAGGCACCACCGCCCCGGAAACCATCGAAc                                                                                                       >  2:215111/1‑143 (MQ=255)
                                     ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACt                                                                                                      <  1:870218/143‑1 (MQ=255)
                                       ccggggCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCAGTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTca                                                                                                    <  2:592560/140‑1 (MQ=255)
                                         tGGGTAGCGGTGGCTAGTTCTCCCAGCTCGGCCCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagc                                                                                                  <  1:761585/143‑1 (MQ=255)
                                         tGGGCAGCGGTGGCGAGTTCCCCCAGCGCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCAGGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagc                                                                                                  <  1:740970/143‑1 (MQ=255)
                                               gCGGTGGCGAGTTCTCCCCGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAc                                                                                            <  2:289889/143‑1 (MQ=255)
                                                   tGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAACGTCGCGTGCCGGTGCTGATCCGCCCCGGCTGCACCAACGCACGGGACACCATCGGACCAAACCAGCAcgcac                                                                                        >  2:296512/1‑141 (MQ=255)
                                                      cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCgcagc                                                                                     >  1:79428/1‑143 (MQ=255)
                                                               ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCAGCAACGCCCGGGAAGCCACCGACCTCAACCCGCACCCGCCACAggcgggcg                                                                            >  2:65448/1‑143 (MQ=255)
                                                                      ggCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCGGAGGCATCTTTCGAACTCTGCCAGCACGCGCCGCAGGCGGGCGCGGACgg                                                                     >  1:289404/1‑143 (MQ=255)
                                                                       gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCGTCGAACTCCGCCAGCACGCGCAGCAGGCGGGCGAGGCCGGc                                                                    >  2:281875/1‑143 (MQ=255)
                                                                       gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTACTGATCGGCACCGGCGGCACCCACGCCCGGGAAACCCTCGGACTCCGCCCGCACGCGCAGCAGGCGGGCGCGGACGGc                                                                    >  2:481734/1‑143 (MQ=255)
                                                                         gcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGCACTCAGCCAGCACGCGCGGCGGGCGCGGGCGGACGGCAt                                                                  >  1:220593/1‑143 (MQ=255)
                                                                               gagCGTAAAGCCATGGCCCGCGTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCAGCGGCACCAACGCCCGGGAAACCATCGAGCTCAGCCAGCACGCGCAGCCGGTGGTCGCGGACGGCATCgtggt                                                            <  2:647614/143‑1 (MQ=255)
                                                                                 gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGa                                                          <  2:197779/143‑1 (MQ=255)
                                                                                    tAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCTGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGAGGGCGCGGACGGCATCGTGGTGATCa                                                       >  1:247231/1‑143 (MQ=255)
                                                                                    tAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGCTCGGCACCGGCGGCACCAACGCCCGGGCAACCATCGGGCTCAGCCAGCACGCGCAGCAGGCGGGCGCGCACGGCGTCgcggtgcgca                                                       >  2:617961/1‑139 (MQ=255)
                                                                                     aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGCTCGGCACCGGCGGCACCAACGCCCGGGAGACCCTCGGGCTCAGCCAGCCCGCGCAGCAGGCGGGCGCGGGCGGCATCgggggggacaa                                                      >  1:25993/1‑136 (MQ=255)
                                                                                           aTTGGCCGCTGTGCTGTCGCTCGTGTCGATCGTGGCGTGCCGGTGTTGGTCGGCACCGGCGGCACCAACGCCCGGGACACCATGGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCcta                                                <  2:545222/143‑1 (MQ=255)
                                                                                                   cTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCACCGCCCGGGACACCTTCGAACTCAGCCAGCACGCACAGCAGGAGAGCACGGACGGCATCATAGGGATCAACAGATACTACTgaa                                        >  1:322172/1‑141 (MQ=255)
                                                                                                    tttGCTATCGGTCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                       <  1:602454/143‑1 (MQ=255)
                                                                                                    tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCATAATAATGGac                                       >  1:765008/1‑142 (MQ=255)
                                                                                                     ttGCTAGCGAGCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaaa                                      <  2:662293/143‑1 (MQ=255)
                                                                                                          atcgatcgTGCCGATCTTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                 <  2:731357/143‑1 (MQ=255)
                                                                                                               tCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCCGCACGCGCAGCAGGCGGGCGCGGCCGGCATCGTGGTGATCAACCCCTCCTACTGGAAAGTGTCGGagg                            >  1:1033944/1‑141 (MQ=255)
                                                                                                               tCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAg                            <  1:916509/143‑1 (MQ=255)
                                                                                                                    tCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCGGCAGGCGGGCGCCGACGGCATCGTGGTGACCCACCCCGACTACCGTAAAGGGTGCGAAggggac                       >  2:555188/1‑138 (MQ=255)
                                                                                                                          gtcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATc                 >  1:203479/1‑143 (MQ=255)
                                                                                                                            cgcgTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCg               <  1:993738/143‑1 (MQ=255)
                                                                                                                              cgTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCCTCGGGGTGATCAACCCCTACGCCTGGGAAGGGACGGAAGCGCAACTGacgcgcc             >  2:658655/1‑137 (MQ=255)
                                                                                                                                cgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGAGCAACCCCTACCACTGGAAAGTGTCGGAAGCGACCCTGATCCGccac           >  2:160932/2‑140 (MQ=255)
                                                                                                                                cgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGACGGGGTTGGCATCCTGGAGAGCAACCCCTACTACTGGAAAGTGTCGGAAGCGACCCTGGTGCGCCAt           >  1:289760/2‑143 (MQ=255)
                                                                                                                                  gcGGTGCTGATCGGCACCGGCGGCACCAACGGCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTAttt         <  2:569811/142‑1 (MQ=255)
                                                                                                                                        cTGATCGGGACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCa   <  2:414753/143‑1 (MQ=255)
                                                                                                                                        cTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTCCTGGAAAGTGTCGGAAGCGAACCTGATCCGCTGTTTCGAgta   >  1:135700/1‑141 (MQ=255)
                                                                                                                                         tGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGAGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGGAGCGAACCTGATCCGCCATCTCGAGCAg  >  2:529894/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                         
CGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAG  >  NC_000913/282376‑282655

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