Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 935,434 A→G E737E (GAA→GAG)  ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913935,4340AG98.2% 202.9 / ‑6.4 57E737E (GAA→GAG) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (36/20);  minor base C (0/1);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CAGTTTGCGCAAACTCAACAACAACGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCG  >  NC_000913/935288‑935581
                                                                                                                                                  |                                                                                                                                                   
cAGTTTGCGCAAACTCAACTACAACGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGACCCTg                                                                                                                                                 <  2:19004/151‑1 (MQ=255)
          aaaCTCAACAACAACGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGcc                                                                                                                                       >  1:389664/1‑151 (MQ=255)
               caacaacaaCGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGcaac                                                                                                                                  >  2:133987/1‑151 (MQ=255)
                aacaacaacGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGcaaca                                                                                                                                 <  1:379911/151‑1 (MQ=255)
                    acaacGTTATTCCGGCGAACAACCGGCTGCGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGACCCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCg                                                                                                                             <  1:708139/151‑1 (MQ=255)
                     caacGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCgg                                                                                                                            >  2:405452/1‑151 (MQ=255)
                      aacGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGt                                                                                                                           >  2:278273/1‑151 (MQ=255)
                       acGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGtt                                                                                                                          <  2:925315/151‑1 (MQ=255)
                       acGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGtt                                                                                                                          >  2:474383/1‑151 (MQ=255)
                         gTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGc                                                                                                                        <  1:116668/151‑1 (MQ=255)
                            aTTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCAcc                                                                                                                     >  2:216131/1‑151 (MQ=255)
                             ttCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCg                                                                                                                    >  2:77229/1‑151 (MQ=255)
                              tCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCgc                                                                                                                   >  2:26051/1‑151 (MQ=255)
                              tCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCg                                                                                                                    <  2:335501/150‑1 (MQ=255)
                               ccGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCgca                                                                                                                  <  2:1060691/151‑1 (MQ=255)
                                   cGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCgcagc                                                                                                                >  1:1008605/1‑149 (MQ=255)
                                     aacaacCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCgcagcagca                                                                                                            <  1:29037/151‑1 (MQ=255)
                                        aacCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAat                                                                                                          >  1:827603/1‑150 (MQ=255)
                                           cGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATca                                                                                                      <  2:157534/151‑1 (MQ=255)
                                               tGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATcagc                                                                                                    >  2:317587/1‑149 (MQ=255)
                                                ggggCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATcagca                                                                                                   <  2:722977/149‑1 (MQ=255)
                                                 gggCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGcc                                                                                                >  1:556620/1‑151 (MQ=255)
                                                             tctcGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGt                                                                                    >  1:428126/1‑151 (MQ=255)
                                                             tctcGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGt                                                                                    >  2:884730/1‑151 (MQ=255)
                                                                cGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTc                                                                                  <  1:301574/150‑1 (MQ=255)
                                                                 gCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTccg                                                                                >  1:658012/1‑151 (MQ=255)
                                                                    ggATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTccgcc                                                                              <  1:77230/150‑1 (MQ=255)
                                                                        gatTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagc                                                                         >  1:9436/1‑151 (MQ=255)
                                                                         atTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagcc                                                                        >  1:830798/1‑151 (MQ=255)
                                                                          ttttGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagccg                                                                       >  1:960818/1‑151 (MQ=255)
                                                                              gAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGt                                                                   >  2:494224/1‑151 (MQ=255)
                                                                                aTTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTAt                                                                 <  1:929540/151‑1 (MQ=255)
                                                                                   ttCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATc                                                                >  2:361624/1‑149 (MQ=255)
                                                                                    tCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagc                                                             >  1:169723/1‑151 (MQ=255)
                                                                                     cGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagca                                                            <  2:980308/151‑1 (MQ=255)
                                                                                       ccAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCagcagc                                                          <  1:309322/151‑1 (MQ=255)
                                                                                        cAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGcc                                                         >  1:260952/1‑151 (MQ=255)
                                                                                                  gTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCg                                               >  2:932774/1‑151 (MQ=255)
                                                                                                       tGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                          >  1:714877/1‑151 (MQ=255)
                                                                                                                  ccACACGAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCcac                               <  1:317597/151‑1 (MQ=255)
                                                                                                                   cacaCGATCCGTGGTTTACGCAAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCCCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                 <  2:623079/148‑1 (MQ=255)
                                                                                                                     cgccaaCCGTTTTTTACGCCAATTGTTGCGCTTGTACGGCAGCGGAAACAACCGGTTGCACGGAAGCAGAAATATCCGCAGCCGCAACAACCCGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCcaca                              <  2:650154/145‑1 (MQ=18)
                                                                                                                        gAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATc                         >  1:359678/1‑151 (MQ=255)
                                                                                                                        gAACCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATc                         >  1:926072/1‑151 (MQ=255)
                                                                                                                         aaCCGTTGTTTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATc                         >  2:754369/1‑150 (MQ=255)
                                                                                                                                gttTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagcag                    >  1:1000464/1‑148 (MQ=255)
                                                                                                                                gttTACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGcc                  >  2:45425/1‑150 (MQ=255)
                                                                                                                                  ttACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGc                >  2:763146/1‑150 (MQ=255)
                                                                                                                                  ttACGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGc                >  2:214053/1‑150 (MQ=255)
                                                                                                                                    aCGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaac             >  1:296493/1‑151 (MQ=255)
                                                                                                                                    aCGCCAATTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaac             <  1:466001/151‑1 (MQ=255)
                                                                                                                                         aaTTGTTGAGCCTGTACAGCACCCGCAACACCCGGTTCCACCGCAGCAGAAATATCAGCCGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCg        <  2:592593/151‑1 (MQ=21)
                                                                                                                                          aTTGTTGAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCgg       <  1:326723/151‑1 (MQ=255)
                                                                                                                                           ttgttgAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt      >  1:131125/1‑151 (MQ=255)
                                                                                                                                            tgttgAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGtt     >  1:888957/1‑151 (MQ=255)
                                                                                                                                              ttgAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgc   >  2:370730/1‑151 (MQ=255)
                                                                                                                                               tgAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgcg  >  1:752072/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGTTTGCGCAAACTCAACAACAACGTTATTCCGGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCG  >  NC_000913/935288‑935581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑