Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,207,404 T→C A570A (GCT→GCC yehM → putative YehN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,207,4040TC100.0% 213.5 / NA 61A570A (GCT→GCCyehMputative YehN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/37);  total (24/37)

TATGGCGGCTGGATGAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCTCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGCCCCGGCGC  >  NC_000913/2207274‑2207537
                                                                                                                                  |                                                                                                                                     
tATGGCGGCTGGATGAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCAt                                                                                                                              >  2:58272/1‑140 (MQ=255)
 aTGGCGGCTGGATGAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATa                                                                                                                             >  2:459909/1‑140 (MQ=255)
     cggcTGGGTGAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTg                                                                                                                         >  1:962668/1‑140 (MQ=255)
       gcTGGATGAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGc                                                                                                                       >  2:517573/1‑140 (MQ=255)
            aTGAAATTAGCGGTATGCAAGGCGCGC‑‑ATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCt                                                                                                                <  2:655336/140‑1 (MQ=255)
              gAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGAACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCt                                                                                                                <  2:278883/140‑1 (MQ=255)
                   tAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGcc                                                                                                           <  1:459909/140‑1 (MQ=255)
                            gCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  1:110678/140‑1 (MQ=255)
                            gCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  2:791173/140‑1 (MQ=255)
                            gCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  2:399512/140‑1 (MQ=255)
                            gCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  1:210158/140‑1 (MQ=255)
                            gCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  1:103080/140‑1 (MQ=255)
                             cAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGTCATATTCTgcg                                                                                                 <  2:962668/140‑1 (MQ=255)
                             cAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgcg                                                                                                 <  1:517573/140‑1 (MQ=255)
                             cAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgcg                                                                                                 <  1:374632/140‑1 (MQ=255)
                              aaGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgcgc                                                                                                <  1:482762/140‑1 (MQ=255)
                              aaGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgcgc                                                                                                <  1:845302/140‑1 (MQ=255)
                                  cgcgcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  >  1:445614/1‑134 (MQ=255)
                                  cgcgcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  >  1:130484/1‑134 (MQ=255)
                                  cgcgcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  2:445614/134‑1 (MQ=255)
                                  cgcgcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgc                                                                                                  <  2:130484/134‑1 (MQ=255)
                                    cgcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGcc                                                                                                           <  1:207/123‑1 (MQ=255)
                                    cgcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGcc                                                                                                           >  2:207/1‑123 (MQ=255)
                                     gcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGcgcg                                                                                               <  2:719306/134‑1 (MQ=255)
                                     gcAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGcgcg                                                                                               >  1:719306/1‑134 (MQ=255)
                                      aaGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGTGCGATCTACa                                                                                        >  2:919478/2‑140 (MQ=255)
                                       aGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgcgc                                                                                                >  1:363035/1‑131 (MQ=255)
                                       aGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTgcgc                                                                                                <  2:363035/131‑1 (MQ=255)
                                        gATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAg                                                                                     <  2:81697/140‑1 (MQ=255)
                                         atatTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAg                                                                                     <  2:950790/140‑1 (MQ=255)
                                                  acgacgTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTaa                                                                            >  1:444718/1‑140 (MQ=255)
                                                  acgacgTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAACGGCGTTaa                                                                            >  1:747249/1‑140 (MQ=255)
                                                   cgacgTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAAt                                                                           >  1:168129/1‑140 (MQ=255)
                                                            gcgcGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGc                                                                  <  2:168129/140‑1 (MQ=255)
                                                                gACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGttt                                                              >  2:366116/1‑140 (MQ=255)
                                                                  cTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTtct                                                            >  2:504889/1‑140 (MQ=255)
                                                                    atcgatcgCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTctct                                                          <  2:756578/140‑1 (MQ=255)
                                                                         tcgCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTg                                                     <  2:302269/140‑1 (MQ=255)
                                                                            cACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTg                                                     <  2:534130/137‑1 (MQ=255)
                                                                            cACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTg                                                     >  1:534130/1‑137 (MQ=255)
                                                                              cGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTACGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTgcggcg                                                <  1:58272/140‑1 (MQ=255)
                                                                               gctgtggctgtgTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGcggcgg                                               >  2:554178/1‑140 (MQ=255)
                                                                                  gtggctgtgTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTATCTTTCTGCGGCGGTAg                                            <  1:161443/140‑1 (MQ=255)
                                                                                   tggctgtgTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTg                                                     <  2:683777/130‑1 (MQ=255)
                                                                                   tggctgtgTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTg                                                     >  1:683777/1‑130 (MQ=255)
                                                                                           tgAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTgg                                   >  2:700452/1‑140 (MQ=255)
                                                                                             aaTCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCCTACTGGAg                                 >  2:660378/1‑140 (MQ=255)
                                                                                                cTAACGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGa                              <  1:504889/140‑1 (MQ=255)
                                                                                                   aCGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAg                                            <  2:848062/123‑1 (MQ=255)
                                                                                                   aCGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAg                                            >  1:848062/1‑123 (MQ=255)
                                                                                                   aCGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGc                           >  1:514199/1‑140 (MQ=255)
                                                                                                   aCGGCAGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGc                           >  1:346843/1‑140 (MQ=255)
                                                                                                        aGACCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCa                      <  1:554178/140‑1 (MQ=255)
                                                                                                          aCCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAgg                    <  1:184185/140‑1 (MQ=255)
                                                                                                          aCCGGATGAAAAGGAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAgg                    <  1:700452/140‑1 (MQ=255)
                                                                                                                      ggAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGccc        <  2:514199/140‑1 (MQ=255)
                                                                                                                      ggAGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGccc        <  2:482680/140‑1 (MQ=255)
                                                                                                                        aGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGccc        <  1:202434/138‑1 (MQ=255)
                                                                                                                        aGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGccc        <  2:23081/138‑1 (MQ=255)
                                                                                                                        aGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGccc        >  1:23081/1‑138 (MQ=255)
                                                                                                                        aGTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGccc        >  2:202434/1‑138 (MQ=255)
                                                                                                                         gTTTCACGCCCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGTGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGCCCCgg     <  1:919478/140‑1 (MQ=255)
                                                                                                                            tCACGCCCACCTGCATAGTTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGCCCCGgcgc  >  1:754285/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                     
TATGGCGGCTGGATGAAATTAGCGGTATGCAAGGCGCGCAGATATTACAAACGACGTTATGCGCGACTATCGATCGCACGCTGTGGCTGTGTGAATCTAACGGCAGACCGGATGAAAAGGAGTTTCACGCTCACCTGCATAGCTGGCAAGCGCTTTGCCATATTCTGCGCGATCTACATAGCGGCGTTAATTTACCCGGCGTTTCTCTTTCTGCGGCGGTAGCCTTACTGGAGCGACGCAGTCAGGCAATTCATGCCCCGGCGC  >  NC_000913/2207274‑2207537

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