Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,416,884 A→G V168A (GTT→GCT)  recE ← Rac prophage, exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,416,8840AG96.5% 190.4 / ‑2.7 57V168A (GTT→GCT) recERac prophage, exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (29/26);  total (30/27)
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.80e-01

CTTCGAAAAAAGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAG  >  NC_000913/1416756‑1417019
                                                                                                                                |                                                                                                                                       
cTTCGAAAAAAGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGc                                                                                                                               >  1:576547/1‑139 (MQ=255)
     aaaaaaGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTaa                                                                                                                          >  2:401177/1‑139 (MQ=255)
     aaaaaaGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTAGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGAGGGAAACAGCCAGGGCAGTaa                                                                                                                          >  2:1169447/1‑139 (MQ=255)
       aaaaGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACa                                                                                                                        <  1:1607924/139‑1 (MQ=255)
          aGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGcc                                                                                                                     <  2:1329842/139‑1 (MQ=255)
               ttATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAg                                                                                                                >  1:503675/1‑139 (MQ=255)
                   cAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGttct                                                                                                            <  2:806472/139‑1 (MQ=255)
                   cAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGttct                                                                                                            <  1:401177/139‑1 (MQ=255)
                         cAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGcc                                                                                                      >  1:1235410/1‑139 (MQ=255)
                           ggTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGccc                                                                                                     <  2:1050364/138‑1 (MQ=255)
                           ggTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCg                                                                                                    <  2:1080961/139‑1 (MQ=255)
                           ggTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTAGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGccc                                                                                                     >  1:1050364/1‑138 (MQ=255)
                           ggTTTGAATTACCAGCATTAGGGAAAACTTTGTCAGTGTCACGAACCGGTTTGTGGTGTTCGCGAATTTCTAGCGGGTCGAGCAGGCTGGCTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCGTGCGCCCg                                                                                                    >  2:1602848/1‑139 (MQ=255)
                            gTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGa                                                                                                   <  1:1169447/139‑1 (MQ=255)
                                gAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAagcagc                                                                                                                                     >  2:842383/1‑101 (MQ=255)
                                gAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAagcagc                                                                                                                                     <  1:842383/101‑1 (MQ=255)
                                gAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGttcttc                                                                                                          >  1:1455859/1‑128 (MQ=255)
                                gAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGttcttc                                                                                                          <  2:1455859/128‑1 (MQ=255)
                                      ccAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGcc                                                                                                                                    >  2:6791/1‑96 (MQ=255)
                                      ccAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGcc                                                                                                                                    <  1:6791/96‑1 (MQ=255)
                                      ccAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATgtgt                                                                                          >  2:716606/1‑138 (MQ=255)
                                      ccAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATgtgt                                                                                          <  1:716606/138‑1 (MQ=255)
                                      ccAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTa                                                                                         >  2:1437342/1‑139 (MQ=255)
                                              aGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTgg                                                                                 >  2:97205/1‑139 (MQ=255)
                                                ggAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCaa                                                                                               >  1:613079/1‑123 (MQ=255)
                                                ggAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCaa                                                                                               <  2:613079/123‑1 (MQ=255)
                                                 gAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGttcttc                                                                                                          <  1:767876/111‑1 (MQ=255)
                                                 gAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGttcttc                                                                                                          >  2:767876/1‑111 (MQ=255)
                                                    aaCTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAggg                                                                                                                                <  1:1085011/86‑1 (MQ=255)
                                                    aaCTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCg                                                                           >  2:944453/1‑139 (MQ=255)
                                                    aaCTTTGACAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAggg                                                                                                                                >  2:1085011/1‑86 (MQ=255)
                                                         tGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTcgcg                                                                      >  1:1564395/1‑139 (MQ=255)
                                                         tGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTcgcg                                                                      >  1:1130792/1‑139 (MQ=255)
                                                          gTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAagcagc                                                                                                                                     <  1:1355414/75‑1 (MQ=255)
                                                          gTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAagcagc                                                                                                                                     >  2:1355414/1‑75 (MQ=255)
                                                               tgtCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAAc                                                                >  2:641678/1‑139 (MQ=255)
                                                                 tCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGt                                                              >  1:99013/1‑139 (MQ=255)
                                                                           gTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCgg                                                    >  2:293586/1‑139 (MQ=255)
                                                                                tgGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGttttt                                               <  2:1532143/139‑1 (MQ=255)
                                                                                  gAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCg                                             <  1:1295898/139‑1 (MQ=255)
                                                                                   aGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGc                                            >  2:1078665/1‑139 (MQ=255)
                                                                                     ttCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTg                                          >  1:165966/1‑139 (MQ=255)
                                                                                     ttCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTg                                          <  2:1130792/139‑1 (MQ=255)
                                                                                      tCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTcc                                                                            <  1:1292561/104‑1 (MQ=255)
                                                                                      tCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTcc                                                                            >  2:1292561/1‑104 (MQ=255)
                                                                                              ttCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGa                                 <  1:97205/139‑1 (MQ=255)
                                                                                               tCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAg                                <  1:641678/139‑1 (MQ=255)
                                                                                                ccAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGc                               <  2:1564395/139‑1 (MQ=255)
                                                                                                 cAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCa                              <  2:165966/139‑1 (MQ=255)
                                                                                                     gggTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGcc                          <  1:1437342/139‑1 (MQ=255)
                                                                                                      ggTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCa                         >  1:526582/1‑139 (MQ=255)
                                                                                                      ggTCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCa                         >  1:345036/1‑139 (MQ=255)
                                                                                                        tCGAGCAGGCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAAc                       >  2:171010/1‑139 (MQ=255)
                                                                                                               ggCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATggg                <  2:45650/139‑1 (MQ=255)
                                                                                                               ggCTGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATggg                <  2:99013/139‑1 (MQ=255)
                                                                                                                  tGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAg                                                                                                  <  1:1322102/54‑1 (MQ=255)
                                                                                                                  tGGTTTTGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAg                                                                                                  >  2:1322102/1‑54 (MQ=255)
                                                                                                                     ttttGTGGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTc          <  2:345036/139‑1 (MQ=255)
                                                                                                                          tgGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCgtg                                     >  2:621315/1‑107 (MQ=255)
                                                                                                                          tgGGAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCgtg                                     <  1:621315/107‑1 (MQ=255)
                                                                                                                            ggAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCa   >  2:1398209/1‑139 (MQ=255)
                                                                                                                             gAAGCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAg  <  2:1311142/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                       
CTTCGAAAAAAGCAGTTATCAGTCCCAGGTTTGAATTACCAGGATTAGGGAAAACTTTGTCAGTGTCACGAACCAGTTTGTGGAGTTCGCGAATTTCCAGCGGGTCGAGCAGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAG  >  NC_000913/1416756‑1417019

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