Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,164	0	A	C	46.2%	75.7 / 22.0	52	T274P (ACC→CCC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (27/1);  new base C (0/24);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 5.86e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAA  >  NC_000913/3001020‑3001309                                                                                                                                                          |                                                                                                                                                   gTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGt                                                                                                                                               >  1:180058/1‑150 (MQ=255)     gAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑Accc                                                                                                                                          >  1:1068277/1‑150 (MQ=255)             aaaCCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑TCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCtt                                                                                                                                 >  1:63555/1‑151 (MQ=255)                  aTTTGTCGGTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTAcc                                                                                                                            >  1:378912/1‑151 (MQ=255)                         ggTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTTCTGACCAGCAAGCTTGGCGGCATT‑CCGGGG‑AAAGTTT‑TCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGAc                                                                                                                      >  2:965656/1‑150 (MQ=255)                         ggTGGCGGTTTTGAAAATAAA‑AGGATTGTGTGGGAGGACCAAAT‑GCCGGCTTTCCGACCCAGCATTTTTGGCGGC‑TT‑CCGGGT‑AAAGTTT‑CCC‑CT‑GTCCGGTG‑AGGTTTTTTCCCCCGAACCCCGT‑CCCCGCCACGCTTTTACCATTGACg                                                                                                                     <  1:743004/151‑1 (MQ=255)                           tGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGC‑TTTACCATTGACggg                                                                                                                   >  1:922023/1‑150 (MQ=255)                              cGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCa                                                                                                                 >  1:813177/1‑150 (MQ=255)                                 ttttGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATg                                                                                                             >  2:718531/1‑151 (MQ=255)                                  tttGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATgg                                                                                                            >  1:438722/1‑151 (MQ=255)                                    tGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                          >  1:1020040/1‑151 (MQ=255)                                     ggTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCg                                                                                                         >  2:142443/1‑151 (MQ=255)                                     gaaaaaaaaCGGGA‑TGTACGGGAAGACCCAAT‑GGCGGCTTTCCTGCCCAGCAAGTGGGGCGGCTTT‑CCGGTA‑AAATTTC‑CCC‑TT‑AGCC‑GTG‑AAGGGTTTTTCCCCCCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                          <  1:144655/144‑1 (MQ=255)                                       aaataaACAGGA‑TGTACGGGAAGAGACAAG‑GGGGGCTTTCCGGCCAAGCAGGCTGGGGGGCTTC‑CGGGTA‑AATGTTC‑CCC‑TT‑AGCG‑GTG‑AAGTGTTTTCCCCCGCACCCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                          <  1:179686/147‑1 (MQ=255)                                        aaaaaaCGGGA‑TGTACGGGAAGAGCAATT‑GGCGCTATCCCGGACCAGCAACGTTGG‑GGCATTGCCGGTA‑AAATTTC‑CCC‑TC‑ACCC‑GGG‑AGGTTTTTTCCCCCGCACCCGGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCg                                                                                                         <  1:441290/145‑1 (MQ=255)                                        aaaaaaCAGGA‑TGTACGGGAAGAGCAAAG‑GGCGGCATTCCGGACCAGCAAGCTTGGCGGCTTT‑CGGGTA‑AAATTTC‑CCC‑TT‑AGCC‑GTA‑AAGAGTTTTCCCCCGCACCCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGa                                                                                                      <  1:1126533/148‑1 (MQ=255)                                          taaACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAc                                                                                                    >  2:489532/1‑151 (MQ=255)                                             accaggtttttcTGGAAGAGCAAGG‑GGGGGCATTCCTGACCAGCAAGCTGGGCGGCATC‑CGGGAA‑AAATTTT‑CCT‑TG‑AGCC‑TGA‑AAGGGTTTTCCCTCGCACCCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcg                                                                                                  <  1:1034121/140‑1 (MQ=255)                                                ggA‑TGTACGGGAAGACCCAAG‑GGCGGCTTTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTA‑AATTTTC‑CCC‑TT‑AGCC‑GTG‑AAGAGTTTTTCCTCGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                 <  1:187249/148‑1 (MQ=255)                                                 gtttGTACGGAAAGACCCAGT‑GGCGCCA‑TCCCGACCCGGCAACTTGGGCGCCAT‑CCGGGT‑AAATTTT‑CCC‑TG‑TGCC‑GGGAAAGAGTTTTCCCCCCCACCCCAA‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                 <  2:1097899/145‑1 (MQ=255)                                                      tACTGGAAGAGCAAAG‑GGCGGCTTTCCGGACCAGCAAGCTTGGCGGCATT‑CCGGTA‑AAATTTC‑CCC‑TT‑AGCC‑GTG‑AAGGGTTTTTCCCCGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa                                                                                          <  1:884312/150‑1 (MQ=255)                                                       aCTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                         >  1:27894/1‑150 (MQ=255)                                                         tGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                      >  1:249186/1‑151 (MQ=255)                                                            aaGAGCAAAG‑GGCGGTATTCCTCACCAGCAAGTTTGGGGGCATT‑CGGGTG‑AAAGTTC‑CCC‑TT‑ACCC‑GTA‑AAGAGTTTTCCCCCGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGa                                                                                    <  1:871244/150‑1 (MQ=255)                                                             aaggCCAATGGGCGCCATCCCTGACCAGCAGGTTGGGCGGCTTC‑CGGGTA‑AAATTTC‑CCC‑TT‑CGCC‑GTA‑AAGAGTGTTCCCCCCCACCCCTT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACa                                                                                        <  1:411798/143‑1 (MQ=255)                                                                 gccaggggcgccATCCCGGACCAGAAGGTTTGGCGCCTTT‑CCGGTA‑AAATTTC‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCCCGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                                                                                  <  2:579194/142‑1 (MQ=255)                                                                    aagggcggcTTCCCTGACAAGAAAGTTTGGGGGCTTT‑CGGGTA‑AAATTTC‑CCT‑TT‑AGCG‑GTA‑AAGAGTTTTCCCCCGAACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                            <  1:779240/148‑1 (MQ=255)                                                                     t‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                            >  1:583972/1‑149 (MQ=255)                                                                             aTTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGAt                                                                   >  2:1017655/1‑151 (MQ=255)                                                                                 cTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc                                                               >  2:526867/1‑151 (MQ=255)                                                                                  tGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                              >  1:295689/1‑151 (MQ=255)                                                                                    aCCAGCAGGCTTGGCGTCTTC‑CGGGTG‑AATTTTC‑CCT‑TT‑AGCC‑GGA‑AAGTGTTTTTCCCCGAACCCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                              <  2:1101868/149‑1 (MQ=255)                                                                                       agaaaGCTGGGGGGCATC‑CCGGTG‑AATTTTC‑CCC‑TG‑AGCC‑GTG‑AAGAGTGTTTCCTCGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                              <  2:929307/143‑1 (MQ=255)                                                                                                  ggcggcctgcGGTGAAAAGTTC‑CCC‑TT‑TTCC‑GTG‑AAG‑TTTTTTCCCCGACACCCGTAACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgc                                                 <  1:657803/141‑1 (MQ=255)                                                                                                   ggcgcgattcgggtaaaTTTT‑CCCTTT‑CGCC‑GTG‑AAGGGTTTTTCCCCGCACCCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCgg                                                  <  2:912845/135‑1 (MQ=255)                                                                                                     cATT‑CCGGTG‑AATTTTT‑CCC‑TT‑GGCC‑GTG‑AAGAGTGTTTCCCCGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTAt                                           <  2:1107133/151‑1 (MQ=255)                                                                                                             ggaaaatgTTC‑CCC‑TT‑AGCC‑GTG‑AGGATTTTTCCCCCGCACCCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGc                                         <  2:1130554/140‑1 (MQ=255)                                                                                                               gaaaaaTTT‑CCC‑TTGCGCC‑GAG‑AAGGTTTTTCCCCCCCACCCCGA‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCa                                        <  2:502056/145‑1 (MQ=255)                                                                                                                 aaaGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                 >  2:326499/1‑151 (MQ=255)                                                                                                                  aaGTTTCCCC‑TT‑GGCC‑GGG‑AAGGGTGTTCCCCCCCACCCCGC‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCa                                   <  1:1072780/149‑1 (MQ=255)                                                                                                                  aaGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACgg                                >  2:1035388/1‑151 (MQ=255)                                                                                                                    gTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACggg                               >  2:586424/1‑150 (MQ=255)                                                                                                                              t‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                      >  2:458626/1‑151 (MQ=255)                                                                                                                                  cc‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCt                   >  1:198262/1‑151 (MQ=255)                                                                                                                                  cc‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCt                   >  2:665491/1‑151 (MQ=255)                                                                                                                                         aaggtttttcccccGCACCCCTT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTCCTg               <  1:538515/137‑1 (MQ=255)                                                                                                                                              ttttcccgcGCACCCCTA‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc              <  2:1015718/138‑1 (MQ=255)                                                                                                                                                  tCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGa      >  2:708179/1‑150 (MQ=255)                                                                                                                                                   ccccGCACCCCGT‑CCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggg       <  1:748122/145‑1 (MQ=255)                                                                                                                                                   ccccGCACCCAGT‑GCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg        <  2:267490/144‑1 (MQ=255)                                                                                                                                                      cGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa  >  1:642439/1‑150 (MQ=255)                                                                                                                                                       gCAACACGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa  >  2:389988/1‑149 (MQ=255)                                                                                                                                                          |                                                                                                                                                   GTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGA‑TGTACTGGAAGAGCCAAT‑GGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTG‑AAAGTTT‑CCC‑TT‑AGCC‑GTG‑AAGAGTGTTTCCTCGCAACCCGT‑ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAA  >  NC_000913/3001020‑3001309

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