Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,012,553 Δ1 bp intergenic (+160/+61) yqeA → / ← yqeB putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,5520G.100.0% 271.8 / NA 70intergenic (+159/+62)yqeA/yqeBputative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (36/34);  total (36/34)

CCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGC  >  NC_000913/3012416‑3012676
                                                                                                                                        |                                                                                                                            
ccATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAgca                                                                                                                            >  2:820917/1‑137 (MQ=255)
   tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAAtt                                                                                                                        >  2:1352742/1‑139 (MQ=255)
   tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTCGTAGGCCTGATAAGCGTAGCGCATCAGCAAtt                                                                                                                        >  2:1231916/1‑139 (MQ=255)
     cGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTg                                                                                                                      >  1:175766/1‑139 (MQ=255)
        ccAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCg                                                                                                                   >  1:701179/1‑139 (MQ=255)
            aGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTCGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTg                                                                                                               <  1:1231916/139‑1 (MQ=255)
                tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           <  1:101345/139‑1 (MQ=255)
                        aagtGCGCTTTGTTGATGCCGGATGCGGCGAGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCa                                                                                                   >  1:1037786/1‑139 (MQ=255)
                          gtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAgg                                                                                                 >  2:1216860/1‑139 (MQ=255)
                            gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGtt                                                                                                                 <  2:1545465/121‑1 (MQ=255)
                            gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGtt                                                                                                                 >  1:1545465/1‑121 (MQ=255)
                            gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGcc                                                                                               >  2:1408755/1‑139 (MQ=255)
                            gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGcc                                                                                               >  2:1476723/1‑139 (MQ=255)
                              gcTTTGTTGATGCGGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  2:1373741/139‑1 (MQ=255)
                              gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  1:1408755/139‑1 (MQ=255)
                              gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  1:1216860/139‑1 (MQ=255)
                              gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  2:500602/139‑1 (MQ=255)
                              gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  1:1352742/139‑1 (MQ=255)
                                       aTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                        <  2:411263/135‑1 (MQ=255)
                                       aTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                        >  1:411263/1‑135 (MQ=255)
                                            ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                               >  2:148218/1‑139 (MQ=255)
                                              aTGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATcc                                                                             >  1:1006751/1‑139 (MQ=255)
                                              aTGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc                                                                              <  2:1481786/138‑1 (MQ=255)
                                              aTGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc                                                                              >  1:1481786/1‑138 (MQ=255)
                                                                 tCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                          <  2:1037786/139‑1 (MQ=255)
                                                                    gTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACt                                                       <  1:820917/139‑1 (MQ=255)
                                                                      ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTc                                                     <  1:9533/139‑1 (MQ=255)
                                                                       cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           <  2:143731/84‑1 (MQ=255)
                                                                       cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           >  1:143731/1‑84 (MQ=255)
                                                                       cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTcc                                                    >  2:691907/1‑139 (MQ=255)
                                                                       cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTcc                                                    >  1:1331523/1‑139 (MQ=255)
                                                                         aCAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATTAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAg                                                  <  2:1331523/139‑1 (MQ=255)
                                                                         aCAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAg                                                  >  2:1276564/1‑139 (MQ=255)
                                                                          cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGc                                                 >  2:465094/1‑139 (MQ=255)
                                                                            aaaCATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCAc                                               <  2:701179/139‑1 (MQ=255)
                                                                                  cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCttt                                         >  2:762468/1‑139 (MQ=255)
                                                                                  cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCttt                                         >  1:905890/1‑139 (MQ=255)
                                                                                    cAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             >  2:936628/1‑85 (MQ=255)
                                                                                    cAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  1:936628/85‑1 (MQ=255)
                                                                                       aTTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                                                                    <  2:462941/91‑1 (MQ=255)
                                                                                       aTTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                                                                    >  1:462941/1‑91 (MQ=255)
                                                                                             tAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                              <  1:148218/139‑1 (MQ=255)
                                                                                                aTTGCAATGCACTTGTAGGTCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGc                           <  2:905890/139‑1 (MQ=255)
                                                                                                  tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                         >  2:1607936/1‑139 (MQ=255)
                                                                                                   gCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCt                        <  1:1607936/139‑1 (MQ=255)
                                                                                                         cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAgg                                                                                                 <  2:1305794/60‑1 (MQ=255)
                                                                                                         cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAgg                                                                                                 >  1:1305794/1‑60 (MQ=255)
                                                                                                         cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                  >  2:266609/1‑139 (MQ=255)
                                                                                                            ttGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           >  2:1034985/1‑47 (MQ=14)
                                                                                                            ttGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           <  1:1034985/47‑1 (MQ=14)
                                                                                                             tGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa              <  2:175766/139‑1 (MQ=255)
                                                                                                              gTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                         <  1:1449662/127‑1 (MQ=255)
                                                                                                              gTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                         >  2:1449662/1‑127 (MQ=255)
                                                                                                                   ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGcc                          <  1:792537/121‑1 (MQ=255)
                                                                                                                   ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGcc                          >  2:792537/1‑121 (MQ=255)
                                                                                                                    cTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGt                                      <  2:1544017/108‑1 (MQ=255)
                                                                                                                    cTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGt                                      >  1:1544017/1‑108 (MQ=255)
                                                                                                                    cTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCAtt       <  1:266609/139‑1 (MQ=255)
                                                                                                                      gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           <  1:65034/37‑1 (MQ=14)
                                                                                                                      gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                           >  2:65034/1‑37 (MQ=14)
                                                                                                                      gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             >  2:1613558/1‑51 (MQ=255)
                                                                                                                      gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                             <  1:1613558/51‑1 (MQ=255)
                                                                                                                      gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa     <  1:1476723/139‑1 (MQ=255)
                                                                                                                      gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa     >  1:734502/1‑139 (MQ=255)
                                                                                                                         aaGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGc  >  1:717600/1‑139 (MQ=255)
                                                                                                                            cGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGt                                                                                                       <  2:357195/35‑1 (MQ=17)
                                                                                                                            cGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGt                                                                                                       >  1:357195/1‑35 (MQ=14)
                                                                                                                                 cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgg                                                                                         <  1:1220510/44‑1 (MQ=25)
                                                                                                                                 cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgg                                                                                         >  2:1220510/1‑44 (MQ=25)
                                                                                                                                 cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                                                                          <  1:196181/59‑1 (MQ=255)
                                                                                                                                 cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                                                                          >  2:196181/1‑59 (MQ=255)
                                                                                                                                        |                                                                                                                            
CCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGC  >  NC_000913/3012416‑3012676

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