Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,126,6451.C92.6% 161.6 / 7.2 54intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (18/32);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 2.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGC  >  NC_002947/1126498‑1126782
                                                                                                                                                       |                                                                                                                                          
tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGGGGAGGTATAGCCACCAAAGCAATAGCGGGGaaaa                                                                                                                                            >  2:97164/1‑148 (MQ=255)
     gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATTCGGAGGTATAGCCCCCATAGCATTAGCGGGTAAAcccgca                                                                                                                                      >  8:393410/1‑145 (MQ=255)
     gggCCGATGCCGATCTTGAACAGGTCGAAGACGATCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATGCGGAGGTTTAGCCACCATCGCAATAGCGGGTTAAcaccct                                                                                                                                      >  1:55355/1‑144 (MQ=255)
      ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                     <  1:161195/149‑6 (MQ=255)
          gatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTac                                                                                                                                  <  5:386724/149‑1 (MQ=255)
                 atCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATGCGGAGGTATAGGCACCAAAGCATTAGCCGGTAAAACCCCTCCAACACAGAGc                                                                                                                           >  6:195104/1‑149 (MQ=255)
                   cTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCGAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATGCGGGGGTATAGCCACCATAGCATTAGCGGGTTAAGCCGCTGCGACACAGAGccc                                                                                                                         >  5:186225/1‑149 (MQ=255)
                        aCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAc                                                                                                                    <  7:15822/149‑1 (MQ=255)
                               gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGCTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                              >  3:356798/1‑148 (MQ=255)
                                 aGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAg                                                                                                           <  2:366738/149‑1 (MQ=255)
                                 aGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAg                                                                                                           <  5:195104/149‑1 (MQ=255)
                                  gACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGCTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAgg                                                                                                          <  4:356798/149‑1 (MQ=255)
                                   aCGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCTTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGt                                                                                                         <  1:97164/149‑1 (MQ=255)
                                         aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATTGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                     <  1:123718/147‑1 (MQ=255)
                                         aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                     <  4:331766/147‑1 (MQ=255)
                                           tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                 <  4:306157/149‑1 (MQ=255)
                                             aCAGGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                <  7:171911/148‑1 (MQ=255)
                                               aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCCCCATAGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCCCCACCCAACCCAGGTGGCGGGGAggg                                                                                              >  6:218812/1‑148 (MQ=38)
                                               aTGTGGACCTCCTGCCAGGGTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCAGCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGc                                                                                             <  7:393410/149‑1 (MQ=38)
                                                     aCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATGCGGAGGGATAGCCACCATCGCAATAGCGGGTAAACCCGCACCTACACAGATCCCCCCAACCACACCAGGGGGCGGGGAGGGCTATCAg                                                                                       >  3:82657/1‑149 (MQ=255)
                                                      cctcgtGCCAGGTTTTTGGGTGTACATTTTGCCTTTTAGGCGGGGGTTTACCGGCGAATGGGTAGTTATAGCAACAATCGAATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                       <  8:58498/148‑1 (MQ=255)
                                                      cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGGAAACCCCCTCCTACACAGAGCCCGCCCCCCAACCCACGGGGCGGGGTCGGCTAtaatg                                                                                      >  7:125115/1‑145 (MQ=37)
                                                      cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAAACCCCTCCGACACAGAGCCCGCCACCCAACCCACGTGGCCGGGCGGGCTAACAgg                                                                                      >  5:326140/1‑149 (MQ=255)
                                                          ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt                                                                                   <  7:97757/148‑1 (MQ=255)
                                                          ctGCCAGGTTGGTGGGTGGACATGGAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGCAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCg                                                                                    <  2:370429/147‑1 (MQ=38)
                                                             ccAGTTTGTTGGGTGCACATTGAGCCCTGTAGGCGCGGTTTTACCCGCGAATGCGGATGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                 <  3:201301/147‑1 (MQ=38)
                                                               aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCCAATTCGGCGGTATAGCCCCCATCGCATTAGCGGGTAAACCCGCACCTACACAGACCCCGCCACCCAACCCAGGTGGCGGGGAGGGCTATAAGGCGtgcaga                                                                              >  2:47675/1‑143 (MQ=255)
                                                               aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                              <  1:47675/148‑1 (MQ=255)
                                                                    gttgGGTGGACATTGAGCCCTGTAGGGGCGGGTTTACGCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAGcc                                                                         <  6:186225/148‑2 (MQ=37)
                                                                          tGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                 <  2:55355/149‑1 (MQ=255)
                                                                               aTTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATTCGGATGTATTTCCCCCATACCCTTATCGGGGAAAACCCCGACAAAACACAGCCCCCCCCCCCAACCACGGGGCGGGGAGGGCTATAAAGCCTCCTGATATATCACCATCGAccg                                                            >  6:343694/1‑147 (MQ=255)
                                                                                         gTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                   >  6:208827/1‑148 (MQ=255)
                                                                                            ttcgcgGGTAAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcga                                                                   >  3:273037/3‑129 (MQ=38)
                                                                                            ttcgcgGGTAAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcga                                                                   <  4:273037/127‑1 (MQ=38)
                                                                                             gcgcGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                              <  5:208827/149‑1 (MQ=255)
                                                                                                    tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                       <  7:192861/149‑1 (MQ=255)
                                                                                                     ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                       <  4:82657/148‑1 (MQ=255)
                                                                                                         ccGCGAAGGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                  <  6:326140/149‑1 (MQ=255)
                                                                                                                ttcGGAGGTATTGCCACCATCGCATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCGACCCATGGCAGTTGGCGGGGAGGTCTATCAg                                                                                       <  4:351636/88‑1 (MQ=21)
                                                                                                                tGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                       >  3:351636/1‑90 (MQ=255)
                                                                                                                   ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGGTGCGGTCGCCGAACACGtt                        >  4:46327/1‑149 (MQ=255)
                                                                                                                     aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                      >  8:17095/1‑149 (MQ=255)
                                                                                                                     aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAAAGCTCTCGATCGAAGGGCAGGCGCACACCAGATTACGGTCGCCGAACACGATga                      >  7:133073/1‑148 (MQ=39)
                                                                                                                         ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                  >  1:42620/1‑149 (MQ=255)
                                                                                                                          taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                 >  5:385453/1‑149 (MQ=255)
                                                                                                                            gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGCATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                <  1:108875/149‑1 (MQ=255)
                                                                                                                            gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGCATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                <  1:103427/149‑1 (MQ=255)
                                                                                                                            gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg               >  5:354401/1‑149 (MQ=255)
                                                                                                                             ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg              >  3:93995/1‑149 (MQ=255)
                                                                                                                              caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc             >  1:375989/1‑149 (MQ=255)
                                                                                                                                ccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg           <  4:93995/149‑1 (MQ=255)
                                                                                                                                      cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc     <  2:42620/149‑1 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg   <  8:396165/149‑1 (MQ=255)
                                                                                                                                          cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGCTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc  <  3:46327/148‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                 <  3:23599/130‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                 >  4:23599/1‑130 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg               <  4:405423/132‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg               >  3:405423/1‑132 (MQ=255)
                                                                                                                                                       |                                                                                                                                          
TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGC  >  NC_002947/1126498‑1126782

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