Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	44.6%	66.4 / 21.3	56	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (23/8);  new base C (0/25);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 4.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556057‑3556345                                                                                                                                                         |                                                                                                                                           gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACt                                                                                                                                            <  1:677640/151‑1 (MQ=255)   gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTcc                                                                                                                                          <  1:744470/151‑1 (MQ=255)       cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGc                                                                                                                                      <  1:159013/151‑1 (MQ=255)        aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGcc                                                                                                                                     >  2:284919/1‑151 (MQ=255)                           cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                  >  2:973214/1‑151 (MQ=255)                           cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                  >  2:937789/1‑151 (MQ=255)                               tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                              >  1:547348/1‑151 (MQ=255)                                 tccccaccgaaAAAAAGCCGTTGGGTCATATTTTCCCGTCCGACTCCAAGCGAGCCGTTATTACCCGCCCCC‑GGGCCGCGCGGCGGTTTATGAATTTTCG‑TTTCATGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                              <  1:761265/142‑1 (MQ=255)                                  cACCGACAAAAAAAAAGCTTTGGGTAATATTTTCCCTTCCGACTTCAAGCGAGACGTTATTCCCCGCCCCC‑GGGCCGCGCGGCAGTTTATAAATTTTCT‑TTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                            <  2:358814/150‑1 (MQ=255)                                        ccaaaaaaaaCCTTTGGGGGATATTTTCCCTTCCGCCTCAAGGCGAAGCGGTTTTCCCCGGCCCG‑TGGTGCCGGGGCAGGTTA‑GAATTTTTTGGTGCATGCGGAGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  1:684083/145‑1 (MQ=255)                                           cAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                  >  1:1019160/1‑151 (MQ=255)                                             aaaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                 >  1:905169/1‑150 (MQ=255)                                                  gCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                            >  2:898057/1‑150 (MQ=255)                                                   cGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                            >  2:99538/1‑149 (MQ=255)                                                    gTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                         >  1:651304/1‑151 (MQ=255)                                                     ttgggTAATATTTTCACTTTCGCCTCCAAGCGAGACGTTATTCCCCGGCCCC‑TGGTCGCGGGGCAGTTTATGAATATTCT‑GTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                          <  2:635694/146‑1 (MQ=255)                                                     ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                        >  2:303602/1‑151 (MQ=255)                                                       cGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                       >  2:897420/1‑150 (MQ=255)                                                        gttGATATTTTCTCTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCC‑GGGTCGCGCGGCAGTTTATGAATATTCT‑GTTCATGCGGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                       <  2:773286/147‑1 (MQ=255)                                                        ggTGATTTTTTCCCTTCCGACTTCAAGCGAGACGGTATTCCCCGGCCCC‑TGGTCGCGCGGCAGTTTATGAATTTTCT‑GTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                         <  2:13122/147‑1 (MQ=255)                                                         gtatttttttCACTTTCGACTTCAGGCAAGCCGGTTTCCCCCGGCCCC‑GGGTCGCGCGGAGGGTTATAAATATTCT‑GTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                      <  2:869424/143‑1 (MQ=255)                                                             ttttttCTCTTTCGCCTTCAAGCGAGACGTTATTCCCCGGCCCC‑TGGTCGCGCGGCAGTTTATGAATTTTCT‑GTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                   <  2:888057/146‑1 (MQ=255)                                                             ttttttCACTTGCGACTACAAGCGAGACGGTATTCCCCGCCCCC‑GGGCCGCGGGGCGGGTTATAAATATTCT‑TTTCAGGGAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                  <  1:238735/147‑1 (MQ=255)                                                               ttttCCCTTTCGACTACAAGCGAGACGTTATCCCCCGGCCCC‑GGGCCGCGCGGCGGGTTATAAATATTTT‑TTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                <  2:603108/149‑1 (MQ=255)                                                               ttttCACTTTCGACTCCAAGCGAGACGGTATTCCCCGGCCCC‑TGGTCGCGCGGCAGTTTATGAATATTCT‑GTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                <  2:518827/149‑1 (MQ=255)                                                                tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                             >  2:512200/1‑151 (MQ=255)                                                                 ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                            >  2:350965/1‑151 (MQ=255)                                                                      tttCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       >  2:562181/1‑151 (MQ=255)                                                                       ttCGTCTACAAGCGAGTCGTTTTTCCCCGCCCCC‑TGGCCGCGCGGGAGGTTATGAATTTTCT‑TTTAATGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       <  1:831033/150‑1 (MQ=255)                                                                         cGACTTCAAGCGAGACGTTA‑TCCCCGGCCCCTTGGTCGCGCGGCAGTTTATAATTTTTTT‑TTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       <  1:797829/148‑1 (MQ=255)                                                                         cGACTCCAAGCGAGACGGTATTCCCCGCCCCC‑GGGCCGCGCGGCAGGTTATGAATTTTCT‑GTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                      <  2:510241/149‑1 (MQ=255)                                                                               cAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                              >  1:101399/1‑151 (MQ=255)                                                                                aaGCGAGACGGTATTCCCCGCCCCG‑TGGGCGCGCGGCAGTTTATGAATATTCT‑GTTCATGAAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                             <  1:167796/151‑1 (MQ=255)                                                                                    gagaCGTTATTCCCCGCCCCC‑TGTCCGCGCGGCAGTTTATGAATTTTTT‑TTTCATGAAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                           <  1:445317/149‑1 (MQ=255)                                                                                    gagaCGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                         >  2:464192/1‑151 (MQ=255)                                                                                    gagaCGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                         >  1:653045/1‑151 (MQ=255)                                                                                           tATTCCCCGCCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                    <  1:336425/149‑1 (MQ=255)                                                                                              tACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                               >  1:538304/1‑151 (MQ=255)                                                                                               tcccggcacccggcgcgcgcgGGCGGTTATTACTTTTCT‑GTTCAGGAAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                                <  2:72982/135‑1 (MQ=255)                                                                                               ccccGCCCCC‑TGGTCGCGGGGCATTTTATAAATTTTCT‑TTTCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                              <  1:1024485/150‑1 (MQ=255)                                                                                                cccGGCCCC‑TGGTCGCGGGGCAGGTTTTGAATTTTCT‑GTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                               <  1:299739/149‑1 (MQ=255)                                                                                                cccGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                             >  2:1007204/1‑151 (MQ=255)                                                                                                cccGGCCCC‑TGGGCGCGCGGCAGTTTATGAATATTCT‑GTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                               <  1:978486/149‑1 (MQ=255)                                                                                                     cccc‑GGGTCGCGCGGCAGTTTATGAATATTCT‑GTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCa                                           <  1:388357/148‑1 (MQ=255)                                                                                                        c‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                     >  1:552778/1‑151 (MQ=255)                                                                                                              cgggcgcaagtttATGAATTTTCT‑TTTCGTGCAGAAAAACCCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCaa                                          <  1:31664/130‑1 (MQ=255)                                                                                                               gcgcgGCAGTTTATAAATTTCCT‑GTCCATGGGGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                                <  1:62794/150‑1 (MQ=255)                                                                                                               gcgcgGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                               >  1:130083/1‑151 (MQ=255)                                                                                                                     cAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                         >  1:812972/1‑151 (MQ=255)                                                                                                                         ttATGAATATTCT‑GTTCATGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      <  2:424323/150‑1 (MQ=255)                                                                                                                             aaaTATTTT‑GTGCATAAAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       <  1:973193/144‑1 (MQ=255)                                                                                                                                tgttCT‑TTCCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg               <  1:169537/148‑1 (MQ=255)                                                                                                                                  ttCT‑GTTCATGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             <  2:131143/150‑1 (MQ=255)                                                                                                                                        ttCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag          <  2:1015783/148‑1 (MQ=255)                                                                                                                                         tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt      >  1:474360/1‑151 (MQ=255)                                                                                                                                            ggCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg            <  1:354114/141‑1 (MQ=255)                                                                                                                                                  aaaaacCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  1:917553/145‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                           GTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCC‑TGGTCGCGCGGCAGGTTATGAATATTCT‑GTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556057‑3556345

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: ‑