Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	44.4%	67.6 / 15.2	54	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (10/20);  new base G (24/0);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 9.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508508‑3508806                                                                                                                                                       |                                                                                                                                                      cAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGaa                                                                                                                                                      <  1:332796/151‑1 (MQ=255)     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAAt                                                                                                                                                  <  2:215512/151‑1 (MQ=255)      cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGC‑GGGAGGgaaaa                                                                                                                                                   >  2:266406/1‑149 (MQ=255)            gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGAAAAAATcacac                                                                                                                                            >  1:860253/1‑150 (MQ=255)                ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGT‑GCGCGGGAAAACTcaccatc                                                                                                                                          >  1:736600/1‑144 (MQ=255)                 tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGGCCGGGAAAAAttac                                                                                                                                               >  1:791402/1‑141 (MQ=255)                  aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTa                                                                                                                                      <  1:346311/150‑1 (MQ=255)                   aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAg                                                                                                                                     >  1:259765/1‑150 (MQ=255)                         gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCcac                                                                                                                              <  1:437006/151‑1 (MQ=255)                           gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                            <  1:429686/151‑1 (MQ=255)                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAAATCACACACGGTAGCCCCcac                                                                                                                              >  1:546771/1‑147 (MQ=255)                                      gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                 <  1:32004/151‑1 (MQ=255)                                              tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCAAACACGGGAGCCCCCACATCAGAAAGCGTGGGgcgccc                                                                                                         >  1:702503/1‑149 (MQ=255)                                              tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCCGCGAAAAAAAAACACACCGGAGACCCCCCATCAAAAAGCgggggggg                                                                                                            >  2:869994/1‑144 (MQ=255)                                              tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGG‑GGGCGGAAAAAATCACACACGGGGGACCCCACATAAGAAAACGGGGGGCGca                                                                                                          >  1:28009/1‑150 (MQ=255)                                               cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAATTCACACACAGGAGCCCCCACATCAGAAAGAGTGGGGGGca                                                                                                          >  2:631039/1‑149 (MQ=255)                                               cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGcac                                                                                                         >  2:400112/1‑150 (MQ=255)                                               cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGca                                                                                                          <  2:364310/149‑1 (MQ=255)                                                       aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCACAcc                                                                                                      >  1:866877/1‑145 (MQ=255)                                                       aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCAAACACAGGAGACCCCACATCAGAAAGAGGGGGGGGCGCAcc                                                                                                      >  2:558125/1‑145 (MQ=255)                                                        gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                               <  2:572593/151‑1 (MQ=255)                                                             tAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GGGCGGGAAAAATCACACACGGTAGGCCCCCCATCAGAAAGCGGGGGGCGCACAcc                                                                                                      >  1:93083/1‑139 (MQ=255)                                                             tAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGG‑GCGCGGGAAAAATTACACCCCGGAGGCCCCCCATCCGAAAACGGGGTGCGCCCACCGGGCcgcg                                                                                              >  2:784342/1‑147 (MQ=255)                                                                    tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGatca                                                                                     >  1:171791/1‑149 (MQ=255)                                                                       gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACAACCAGCTGCGGAGCGATCATccc                                                                                 >  2:352215/1‑150 (MQ=255)                                                                          cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                              <  1:689798/150‑1 (MQ=255)                                                                            aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATTAAACACGGTAGCCCCCACATCAGAAAGCGGGGGGCGCCCACCCAGCTGCGTAGCGATCATcccct                                                                               >  1:146054/1‑145 (MQ=255)                                                                                tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGC‑GCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGAGTGGGGGGCACACCGAGCTGCGGAGCGGTCAACCCctgctggt                                                                         >  2:255574/1‑147 (MQ=255)                                                                                    aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                    <  1:151576/150‑1 (MQ=255)                                                                                     gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAg                                                                  <  1:560079/151‑1 (MQ=255)                                                                                         gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                               <  2:449000/150‑1 (MQ=255)                                                                                            tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAAAGGGGGGGGGACACCGAGCCGCGTAGCGATCATCCCTTGCTGATAAACCAGCAccg                                                             >  2:678111/1‑147 (MQ=255)                                                                                            tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAACACTgg                                                            >  1:810865/1‑150 (MQ=255)                                                                                               tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGGGGGAAAAATTACAAAAGGGAGCCCCCCCAACCAAAAGCGTGGTGGGGCCACCCCGCCGCGGAGCGGTCCTCCCCTTCTGGTAAAACCGCAcctgcccc                                                        >  2:73626/1‑144 (MQ=255)                                                                                                gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGC‑GCGCGGGAAAAATCACACACGGGAGCCCCCCCCACAGAAAACGGGGGGCGCCCACCGAGCCGCGGAGCGATCATCCCCTGCTTATAAACCAacacccggacc                                                        >  2:81356/1‑143 (MQ=255)                                                                                                  aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGC‑GCGCGGGAAAAATCACACACGGTAACCCCCACATCAGAAAGCGGGGGGCGCACACCGCGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTGGcccca                                                       >  2:760408/1‑149 (MQ=255)                                                                                                   aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGC‑GCGCGGGAAAATTTACACACGGTAGCCCCCCCACCAGAAAACGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGcaccac                                                      >  2:789132/1‑149 (MQ=255)                                                                                                    aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACCGTAGCCCCCCCATCAGAAAGCGCGGGGGGCACCCCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGcaccccct                                                    >  1:49679/1‑146 (MQ=255)                                                                                                       ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                <  1:313539/151‑1 (MQ=255)                                                                                                        gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                               <  2:231037/151‑1 (MQ=255)                                                                                                            gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAAGGTGGGGGGGACACCCAGCTGCGTAGCGGTCATCCCCTTCTGATAAAACAGCACCGGCACCACTTCGCggg                                              >  2:189500/1‑148 (MQ=255)                                                                                                              ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                         <  2:383633/151‑1 (MQ=255)                                                                                                               gaaagaaaAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGTGGGCCCACCGCGCTGCGTAGGGGTCATCCCTTTCTGATAAAACAACACTGGCACCACTTCGCggg                                              >  1:552476/1‑145 (MQ=255)                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCCCCACCTCCCGGGGAtggggg                                      >  1:859404/1‑149 (MQ=255)                                                                                                                          gCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  1:69351/151‑1 (MQ=255)                                                                                                                          gCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  1:352223/151‑1 (MQ=255)                                                                                                                            gACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                           <  1:771746/151‑1 (MQ=255)                                                                                                                             aCCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                          <  1:253852/151‑1 (MQ=255)                                                                                                                                 ttttGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGGAGACCCCCCCCCAGAAAACGGGGGGGGCACCACGAGCTGCGTCGCGATCATCCCCTGCTGATAAAACCGCACCGGCCCCACTTCGCCGGGATTGTGGCTGTGGCCAATgg                         >  2:478078/1‑148 (MQ=255)                                                                                                                                  tttGTGGCGGTTGG‑GCGCGGGGAAAATCACACAAGGTAGCCCCCACATCAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACTTCGCGGGGATGGGGGCGGTGGCCCATTGGCgg                     >  2:879936/1‑151 (MQ=255)                                                                                                                                      tggcggtggc‑gCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 >  1:662030/1‑151 (MQ=255)                                                                                                                                       ggcggtggc‑gCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                >  2:292527/1‑151 (MQ=255)                                                                                                                                             ggc‑gCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt          <  2:866858/151‑1 (MQ=255)                                                                                                                                                 gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt       <  1:442572/151‑1 (MQ=255)                                                                                                                                                 gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt       >  2:50657/1‑151 (MQ=255)                                                                                                                                                    cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg    >  1:387234/1‑151 (MQ=255)                                                                                                                                                       aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  2:218004/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508508‑3508806

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