Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,214,846	0	C	G	55.1%	33.7 / 106.0	69	H600Q (CAC→CAG)	rpoD	RNA polymerase sigma factor RpoD
Reads supporting (aligned to +/- strand):  ref base C (17/14);  new base G (18/20);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCC  >  NC_000913/3214576‑3215131                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                               gATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAg                                                                                                                                                                                                                                                                                >  1:247405/1‑287 (MQ=255)        gatgatgaAGGTTCGCATCTGGGGGATTTCATCGAGGATACCCCCCTCGTGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCg                                                                                                                                                                                                                                                                         <  1:308829/287‑1 (MQ=255)               aaGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCTTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTTCGTAGCt                                                                                                                                                                                                                                                                    >  1:346010/1‑285 (MQ=255)                               ggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCgg                                                                                                                                                                                                                                                  <  1:444422/287‑1 (MQ=255)                                  tttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGCTATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATTGGTAg                                                                                                                                                                                                                                               >  1:722433/1‑287 (MQ=255)                                  tttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGtgctgc                                                                                                                                                                                                                                                                          >  1:806540/1‑260 (MQ=255)                                                     cTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCg                                                                                                                                                                                                                                  <  1:223130/281‑1 (MQ=255)                                                      tCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGGAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTCGc                                                                                                                                                                                                                                                                     >  1:121922/1‑245 (MQ=255)                                                       cGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgccg                                                                                                                                                                                                                          >  1:97260/1‑287 (MQ=255)                                                                  tGGATTCTGCGACCACCGGAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAg                                                                                                                                                                                                                                                                                           <  1:90255/211‑1 (MQ=255)                                                                     aTTCTGCGACCACCGAAAGCCGGCGAGCGGCAACGCACGACGTGCTGGCTTGCCTGACCGCGCGGGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAGCGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAg                                                                                                                                                                                                                                                                                           <  1:389079/208‑1 (MQ=255)                                                                                            cgtgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCCCTAg                                                                                                                                                                                                             >  1:814143/1‑263 (MQ=255)                                                                                                   caccgcagggcaggcttgctgGCCTGCCCGGGCGGTAAGCAAACGTTCTGCGTATGGTTTTCGGTATCGAGATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGc                                                                                                                                                                                            <  2:97260/260‑1 (MQ=255)                                                                                                            acgTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTAt                                                                                                                                                                   <  1:596990/288‑1 (MQ=255)                                                                                                                gctggctgGCCTGACCGCGCGTTAAGCTAGAGTTCTGCATATGCGTTTCTGTATCGATATGAGCACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTCCCCGGGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAAGTGCGTCGCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTg                                                                                                                                                                                                                                                               <  1:757232/193‑1 (MQ=255)                                                                                                                gctggctgGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGGTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAAGGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAg                                                                                                                                                                                                                                                                      <  1:786744/186‑1 (MQ=255)                                                                                                                 ctggctggCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGcc                                                                                                                                                               >  1:461833/1‑287 (MQ=255)                                                                                                                 ctggctggCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGcc                                                                                                                                                               >  1:752739/1‑287 (MQ=255)                                                                                                                   ggctggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAg                                                                                                                                                                                                             >  1:92236/1‑240 (MQ=255)                                                                                                                   ggctggcCTGAACGCGCGGGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCa                                                                                                                                                             <  1:855781/287‑1 (MQ=255)                                                                                                                    gctggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACg                                                                                                                                                           <  1:153251/288‑1 (MQ=255)                                                                                                                     ctggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTCGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACCCCGCACCGGGCACTAGGCCCCCCGCACaactctctcacttg                                                                                                                                                                                    >  2:153251/1‑251 (MQ=255)                                                                                                                       ggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAg                                                                                                                                                                                                                                                                                           >  2:593691/1‑158 (MQ=255)                                                                                                                       ggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAg                                                                                                                                                                                                                                                                                           <  1:593691/158‑1 (MQ=255)                                                                                                                         ccTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTATGCCGGATCAGGCGTTACGCCGCCCCCGGCACTAGGCCCCCTGCCAAAACTGACCCTTTTCGGGGGGGGTTTTTTTTCGcgcctcccc                                                                                                                                                        >  2:342327/1‑279 (MQ=255)                                                                                                                         ccTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGc                                                                                                                                                                                                                  <  2:424695/229‑1 (MQ=255)                                                                                                                         ccTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGc                                                                                                                                                                                                                  >  1:424695/1‑229 (MQ=255)                                                                                                                                   cgTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGtg                                                                                                                                                                                                                                                                              <  1:517492/159‑1 (MQ=255)                                                                                                                                   cgTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGtg                                                                                                                                                                                                                                                                              >  2:517492/1‑159 (MQ=255)                                                                                                                                   cgTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATc                                                                                                                                                                                                                                       <  2:130655/198‑1 (MQ=255)                                                                                                                                   cgTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATc                                                                                                                                                                                                                                       >  1:130655/1‑198 (MQ=255)                                                                                                                                            aaaGTGCTGGGTATGCGTTTCGGTATCGATATGATCACCGACTACACGCTGGAAGAAGTGGGTAAACAGTGCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGGGCCGTTCTGAGGTGCTGCGTAGCTTCCTGGAc                                                                                                                                                                                                                                                            <  2:454069/168‑1 (MQ=255)                                                                                                                                                  cTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAACCGCCACCTTTTC‑GGTGGCGTTTTTTATCTCCCACGCACTACCAACCCCTGGTCCAGCACGCGATa                                                                                                                             >  1:735885/1‑288 (MQ=255)                                                                                                                                                     cGTATGCGTTTCGGTATCGATATGAACACTGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTACCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCg                                                                                                                                                                                                                                                   >  1:454069/1‑168 (MQ=255)                                                                                                                                                     cGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                                                                                                                                                      >  1:109012/1‑197 (MQ=255)                                                                                                                                                     cGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                                                                                                                                                      <  2:109012/197‑1 (MQ=255)                                                                                                                                                               tCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAg                                                                                                                                                                                                                                                                                           <  1:595777/118‑1 (MQ=255)                                                                                                                                                               tCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAg                                                                                                                                                                                                                                                                                           >  2:595777/1‑118 (MQ=255)                                                                                                                                                                cGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGt                                                                                                                >  1:191617/1‑287 (MQ=255)                                                                                                                                                                         tatGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCa                                                                                                                                                   <  1:416758/243‑1 (MQ=255)                                                                                                                                                                             aaCACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTg                                                                                                                                                                                                                                                                                   <  2:590461/112‑1 (MQ=255)                                                                                                                                                                             aaCACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTg                                                                                                                                                                                                                                                                                   >  1:590461/1‑112 (MQ=255)                                                                                                                                                                             aaCACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCg                                                                                                                                                                                                                                                                                        >  1:11411/1‑107 (MQ=255)                                                                                                                                                                             aaCACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCg                                                                                                                                                                                                                                                                                        <  2:11411/107‑1 (MQ=255)                                                                                                                                                                                  cGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAg                                                                                                                                                                                                                                                                                           >  2:553482/1‑99 (MQ=255)                                                                                                                                                                                  cGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAg                                                                                                                                                                                                                                                                                           <  1:553482/99‑1 (MQ=255)                                                                                                                                                                                  cGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCg                                                                                                                                                                                                                                  >  1:280444/1‑156 (MQ=255)                                                                                                                                                                                  cGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCg                                                                                                                                                                                                                                  <  2:280444/156‑1 (MQ=255)                                                                                                                                                                                   gACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATCCGCTTCAACCCGTTTCTCCAGTGAAACGCGACt                                                                                            >  1:507623/1‑288 (MQ=255)                                                                                                                                                                                    aCTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGAt                                                                                                                              >  1:599488/1‑253 (MQ=255)                                                                                                                                                                                    aCTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTGACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGAt                                                                                                                              <  2:599488/253‑1 (MQ=255)                                                                                                                                                                                         acGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGc                                                                                                                                                                                                           >  1:797836/1‑172 (MQ=255)                                                                                                                                                                                                       gggTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGcc                                                                                                                                                               >  2:414171/1‑201 (MQ=255)                                                                                                                                                                                                       gggTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGcc                                                                                                                                                               <  1:414171/201‑1 (MQ=255)                                                                                                                                                                                                                  tCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCaa                                                                >  1:299245/1‑285 (MQ=255)                                                                                                                                                                                                                  tCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGc                                                                                                                                                                                                                                             <  1:625617/113‑1 (MQ=255)                                                                                                                                                                                                                  tCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGc                                                                                                                                                                                                                                             >  2:625617/1‑113 (MQ=255)                                                                                                                                                                                                                                gAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCa                                                                                                                                                   >  2:823036/1‑188 (MQ=255)                                                                                                                                                                                                                                gAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCa                                                                                                                                                   <  1:823036/188‑1 (MQ=255)                                                                                                                                                                                                                                      aTCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGTACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGgca                                                                        <  1:545355/257‑1 (MQ=255)                                                                                                                                                                                                                                      aTCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTGTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAcg                                                                                                 <  1:556565/232‑1 (MQ=255)                                                                                                                                                                                                                                       tCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAg                                         <  1:632611/287‑1 (MQ=255)                                                                                                                                                                                                                                                  gaagctaagGCTCAGAGCAAACCGATTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTACTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACGAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCTACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTccc                             <  1:641088/288‑1 (MQ=255)                                                                                                                                                                                                                                                                           gTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                                                                                                                                                      <  1:467991/79‑1 (MQ=255)                                                                                                                                                                                                                                                                           gTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                                                                                                                                                      >  2:467991/1‑79 (MQ=255)                                                                                                                                                                                                                                                                           gTCAGCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAAt                                               >  1:513519/1‑245 (MQ=255)                                                                                                                                                                                                                                                                           gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                                                                                                                                                      <  2:503485/79‑1 (MQ=255)                                                                                                                                                                                                                                                                           gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCAcc                                                                                                                                                                                                                      >  1:503485/1‑79 (MQ=255)                                                                                                                                                                                                                                                                              aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTGGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATcc  <  1:342327/287‑1 (MQ=255)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                               GATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTC‑GGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCC  >  NC_000913/3214576‑3215131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑