Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,822 T→A 7.2% L50Q (CTG→CAG)  yegD → Hsp70 chaperone family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,8220TA 7.2% 147.6 / 2.6 58L50Q (CTG→CAG) yegDHsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base A (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGA  >  NC_000913/2147674‑2147969
                                                                                                                                                    |                                                                                                                                                   
gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCt                                                                                                                                                     <  4:329515/149‑1 (MQ=255)
   tttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTa                                                                                                                                                  <  5:270415/149‑1 (MQ=255)
   tttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTa                                                                                                                                                  >  2:84092/1‑149 (MQ=255)
           ttttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcat                                                                                                                                           >  6:169828/1‑148 (MQ=255)
               gATTACGATACAGCAAAGTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTTAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGAGAATGGCTGTACCGCCATCatg                                                                                                                                       <  2:49254/148‑1 (MQ=255)
                aTTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAACACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgat                                                                                                                                     <  6:11451/149‑1 (MQ=255)
                 ttACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgatg                                                                                                                                    <  8:101335/149‑1 (MQ=255)
                   aCGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGtt                                                                                                                                  >  1:342331/1‑149 (MQ=255)
                      gTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCg                                                                                                                               >  3:315064/1‑149 (MQ=255)
                       tACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCgg                                                                                                                              <  2:342331/149‑1 (MQ=255)
                          aGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgatg                                                                                                                                    >  4:166969/1‑140 (MQ=255)
                          aGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgatg                                                                                                                                    <  3:166969/140‑1 (MQ=255)
                            cAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAgac                                                                                                                         >  7:344477/1‑149 (MQ=255)
                             aaaCTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAgacg                                                                                                                        <  1:84092/149‑1 (MQ=255)
                                 tGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGAcgacga                                                                                                                    <  1:235857/149‑1 (MQ=255)
                                   ttCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGAcgacga                                                                                                                    <  8:344477/147‑1 (MQ=255)
                                       gTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAAc                                                                                                              >  8:174897/1‑149 (MQ=255)
                                                aTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCt                                                                                                     >  7:112260/1‑149 (MQ=255)
                                                 tGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCAAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTg                                                                                                    >  1:392053/1‑149 (MQ=255)
                                                   cGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGtt                                                                                                  <  2:162140/149‑1 (MQ=255)
                                                   cGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGt                                                                                                   <  3:231098/148‑1 (MQ=255)
                                                            aaaCCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGc                                                                                         >  5:47204/1‑149 (MQ=255)
                                                                cGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGAtt                                                                                     >  7:255331/1‑149 (MQ=255)
                                                                  cATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCg                                                                                   >  7:278306/1‑149 (MQ=255)
                                                                    tttGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGtt                                                                                 <  3:251310/149‑1 (MQ=255)
                                                                    tttGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCAGTACCGCCATCATGATGTTCCGGCAGACGACGATg                                                                                                                  >  7:224915/1‑116 (MQ=255)
                                                                    tttGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCAGTACCGCCATCATGATGTTCCGGCAGACGACGATg                                                                                                                  >  7:219521/1‑116 (MQ=255)
                                                                    tttGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCAGTACCGCCATCATGATGTTCCGGCAGACGACGATg                                                                                                                  <  8:219521/116‑1 (MQ=255)
                                                                    tttGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCAGTACCGCCATCATGATGTTCCGGCAGACGACGATg                                                                                                                  <  8:224915/116‑1 (MQ=255)
                                                                    tttGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGAGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGt                                                                                  >  5:373777/1‑148 (MQ=255)
                                                                          aaaaaTGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAAcc                                                                           >  8:165984/1‑149 (MQ=255)
                                                                           aaaaTGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAAcc                                                                           <  8:112260/148‑1 (MQ=255)
                                                                             aaTGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTCCCGCCATCAGGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAAAcgc                                                                         >  6:253886/1‑148 (MQ=255)
                                                                                   aaaCGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGaagaag                                                                  <  4:315064/149‑1 (MQ=255)
                                                                                       gACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGatat                                                              <  6:373777/149‑1 (MQ=255)
                                                                                        aCAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATc                                                             <  8:255331/149‑1 (MQ=255)
                                                                                          aGCACGCTGCTGCCTACAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGa                                                           >  6:24229/1‑149 (MQ=255)
                                                                                           gCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGAt                                                          <  5:24229/149‑1 (MQ=255)
                                                                                            cACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATg                                                         >  2:66616/1‑149 (MQ=255)
                                                                                               gctgctGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTa                                                      >  3:319645/1‑149 (MQ=255)
                                                                                                   ctgcCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGc                                                  <  8:278306/149‑1 (MQ=255)
                                                                                                        ttCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGaa                                               >  7:218960/1‑147 (MQ=255)
                                                                                                            aTGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCTAATGGCTGTACCGCCATCATGATGTTCCGGAAGAAGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCTCGAAGAAGATATCGATGTTACGGCGAAAAGCGt                                         <  5:253886/149‑1 (MQ=255)
                                                                                                            aTGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGt                                         <  4:319645/149‑1 (MQ=255)
                                                                                                             tGCTTTGCGCGCCAACGCGGGAAGCGGTAAGCGAATGGCGGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTg                                        <  2:252347/149‑1 (MQ=255)
                                                                                                                    cgcgcCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCg                                 <  5:169828/149‑1 (MQ=255)
                                                                                                                       gcCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGCCGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGAGGTTACGGCGAAAAGCGTGCAGTTCGGTc                              >  1:370362/1‑149 (MQ=255)
                                                                                                                             gcgTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCt                        <  4:301060/149‑1 (MQ=255)
                                                                                                                                tGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCAc                     <  2:125577/149‑1 (MQ=255)
                                                                                                                                      ggTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACt                    <  8:171193/144‑1 (MQ=255)
                                                                                                                                      ggTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACt                    >  7:171193/1‑144 (MQ=255)
                                                                                                                                       gTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGca                >  6:158705/1‑147 (MQ=255)
                                                                                                                                         atcatgaTGGCGGTACAGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAg             >  5:156622/7‑148 (MQ=255)
                                                                                                                                              aaTGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACAtt       >  5:86621/1‑149 (MQ=255)
                                                                                                                                                tGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATtga     <  7:165984/149‑1 (MQ=255)
                                                                                                                                                  gCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGGGCAGTTCGGTCTTTCCTCACTGGCACAGTACATtgatg   <  6:47204/149‑1 (MQ=255)
                                                                                                                                                   cTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACgg                                                   <  2:259972/100‑1 (MQ=255)
                                                                                                                                                   cTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACgg                                                   >  1:259972/1‑100 (MQ=255)
                                                                                                                                                   cTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATtgatga  <  2:392053/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTCGTTATAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGA  >  NC_000913/2147674‑2147969

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