Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG50.0% 48.0 / 29.2 50T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (5/20);  new base G (25/0);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 2.25e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366379‑366675
                                                                                                                                                     |                                                                                                                                                      
aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAggggggg                                                                                                                                                     >  2:328146/1‑150 (MQ=255)
 gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTTGGGGCCCgggggggt                                                                                                                                                    >  1:954399/1‑149 (MQ=255)
 gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCGCCCgggggggt                                                                                                                                                    >  1:687568/1‑149 (MQ=255)
          gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTATTGGGCGGCCGGGGGGGttttcttttt                                                                                                                                           >  2:1196149/1‑150 (MQ=255)
           ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGTTTTGGGGCCCGGGGGGGtttttttttt                                                                                                                                            >  1:75068/1‑150 (MQ=255)
           ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGCTTTGGGGGCCGGGGGGGtttttttttta                                                                                                                                           >  1:1261256/1‑150 (MQ=255)
             tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGCTTTGGCGGCCGGGGGGGttttctttttca                                                                                                                                         >  1:333850/1‑147 (MQ=255)
                tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTGTTGGGGGCCGGGGGGGGTTTTCTTTTCACCAgg                                                                                                                                     >  1:437915/1‑150 (MQ=255)
                  agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                                    <  1:99953/150‑1 (MQ=255)
                   gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGGGGTTTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTgag                                                                                                                                  >  1:795688/1‑151 (MQ=255)
                    tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                 <  2:1017862/151‑1 (MQ=255)
                     gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCCGCTTGGGGGCGACGGGTGGTTTTCTTTTTAccccagagaa                                                                                                                                >  2:994163/1‑143 (MQ=255)
                       gCTAACTCCCATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGAGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAACCCGTGAGACg                                                                                                                               <  2:1003694/150‑1 (MQ=255)
                       gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTTTTGGGGGCCAGGGGGGGTTTTTTTTTCTcccgggggaggg                                                                                                                              >  1:772209/1‑141 (MQ=255)
                          aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGGCTTGTGCGCCACGGTGGGTTTTCTTTTTACCACTGCGGCgggg                                                                                                                            >  2:489003/1‑149 (MQ=255)
                                   tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGTTTTTTCTTTTACCAGGGGGACGGGCAACAGCTGAt                                                                                                                 >  2:556512/1‑151 (MQ=255)
                                   tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCGCCAGGGGGGTTTTTCTTTTCACCAGGGAGACGGGCAACAGCTGAt                                                                                                                 >  2:264559/1‑151 (MQ=255)
                                        gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTTCCAGTGAGAAGGGCAACAGCCTATTGcc                                                                                                             >  2:543405/1‑150 (MQ=255)
                                        gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTGTTGGGGGCCCGGGGGGGTTTTTTTTTCACCCGTGAGACGGGGAACAGGTGGTTggc                                                                                                             >  2:740108/1‑148 (MQ=255)
                                             gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGGGAGGCGGGCCACAGCTGATTGCCCCTct                                                                                                        >  1:805807/1‑149 (MQ=255)
                                              cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                                                                                                       <  2:638416/150‑1 (MQ=255)
                                                cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTATTTGGGGCCAGGGGGGGTTTTTTTTTTCCCCGGGGGACGGGCAAAAGGTGGGTGccccttcccc                                                                                                     >  1:742634/1‑143 (MQ=255)
                                                  cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTTTTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGGAACAGCTGGTTGGCCTTCACcccc                                                                                                   >  1:1122833/1‑150 (MQ=255)
                                                           cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                         <  2:960349/151‑1 (MQ=255)
                                                           cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGGTTTTCTTTTTCCCAGTGGGAAGGGCAACAGCTGGTTGCCCCTTCCCGCCCGGccccggg                                                                                         >  1:666734/1‑147 (MQ=255)
                                                                   tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt                                                                                    <  2:864037/148‑1 (MQ=255)
                                                                        aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGGGAGACCGGCCACAGATGAGTGCCCCTCACCGCCTGGCCCCGAGGGAGGTGCCGCaa                                                                            >  2:802497/1‑151 (MQ=255)
                                                                           ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGGATTTGGGGCCAGGGGGGGTTTTCTTTTCACCACGGGGGAGGGCGAAAGACGGATTGCCCTCCCCGCCCTGCCCCGTGAGAGTTTCGGCCAGCg                                                                         >  1:414624/1‑151 (MQ=255)
                                                                                          gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTGGGGGCCAGGGGGGGTTTTCTTTTCACCCGGGAGAAGGGCCAAAGCTGATTGGCCCTCACCGCCTTTCCCCGTGGGAGGTGGAGCACGGGGTCCCCCCTGGGTTg                                                          >  1:873520/1‑151 (MQ=255)
                                                                                          gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTGGGCGCCCGGGGGGTTTTTCTTTTTACCAGTGGGACGGGCAACAGCTGATTGCCCCTCCCCCCCTGGCCCTGAGGGAGTTGCAGCAAACGGGCCACCCCGGGTTg                                                          >  1:725620/1‑151 (MQ=255)
                                                                                          gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGGTTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGGGAAGGGCAAACGCTGATTGTCCTTCCCCCCCCGGGCCCGGGGGAGGTGCAGCCAACGGTCCACGCTGGTTTg                                                          >  1:278889/1‑151 (MQ=255)
                                                                                          gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCAGGGGGGTTTTTTTTTTCACCAGGGAGACGGGGAACAGCTGATTTCCCCTCCCCCCCCGGCCCTGGGGGGGGTGGAGCAAGCGGGCCAccccgggttt                                                          >  1:735259/1‑142 (MQ=255)
                                                                                          gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGGGCCAGGGGGGGTTTTTTTTTTACCCGGGAGACGGGGAAAAGCTGAATGGCCTTCACCCCCCGGCCCTGCGGGAGGTGCAGCAAGCGGGCCCCGCTGGGTTg                                                          >  1:834801/1‑151 (MQ=255)
                                                                                             ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTGGGGGGCCGGGGGGGTTTTTTTTTCTCCACCGGGACGGGGAACAGCTGGTTGCCCTTCCCCCCCCTGCCCCCAGAGGGGTGCAGCAAGCCGCCCCCGCCGGGTTGcc                                                        >  2:893192/1‑150 (MQ=255)
                                                                                                        gCCAACGCGCGGGGAGAGGGGGTTGGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                             <  1:781186/150‑1 (MQ=255)
                                                                                                        gCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGGGAGACGGGCAACAGCTGGTTGCCCTTCCCCGCCCGGCCCTGAGGGAGTTGCAGCAAACGGGCCACGCTGGTTTGCCCCACCAGGCgga                                            >  1:189500/1‑149 (MQ=255)
                                                                                                                 cgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                   <  2:1122407/151‑1 (MQ=255)
                                                                                                                   gggAGAGGCGGTTTTTGTATTTGGCGCCCGGGGGGGTTTTTTTTTCACCCGTGAGAAGGGGAAAAGCTGATTGCCCCTTACCGCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGTCCCGCTGGTTTGCCCCAGCAGGGGAAAATTCTGtttt                                 >  2:865771/1‑150 (MQ=255)
                                                                                                                        agGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                             <  2:771893/150‑1 (MQ=255)
                                                                                                                             gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                       >  2:903264/1‑151 (MQ=255)
                                                                                                                             gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                       >  1:485810/1‑151 (MQ=255)
                                                                                                                               ttGCGTATTGGGCGCAAGGGTGGTTTTTCTTTTCTCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                     <  2:950736/151‑1 (MQ=255)
                                                                                                                                  cGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                   <  1:829218/150‑1 (MQ=255)
                                                                                                                                     tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                >  1:406436/1‑151 (MQ=255)
                                                                                                                                           gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat           <  1:828688/150‑1 (MQ=255)
                                                                                                                                           gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat           <  2:518352/150‑1 (MQ=255)
                                                                                                                                           gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat           <  2:42403/150‑1 (MQ=255)
                                                                                                                                              ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTTACCGCCTGGCACTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATAAAAc        <  2:1164990/150‑1 (MQ=255)
                                                                                                                                              ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCCCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc        <  2:687259/150‑1 (MQ=255)
                                                                                                                                              ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc        <  1:1196643/150‑1 (MQ=255)
                                                                                                                                              ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa       <  2:954050/151‑1 (MQ=255)
                                                                                                                                              ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa       <  1:729202/151‑1 (MQ=255)
                                                                                                                                              ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa       <  1:676868/151‑1 (MQ=255)
                                                                                                                                                   ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt      <  2:156845/148‑1 (MQ=255)
                                                                                                                                                     tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  <  2:660553/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366379‑366675

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