Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,328	0	C	A	48.3%	77.0 / 33.6	58	intergenic (+337/‑76)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (24/6);  new base A (0/28);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 4.79e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGC  >  NC_000913/345179‑345476                                                                                                                                                      |                                                                                                                                                     aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                                                     <  2:253202/151‑1 (MQ=255) aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                                                     <  1:343602/151‑1 (MQ=255)  aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaa                                                                                                                                                    <  1:761257/151‑1 (MQ=255)  aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                                                     >  1:302139/1‑150 (MQ=255)   aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaa                                                                                                                                                    >  1:707915/1‑150 (MQ=255)    aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaa                                                                                                                                                   >  2:164514/1‑150 (MQ=255)      aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAAc                                                                                                                                                >  1:1017919/1‑151 (MQ=255)           ctcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGcc                                                                                                                                            >  1:997226/1‑150 (MQ=255)              aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                                         >  2:110024/1‑150 (MQ=255)               tGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                                         >  1:11410/1‑149 (MQ=255)                gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                                                      >  1:810432/1‑151 (MQ=255)                  aGATGTTTTGTGAATGTGTTGGCTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGACTCATATCAACAACCGCCATTAGCtc                                                                                                                                    >  2:363823/1‑151 (MQ=255)                   gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc                                                                                                                                    >  2:891184/1‑150 (MQ=255)                                 gttttggttgACATTAATATAAAAAAAATCATAATCCCATCATGTTTGTGGGGGTTTTTTTCTGTTTCTTACGCACCGCCCCCAAGCCTAACAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAg                                                                                                                      <  2:193198/147‑1 (MQ=255)                                  tgttgtTTAATATTAATAAAAAAAAAATCATAATCCCATCTTGTTTGTGTGGGTTTTTTTTTGTATCTTCCGCATCGCACCCAAGCCTAAAAAAAAATATGCTGAAAGTCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAga                                                                                                                     <  2:190872/150‑1 (MQ=255)                                    ttgtttgACATTAAAAAAAAAAAAATCATAATCCCATCTTGTTTGTGGGGGTTTTTTTCTTTATCTTACCCACCGCACTAAAGCCTAAAAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                                    <  1:741842/149‑1 (MQ=255)                                       gttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                                >  1:768366/1‑150 (MQ=255)                                             aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                          >  1:90726/1‑150 (MQ=255)                                              ttaataaaaaaaaaaaCATATTCCCATCATGTTTGTGGGGGGTTTTTTCTTTATCTCCCCCACCGCACCCAAGCCAAAAAAAAAATAGGCTGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAg                                                                                                           <  2:181245/140‑1 (MQ=255)                                                   ataaaaaaaaCCAAAATCCCATCATTTTTGTGTGGGGTTTTTTCTTTATCTTACGCACCGCACCCAAGCCTGAAAAAAAATATGCTGAAGGGCTCAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                                       <  2:947013/144‑1 (MQ=255)                                                        aaaaaTAATAATCCCATCATTTTTGTGGGGGGTTTTTTCTTTATTTTCCGCACCGCACACAAGCCTGAAAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                              <  1:961845/151‑1 (MQ=255)                                                                 attcccATCATGTTTGTGGGGGTTTTTTTCTGTATCTTCCGCACCGCACTCAACCCTGAAAAAAAATATGTTGTAAGGTTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                     <  1:406662/147‑1 (MQ=255)                                                                 aaTCCCATCTTTTTTGTGGGGGTTTTTTTTTTTATTTTCCCCACCCCACCAAACCCTAAAAAAAAATATGCTGTAAGGCTAATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                                        <  2:284122/148‑1 (MQ=255)                                                                       atcagttttGGGGGGGGTTTTTTTTTTTTCTTCCGCACCCCACCAAAGCCTAAAAAAAAATATGTTGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                                                <  2:655686/144‑1 (MQ=255)                                                                           tGTTTGTGTGGGGTCTTTTCTTTATCTTCCGCACCGCACTCACGCCGGACAAAAAATATGCTGTAAGCCTCAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                            <  2:30150/150‑1 (MQ=255)                                                                             tttGTGTGGGTTTTTTTCTGTATCTTCCGCACCCCACCCAAGCCAAACAAAAAATATGCTGAAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                           <  2:43858/149‑1 (MQ=255)                                                                             tttGTGGGGGTTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                           <  1:886221/149‑1 (MQ=255)                                                                             tttGTGGGGGGTTTTTTCTGATTCTTACCCACCCCACTAAACCCAAAAAAAAATTAGGCTGTAGGGTTCATTAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTca                                                                           <  1:324261/149‑2 (MQ=255)                                                                                 tgtgGGTTTTTTTCTATATCTTCCGCCCCGCACAAAACCCAGAAAAAAAAAATTCTGTAGGGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                       <  2:464781/149‑1 (MQ=255)                                                                                    ggggTTTTTTCTGTATCTTCCGCACCCCCCTCAAGCCGGAAAAAAATTATGTGGTAAGGTTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                       <  2:464894/146‑1 (MQ=255)                                                                                     ggtttttttCTTTATCTTACGCATCGCACTCAAGCCAGAAAAAAAATATGCTGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcc                                                                 <  2:691264/146‑1 (MQ=255)                                                                                        tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                              >  1:660403/1‑151 (MQ=255)                                                                                         cTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                             >  2:709414/1‑151 (MQ=255)                                                                                          ttttttttATCTTACGCACCCCCCTCAAGCCTGACAAAAAATATGCTGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                             <  2:201165/150‑1 (MQ=255)                                                                                          ttttCTGAATCTTACGCACCCCACACAAGCCTGAAAAAAAATAGGCTGAAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                                <  2:1017870/147‑1 (MQ=255)                                                                                           tttCTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                            <  2:678049/150‑1 (MQ=255)                                                                                                gTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaa                                                       >  2:188475/1‑150 (MQ=255)                                                                                                     ttACGCATCGCACCCAACCTTGAAAAAAAATTTGTTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                                                   <  1:863386/149‑1 (MQ=255)                                                                                                      tACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                >  2:841249/1‑151 (MQ=255)                                                                                                         gCATCGCCCTCAACCCTGAAAAAAAAAATGCGGTAAGGCTAAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                  <  1:616894/146‑1 (MQ=255)                                                                                                           acccccctaaacccTGAAAAAAAATATCCTGAAGGACTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                            <  1:416313/138‑1 (MQ=255)                                                                                                           aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                            >  1:178564/1‑150 (MQ=255)                                                                                                              gCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                         >  1:796579/1‑150 (MQ=255)                                                                                                                aCTCAAGCCAGAAAAAAAATATGCTGTAAGCTTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                        <  2:268883/149‑1 (MQ=255)                                                                                                                 ctaaaGCCTGAAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                                      <  1:164524/147‑1 (MQ=255)                                                                                                                   cAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                   >  1:239221/1‑151 (MQ=255)                                                                                                                      gCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCt                                >  2:599456/1‑151 (MQ=255)                                                                                                                        cTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAACAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                               >  2:830151/1‑150 (MQ=255)                                                                                                                         tGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                             >  1:414802/1‑151 (MQ=255)                                                                                                                             aaaaaataattagTAAGTCTCTTATCAAATACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCATATCCCCttatt                            <  2:715703/136‑1 (MQ=255)                                                                                                                                  atatGCTGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATccc                    <  2:422868/151‑1 (MQ=255)                                                                                                                                   tatGCTGAAGGTCTAATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATccc                    <  2:733035/150‑1 (MQ=255)                                                                                                                                    atGCTGAAGGGCCCAAAGAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                   <  1:727624/150‑1 (MQ=255)                                                                                                                                         gTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa               >  1:703599/1‑149 (MQ=255)                                                                                                                                           aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt            <  2:1020936/150‑1 (MQ=255)                                                                                                                                           aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg           >  2:883629/1‑151 (MQ=255)                                                                                                                                           aaGGCTCAAATCAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg             <  1:145496/149‑1 (MQ=255)                                                                                                                                             ggCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg          <  1:978065/150‑1 (MQ=255)                                                                                                                                                      cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGc  >  2:56472/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGC  >  NC_000913/345179‑345476

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