Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,041,632 T→A 100% V107E (GTG→GAG)  ygfZ → folate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,6320TA99.7% 157.6 / ‑7.1 58V107E (GTG→GAG) ygfZfolate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (32/23);  minor base C (1/0);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGT  >  NC_000913/3041493‑3041772
                                                                                                                                             |                                                                                                                                            
gccCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGAc                                                                                                                                           >  1:525992/1‑143 (MQ=255)
gccCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGAc                                                                                                                                           >  1:220750/1‑143 (MQ=255)
             ccAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGCCTGAACTGAAAAAATATGCGGTATTCTATAAAGAGACCATCGCGCCAgac                                                                                                                              >  1:15201/1‑143 (MQ=255)
              cAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAgacg                                                                                                                             >  1:640213/1‑143 (MQ=255)
              cAAAGGTAAAATGTGGAGCAATTTACGTCTGTGCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAgacg                                                                                                                             <  1:437608/143‑1 (MQ=255)
              cAAAGGTAAAATGTGGAGCAATTTACGTCTGCTGGATATGCTCTATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAgacg                                                                                                                             <  2:866337/143‑1 (MQ=255)
               aaaGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGACGGCTACCACGGCCAGAAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAgacga                                                                                                                            <  1:888002/143‑1 (MQ=255)
                   gTAACATGTGGAGCAATTTCTGTCTGTTCCGCGACGGCGATTGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAg                                                                                                                        <  2:46115/143‑1 (MQ=255)
                      aaaTGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTACAGCGACCATCGCGCCAGACGAGGAGCgt                                                                                                                     >  1:102605/1‑143 (MQ=255)
                      aaaTGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgt                                                                                                                     >  1:895154/1‑143 (MQ=255)
                           tgGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGcgct                                                                                                                >  2:309731/1‑143 (MQ=255)
                                  ttttACTTTTTTTCCGCCACGGCTATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgt                                                                                                         <  1:839241/142‑1 (MQ=255)
                                   tttACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgtg                                                                                                        >  1:7553/1‑143 (MQ=255)
                                     tACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGCACTGCAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGGGCTGCTTGGCGtt                                                                                                      >  2:524361/1‑143 (MQ=255)
                                      aCGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTg                                                                                                     >  1:217295/1‑143 (MQ=255)
                                        gTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAACAATATCCGGCATTCCCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGcc                                                                                                   >  2:632386/1‑143 (MQ=255)
                                         tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                                  >  1:554518/1‑143 (MQ=255)
                                         tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                                  >  1:554574/1‑143 (MQ=255)
                                         tCTGCTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                                  <  1:284677/143‑1 (MQ=255)
                                           tGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGATGAGCGTGTGCTGCTTGGTGTTGCCggg                                                                                                >  2:842281/1‑142 (MQ=255)
                                           tGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGCGCGTGTGCTGCTTGGTGTGGCCggc                                                                                                >  2:269572/1‑142 (MQ=255)
                                                          gATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTCTTCAGgcgcgcgc                                                                                  >  1:144049/1‑143 (MQ=255)
                                                                 ttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTCTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTg                                                                           >  2:222641/1‑143 (MQ=255)
                                                                     aTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                                                                       >  1:57044/1‑143 (MQ=255)
                                                                      tGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaaa                                                                      >  1:208299/1‑143 (MQ=255)
                                                                       ggATTGGACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAt                                                                     <  2:709316/143‑1 (MQ=255)
                                                                        gATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGACAAAtc                                                                    <  1:60017/143‑1 (MQ=255)
                                                                          ttGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCCTGGTGTTGCCGGTCTTCCGGCGCCCGCCGCCCTGCCaaagccc                                                                  >  2:648207/1‑139 (MQ=255)
                                                                           tGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGGTTTCAGGCGCGCGTCGCGCTGGCAActctct                                                                 >  2:684941/1‑143 (MQ=255)
                                                                              aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCttt                                                               <  1:369738/143‑1 (MQ=255)
                                                                               aCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTa                                                              <  1:427085/143‑1 (MQ=255)
                                                                               aCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTGTa                                                              >  2:611847/1‑143 (MQ=255)
                                                                                cGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                                                             <  1:674491/143‑1 (MQ=255)
                                                                                 ggCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCACTAAAGAGACCATCGCAGCCGACGACGGGCGTGAGCGGCCTGGTGTTGCCGGATATCAGGCGCGCGCCGCGCCGGCAAAACACATTAGc                                                            >  2:775815/1‑143 (MQ=255)
                                                                                   cgcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa                                                          >  1:350310/1‑143 (MQ=255)
                                                                                    gcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGaa                                                         <  1:366041/143‑1 (MQ=255)
                                                                                     cAGCGTGCGTGAACCGCAGCGGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGGCGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAAc                                                        <  2:551580/143‑1 (MQ=255)
                                                                                            cgtgACCCTCGGCGGACTGAACTGAAATAATATTCGGTATTCTCTAACGAGACCATCGCGCTAGACGACGAGCGTGTGCTGCTTGGGGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTc                                                 <  2:659117/143‑1 (MQ=255)
                                                                                             gtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGGCGACCGCTTTAGCGAACTGCCTTCg                                                >  2:140479/1‑143 (MQ=255)
                                                                                              tgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCga                                               >  2:891850/1‑143 (MQ=255)
                                                                                               gAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaa                                              <  1:320941/143‑1 (MQ=255)
                                                                                                aaCCGCAGCCGACGGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                             <  1:857693/143‑1 (MQ=255)
                                                                                                 aCCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaag                                            <  1:460926/143‑1 (MQ=255)
                                                                                                   cGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTATCGAACTGCCTTCgaaagaa                                          >  1:38673/1‑143 (MQ=255)
                                                                                                   cGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaa                                          >  1:694311/1‑143 (MQ=255)
                                                                                                    gCAGCTGGGTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaaa                                         <  2:648539/143‑1 (MQ=255)
                                                                                                         tgactgaACTGAAAAAATATGCGGTATTCACTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTCTATCGCACTGGCTTCGCACGAACAACAg                                    >  2:163131/1‑143 (MQ=255)
                                                                                                                  tgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCACCAGACGACGCGCGTGTGCTGCTTGGTGTTGCCGGAGTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACCGCCCTCGAAAGACAAACAGGTAGtgaga                           >  2:107966/1‑139 (MQ=255)
                                                                                                                    aaaaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaaga                         >  2:711414/1‑143 (MQ=255)
                                                                                                                     aaaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaa                        <  2:895154/143‑1 (MQ=255)
                                                                                                                      aaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaag                       <  1:520626/143‑1 (MQ=255)
                                                                                                                            tGCGGTAGTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGa                 <  2:679835/143‑1 (MQ=255)
                                                                                                                                    tctcTAAAGAGTTCATCGCGCCAGACTCCCAGCTTGTGCTGCTTGGTGTTGCCTGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg         <  1:876297/143‑1 (MQ=255)
                                                                                                                                      tctAAAGAGACCATCGCGCCGGGCGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCt       <  1:211707/143‑1 (MQ=255)
                                                                                                                                      tctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCt       >  1:819702/1‑143 (MQ=255)
                                                                                                                                      tctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGATTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAAGTAGTCAAAAAAAGCGCAACCAGTTTACt       >  2:23342/1‑143 (MQ=255)
                                                                                                                                         aaaGGGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATg    <  1:269572/143‑1 (MQ=255)
                                                                                                                                           aGTGACCATCGCGCCAGTCGACGTGCGTGTGCTGCTTGGTGTGGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTGGTCAAAGTAGGCGCGACCACTTTGCTATGGt  <  1:233542/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
GCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGT  >  NC_000913/3041493‑3041772

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