Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,168,700 A→G 7.0% intergenic (+217/+36) eda → / ← leuA KHG/KDPG aldolase/2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,168,7000AG 7.0% 159.7 / 2.8 57intergenic (+217/+36)eda/leuAKHG/KDPG aldolase/2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (2/2);  total (31/26)
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.26e-01

GTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGTACCGGGCG  >  NC_002947/1168558‑1168846
                                                                                                                                              |                                                                                                                                                  
gTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGCTCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGGTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGAt                                                                                                                                              >  5:70488/1‑149 (MQ=255)
     gCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAAGAAAAAAAGCCCGCATCGTTTGATGCggg                                                                                                                                         <  1:196968/149‑1 (MQ=37)
         gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCttt                                                                                                                                     >  4:146007/1‑149 (MQ=255)
           ccTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGGGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGc                                                                                                                                        <  4:313310/144‑1 (MQ=255)
           ccTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGGGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGc                                                                                                                                        >  3:313310/1‑144 (MQ=255)
            cTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCGTTTGATGCggg                                                                                                                                         >  4:164490/1‑142 (MQ=255)
            cTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCGTTTGATGCggg                                                                                                                                         <  3:164490/142‑1 (MQ=255)
                   gggTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCt                                                                                                                            >  4:48520/1‑148 (MQ=255)
                           ccGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATttcagg                                                                                                                   >  7:169936/1‑149 (MQ=255)
                           ccGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATttcag                                                                                                                    <  8:89117/148‑1 (MQ=255)
                           ccGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATttcag                                                                                                                    >  7:89117/1‑148 (MQ=255)
                            cGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATttcaggc                                                                                                                  <  8:169936/149‑1 (MQ=255)
                                              tGCGTTGGGTACACTCCGGGGTTTTGACACAGGCTTTAGTTGGAGGGGGTTTACCCGCGAATGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCt                                                                                                 <  2:280281/148‑1 (MQ=255)
                                                cgttGGGAACACACCGGGTTTTGGACACAGGCCTTAGTAGCAGCGGTTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCtt                                                                                                     <  7:274727/139‑1 (MQ=255)
                                                cGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGTAGCGGGTTTACCCGCTAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTc                                                                                               <  5:284636/148‑1 (MQ=255)
                                                        acacacCGGGTGTTGGCCACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCt                                                                                      <  4:338969/149‑1 (MQ=255)
                                                        acacacCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCGCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGc                                                                                       >  5:61048/1‑148 (MQ=255)
                                                          acacCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGCGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCt                                                                                      >  3:212579/1‑147 (MQ=255)
                                                                    ttGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGt                                                                           >  6:49376/1‑148 (MQ=255)
                                                                      ggACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGtt                                                                          >  4:164477/1‑147 (MQ=255)
                                                                       gACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCa                                                                        <  5:49376/148‑1 (MQ=255)
                                                                        acacAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAgc                                                                      >  8:276658/1‑149 (MQ=255)
                                                                           caGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCTCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCa                                                                   <  6:61048/149‑1 (MQ=255)
                                                                            aGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCAc                                                                  >  7:50077/1‑149 (MQ=255)
                                                                              gCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACt                                                                 >  3:101315/1‑148 (MQ=255)
                                                                                     tagGAGCGGGTTTACCCGCGAAAGGGCCAGTACAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAgcg                                                         >  5:24097/1‑149 (MQ=255)
                                                                                       ggAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTc                                                                                                                           >  4:127416/1‑81 (MQ=255)
                                                                                       ggAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTc                                                                                                                           <  3:127416/81‑1 (MQ=255)
                                                                                          gCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTg                                                    <  3:164477/149‑1 (MQ=255)
                                                                                            gggTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCGGCGCACTGAACAGCGCCTTg                                                    >  2:244857/1‑147 (MQ=255)
                                                                                            gggTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTg                                                    >  7:35382/1‑147 (MQ=255)
                                                                                                 tACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTg                                              <  3:48520/148‑1 (MQ=255)
                                                                                                  aCCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTgg                                             <  8:35382/148‑1 (MQ=255)
                                                                                                    ccGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCg                                          >  4:95506/1‑149 (MQ=255)
                                                                                                          aaGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGa                                    <  8:50077/149‑1 (MQ=255)
                                                                                                          aaGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGa                                    >  2:6079/1‑149 (MQ=255)
                                                                                                               ccAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTtc                               >  8:129222/1‑149 (MQ=255)
                                                                                                               ccAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTtc                               <  6:70488/149‑1 (MQ=255)
                                                                                                                  gTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCa                            <  3:95506/149‑1 (MQ=255)
                                                                                                                   tCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTgg                                             <  1:141462/131‑1 (MQ=255)
                                                                                                                   tCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTgg                                             >  2:141462/1‑131 (MQ=255)
                                                                                                                   tCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCAt                           <  1:6079/149‑1 (MQ=255)
                                                                                                                        gCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATa                      >  7:224562/1‑149 (MQ=255)
                                                                                                                          caacaaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATAcc                    >  2:175132/1‑149 (MQ=255)
                                                                                                                          caacaaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATAcc                    >  1:172040/1‑149 (MQ=255)
                                                                                                                           aacaaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGAt                       >  3:111069/1‑145 (MQ=255)
                                                                                                                           aacaaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGAt                       <  4:111069/145‑1 (MQ=255)
                                                                                                                           aacaaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATAccc                   >  6:108985/1‑149 (MQ=255)
                                                                                                                             caaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGaa                                                  <  5:82624/116‑1 (MQ=255)
                                                                                                                             caaAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGaa                                                  >  6:82624/1‑116 (MQ=255)
                                                                                                                               aaaaaaGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGc               <  4:101315/149‑1 (MQ=255)
                                                                                                                                   aaGCCCGCATCGTTTGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATg           >  6:101876/1‑149 (MQ=255)
                                                                                                                                    aGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGt          <  4:212579/149‑1 (MQ=255)
                                                                                                                                      cccGCATCAAACGATGCGGGCTGTATCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGTAc        >  7:220736/1‑149 (MQ=255)
                                                                                                                                         gCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGTACCgg     <  6:24097/149‑1 (MQ=255)
                                                                                                                                          cATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGTACCtgg    <  8:224562/149‑1 (MQ=255)
                                                                                                                                            tCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGTACCGggcg  >  3:121026/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCAGGCTTCAGGCCGCCTTGGCTTCCTGCTGGCTCAGCGAGCGGTTCAGCGCACTGAACAGCGCCTTGAAGCTGGCCGTGGTGATGTTCTCATCGATACCCACGCCATGTACCGGGCG  >  NC_002947/1168558‑1168846

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: ‑