Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,976,539	G→A	L201F (CTT→TTT)	motA ←	motility protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,976,539	0	G	A	96.3%	171.2 / ‑4.4	55	L201F (CTT→TTT)	motA	motility protein A
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (34/18);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCA  >  NC_000913/1976292‑1976771                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         aCCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGACTAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGGGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAacgc                                                                                                                                                                                                                                       >  1:355893/1‑251 (MQ=255)       cgcacATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAg                                                                                                                                                                                                                                 <  1:526204/251‑1 (MQ=255)               ctctTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTcggcgg                                                                                                                                                                                                                         >  1:374546/1‑251 (MQ=255)               ctctTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTcggcgg                                                                                                                                                                                                                         >  1:516139/1‑251 (MQ=255)                         tCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCAGCGCATGTGCGATAACCGCACCCAGCTCGGCGGCAGGACgctc                                                                                                                                                                                                               >  1:2063/1‑248 (MQ=255)                                  gacggacgTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAAcc                                                                                                                                                                                                      >  1:281336/1‑251 (MQ=255)                                    cggccgttcgaaggAATAGAGCGGTTTGCGACCAAACTCACCGGCGCTAGGCGGTGAGAAACCGGTCAGATTAGAAAGCAGAGTGACTTTGACGCAATGCATAATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATCAACGCACCCAGCTCGGCGGCAGGACGATCGGCAGAACCt                                                                                                                                                                                                     >  2:327518/13‑250 (MQ=255)                                    cggacgTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCt                                                                                                                                                                                                     >  2:57037/1‑250 (MQ=255)                                    cggacgTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTa                                                                                                                                                                                                    <  2:102470/251‑1 (MQ=255)                                                 gaagagagCGTTTTGCGACCAATCTCAACGGCGAGAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTGGACTCACTGCATCATTTTGCTGGTTTCGGCGCTTTGCTGACGTAAAAC‑AGTCGCGAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTGCCCGCCAGTGCATGTGCGATAAACTCACCCAGCTCGGCGGCAGGACGATTGGTTGAACCTAACGTGTGAACGAc                                                                                                                                                                                       <  2:355893/247‑1 (MQ=255)                                                     agagCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGAccccc                                                                                                                                                                                   >  1:454803/1‑251 (MQ=255)                                                         gggtgtGCGACCAAACGAAAAGGTGATAGGAGGTGCGTAACCGTTAAGATTAGAAAGAAGAGTGAGTGTGACGAACTGCATCATTTTGCTGGTGTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGACGAGGAAAGTCTCCACCAGCGCATGTGCGATAAGCGCACCCTGCGAGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAATGAccgcaagga                                                                                                                                                                               <  2:232519/246‑8 (MQ=255)                                                         cGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTa                                                                                                                                                                               <  1:153410/251‑1 (MQ=255)                                                         cGTTTTGCGACCAAACGCAACGGCGAGAGGCGGTGCGTAACCGTTCAGAGTAGAAAGCAGAGTGACTTTGACGCACGGCTTCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCTCCACCATCGCATGTGCGATAAGGGCACCCAGCTCGGCGGCAGGGCGATCGGCTGAACCTAACGCGTGAACggcgccgaggg                                                                                                                                                                               <  2:232537/251‑11 (MQ=255)                                                                      aaCTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTGGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCCCCATCGCATGTGCGATAAACGCCCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGTCCCCCATTACAGCCGCAACAAt                                                                                                                                                                  <  1:136928/251‑1 (MQ=255)                                                                            aCGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaaa                                                                                                                                                            <  2:152072/251‑1 (MQ=255)                                                                             cGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAAc                                                                                                                                                           >  1:78568/1‑251 (MQ=255)                                                                                    aGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGCGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGa                                                                                                                                                     >  1:481906/1‑250 (MQ=255)                                                                                     ggCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAAAACCAACCGCCGGa                                                                                                                                                     >  2:184042/1‑249 (MQ=255)                                                                                                  cGTTCAGATTAGAAAGCAGAGTGACTTTGACGCCCTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGAcc                                                                                                                                      <  2:516139/251‑1 (MQ=255)                                                                                                       aGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcgc                                                                                                                                 >  1:378690/1‑251 (MQ=255)                                                                                                                aGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGtt                                                                                                                        >  1:533963/1‑251 (MQ=255)                                                                                                                      gTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCgg                                                                                                                  <  1:57037/251‑1 (MQ=255)                                                                                                                         aCTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGAc                                                                                                               <  1:275928/251‑1 (MQ=255)                                                                                                                          cTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACa                                                                                                                                                                             >  1:293768/1‑188 (MQ=255)                                                                                                                          cTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACa                                                                                                                                                                             <  2:293768/188‑1 (MQ=255)                                                                                                                              gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTg                                                                                                          <  2:162270/251‑1 (MQ=255)                                                                                                                                  cACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTACAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACAACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGGAAGTCCCCGACACGCGCCAGACTGTTCGCCGGGACTTCTGCCTc                                                                                                      >  2:465992/1‑251 (MQ=255)                                                                                                                                   aCTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCCAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCg                                                                                                     >  2:28945/1‑251 (MQ=255)                                                                                                                                             tttGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCg                                                                                           >  1:372179/1‑251 (MQ=255)                                                                                                                                                gCTGGGTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAAGGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACTCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACGAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtct                                                                                        <  1:233722/251‑1 (MQ=255)                                                                                                                                                  tGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCAAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCACGACCAGAGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtctc                                                                                       >  1:213090/1‑250 (MQ=255)                                                                                                                                                              ttttCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAACTCCGTAAGCCAATAAACTGCCGAGGAACGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGc                                                                                                    >  2:256870/1‑225 (MQ=255)                                                                                                                                                              ttttCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGc                                                                                                    <  1:256870/225‑1 (MQ=255)                                                                                                                                                                   tGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCACATAAAATGACGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAAAGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGACGAGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTACCTTTCGTGCGTCTCAATCTCTTcatccacca                                                                      >  2:335049/1‑247 (MQ=255)                                                                                                                                                                     gactAAAACCAGTAGCTACTGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGAGCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCAGCcagcag                                                                     <  2:132419/247‑1 (MQ=255)                                                                                                                                                                           aaC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCACACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTCTCGTGAGTCTCACTCTCTTCATCCATCAGAGCTTCgc                                                             >  2:279624/1‑250 (MQ=255)                                                                                                                                                                             c‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCGCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGAtt                                                           <  1:279624/251‑1 (MQ=255)                                                                                                                                                                               aGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAGTGCCGAGGAAAGTCCCCACCATCGCATGTGCGAAAAACGCACCCAGCTCGGCGGCAGGCCGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGGGTCCCCGACCCGAGCAAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGGGCGTCACAATCTCTTCATCCATCAGAGCTTCGCttt                                                          >  1:486281/1‑251 (MQ=255)                                                                                                                                                                                      aaTGGGGAAATAAATCCGTAAGCCAATACAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGAGAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCAGCAGAGCTTCGATTTCGAAGgt                                                   <  2:2063/251‑1 (MQ=255)                                                                                                                                                                                       aTGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGgtg                                                  >  1:333304/1‑251 (MQ=255)                                                                                                                                                                                         ggggAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGCACCTAACGCGTGAACGACCCCCATTACAGCAGCAACAATACCAAACGCCGGACGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCTATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGtt                                                >  1:252027/1‑251 (MQ=255)                                                                                                                                                                                                 aaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTTTTCATGTGAc                                        >  2:421449/1‑251 (MQ=255)                                                                                                                                                                                                 aaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGAc                                        >  1:66830/1‑251 (MQ=255)                                                                                                                                                                                                 aaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGAGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGAc                                        >  1:244368/1‑251 (MQ=255)                                                                                                                                                                                                     ccGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCt                                    >  1:108240/1‑251 (MQ=255)                                                                                                                                                                                                         aaGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGAGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGAt                                 >  1:463393/1‑250 (MQ=255)                                                                                                                                                                                                           gCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAAt                              >  1:51639/1‑251 (MQ=255)                                                                                                                                                                                                             cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTcccc                                                                                                                                          >  2:108610/1‑141 (MQ=255)                                                                                                                                                                                                             cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTcccc                                                                                                                                          <  1:108610/141‑1 (MQ=255)                                                                                                                                                                                                              aaTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTcatcca                                                                         >  1:60446/1‑205 (MQ=255)                                                                                                                                                                                                                        cGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATa                  >  2:454045/1‑250 (MQ=255)                                                                                                                                                                                                                           ggAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATcg              >  1:276614/1‑251 (MQ=255)                                                                                                                                                                                                                                gTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGAt         >  1:138235/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccaTCGCATGTGCGATAAACGCTCCCAGCCCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGCCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGCCGATAAAATCa  <  2:369149/251‑1 (MQ=255)                                                                                                                                                                                                                                        caTCGCATGAGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGAGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCa  >  2:429378/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         ACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAACGGCGATAGGCGGTGCGTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAAC‑AGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCA  >  NC_000913/1976292‑1976771

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