Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3950CA40.0% 89.8 / 32.8 56S131R (AGC→AGAydcRputative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (28/5);  new base A (0/22);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA  >  NC_000913/1510249‑1510542
                                                                                                                                                     |                                                                                                                                                    
cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa                                                                                                                                                  <  1:99567/151‑1 (MQ=255)
 gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAAc                                                                                                                                                >  2:469308/1‑151 (MQ=255)
 gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAAc                                                                                                                                                <  2:953835/151‑1 (MQ=255)
  agaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAAcc                                                                                                                                               >  1:707828/1‑151 (MQ=255)
    aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGc                                                                                                                                             <  1:305398/151‑1 (MQ=255)
       aaGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCcac                                                                                                                                          >  1:680238/1‑151 (MQ=255)
         gCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCcac                                                                                                                                          >  2:757277/1‑149 (MQ=255)
           aGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCg                                                                                                                                      >  1:38617/1‑151 (MQ=255)
                  atCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCg                                                                                                                               >  1:740570/1‑151 (MQ=255)
                    cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTg                                                                                                                             >  2:855929/1‑151 (MQ=255)
                    cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTCCCCGACCCGCGACTTTCCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGAAAAGCCGCCACAGCGATGAGCg                                                                                                                               <  1:137881/149‑1 (MQ=255)
                        acTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTg                                                                                                                         >  2:130434/1‑151 (MQ=255)
                          ttttttttttGATTTGTTGCGGCCCAGTCGCGACCCGTGGGTCTTCCCTTTTCCCCCGCCTTTCCCCGCCCCGCTATTTTCCCCCCCCCAAAAAAAAACCGGCTCGCGGGCGGGGGTAAAAAAAACCGCCACAGCGATGAGCGt                                                                                                                              <  1:601257/139‑1 (MQ=255)
                          ttATTTTTTTGATAGGTGGCAGCCCGGCCGCGCCCCGTGTGTCGTCCGTTTTCCCCCGCCTTTCCCCGCCCCGCGATTTCCCCCCCCCCAACAACAAACCCGCTCGCTGGCGAAGGAAAGAAAAACCGCCACAGCGATGAg                                                                                                                                 <  1:886062/141‑1 (MQ=255)
                          ttATATTTTTAATATGTTGCAGGCCAGTCGCGATCCGTCGGTCTTTCCTTTTCCCCCGGCCTTCCCCGCCCCGCGATTTTCCCCCCCCCAACAACAAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTg                                                                                                                         <  2:247929/149‑1 (MQ=255)
                             tttttttGATTTGTTGCGGGCCGGTCGCGACCCGTCGGTCTTTCCTTTTGCCTCGGCTTTCCCCGCCCCGCGCTTTTCCCCCCCCCAACAACAAACCCGCTCGTTGGCGCAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTg                                                                                                                         <  1:659192/144‑1 (MQ=255)
                               tttttGATTTGTGGCAGGCCAGGCGCGCCCCGCGGTTCGTTCGGTTTCCCGCGCCTTTCCCCGACCCGCGCTTTCCCCCCCCCAAACAACAAACCCGCTCGCTGGCGCAGGAAACCAAATCCGCCACAGCGATGAGCGTGATTGaaaa                                                                                                                     <  2:613704/148‑1 (MQ=255)
                               ttttgGATTTGTGGCGGGCCAGGCGCGTCCCGTCGTTCGTCCCGTTGCCCCCGCCCTTCCCCGCCCCGCGACTTTCCCCCCTCCAACAACAAACCCGCTCGCGGGCGCGGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGaaaa                                                                                                                     <  2:382237/148‑1 (MQ=255)
                                ttttGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACtt                                                                                                                  >  2:943703/1‑150 (MQ=255)
                                ttttGATATGTTGCAGCCCGGCCGCGATCCGGGGGTCTTTCTTTTTCCCCGGCCCTTCCCCGCCCCGCCTCTTTCCCCCCCCCAAAAAAAAACCCCCCCGCTGGCGCAGGTAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACtt                                                                                                                  <  1:316134/150‑1 (MQ=255)
                                       atGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc                                                                                                            >  2:782535/1‑149 (MQ=255)
                                        tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                                         >  2:822317/1‑151 (MQ=255)
                                        tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGAGTGATTGAAAACTTACCGCCAgg                                                                                                         >  1:850575/1‑151 (MQ=255)
                                           tGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaa                                                                                                       >  1:263547/1‑150 (MQ=255)
                                              aGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                                   >  2:678032/1‑151 (MQ=255)
                                               ggccggccGCGACCCGCGGGTCTTCCCTTTTCCCCCGCCTTTCCCCGCCCCGCGTTTTTCCCCCCCCAAAAAACAAACCCGCCGGCGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg                                                                                                    <  2:21365/149‑1 (MQ=255)
                                                 ccacccgcgccccGTCGTTCGTCCCTTTGCCCTCGGCTTTCCCCGCCCCGCGTTTTTCCCCCCCCCACCAACAAACCCGCTCGCTGGCGGAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                                   <  2:241028/137‑1 (MQ=255)
                                                   aGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAAc                                                                                              >  1:680196/1‑151 (MQ=255)
                                                       gcgATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCg                                                                                          >  1:715351/1‑151 (MQ=255)
                                                       gccacccGTCGGTCTTTCGTTTGGCCTCGGCTTTCCCCCCCCCCCGCTTTTCCCCCCCCCAACAACAAACCCCCTCGCGGGCGCGGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGc                                                                                           <  1:830892/145‑1 (MQ=255)
                                                        gcgccccttggttcttCCCTTTGCCCCCGCCTTTCCCCGCCCCGCGTTTTTCCCCCCCCCAACAACAAACCCGCTCGCGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt                                                                                             <  2:409808/134‑1 (MQ=255)
                                                         gATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc                                                                                        >  1:219915/1‑151 (MQ=255)
                                                          gtcccttcgttcttCCCTTTGCCCCCGCCTTTCCCCACCCCGCGCTTTTCCCCCCCCCAACAACAAACCCCCCCGCGGGCGCAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt                                                                                             <  1:761777/134‑1 (MQ=255)
                                                            ccGTGGGCCTTTCCGTTTCCCCCGCCTTTCCCCGCCCCCCGCTTTTCCCCCCCCAAACAACAAACCCCCCCGGGGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTg                                                                                            <  1:993116/144‑1 (MQ=255)
                                                              gtcggtcgTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt                                                                                   >  1:743935/1‑151 (MQ=255)
                                                                 ggtcgTTCCGTTTGCCTCGGCCTTTCCTGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAtt                                                                                >  1:761692/1‑151 (MQ=255)
                                                                         cctttcgccgcgCCTTTCCCCGCCCCGCGCCTTTCCCCCCCCAAAAAACAAACCCGCTCGCGGGCGCAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTa                                                                                  <  2:671055/131‑1 (MQ=255)
                                                                                   gcccTTCCCCGCCCCGCGCTTTTCCCCCCCCAAAAAACAAACCCGCCCGGGGGCGAGGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCt                                                                <  2:97755/147‑1 (MQ=255)
                                                                                      ctttccccgccccgcGTTTTTCCCCCCCCCAAAAACAAACCCGCTCGCGGGCGGAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGcc                                                                 <  1:831445/141‑1 (MQ=255)
                                                                                      ctttccccgccccgcGCTTTTCCCCCCCCCAACAACAAACCCCCTCGCGGGCGCAGAAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt                                                               <  1:757156/143‑1 (MQ=255)
                                                                                        ttccccgccccgcGTTTTTCCCCCCCCAAAAAACAAACCCGCTCGCGGGCGGGGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa                                                            <  2:606687/146‑1 (MQ=255)
                                                                                                 cccGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc                                                  >  2:368598/1‑151 (MQ=255)
                                                                                                     cgACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAttt                                              >  1:245437/1‑151 (MQ=255)
                                                                                                        cTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc                                           >  1:326265/1‑151 (MQ=255)
                                                                                                         tttcCCCCCCCAAACAACAAACCCCCCCGCGGGGGCGGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAt                                                <  2:126953/145‑1 (MQ=255)
                                                                                                                cccccAACAACAAACCCCCCCGCTGCCGAGGGAAACAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGaa                                   <  2:931363/148‑1 (MQ=255)
                                                                                                                ccTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga                                    >  1:175967/1‑150 (MQ=255)
                                                                                                                    caacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg                               >  1:439043/1‑151 (MQ=255)
                                                                                                                             aacccccccgctGGCGAAGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt                       <  1:782414/143‑1 (MQ=255)
                                                                                                                              aaCCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTg                      >  1:973437/1‑151 (MQ=255)
                                                                                                                                 ccctcgggggCGAAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg                   <  2:570332/143‑1 (MQ=255)
                                                                                                                                     ccctGCCGCAGGTCAATAAATCCGCCACAGCGATGAGCGTGATTGAAATCTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc               <  1:382226/149‑1 (MQ=255)
                                                                                                                                         ggCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt            >  1:303650/1‑150 (MQ=255)
                                                                                                                                         ggCGAAGGTAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt            <  2:1069621/150‑1 (MQ=255)
                                                                                                                                            gcAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTag          >  2:144598/1‑149 (MQ=255)
                                                                                                                                                 aaaacAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCAtt   <  2:884168/147‑1 (MQ=255)
                                                                                                                                                   aGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTa  >  2:377316/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCCAACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA  >  NC_000913/1510249‑1510542

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