Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,332,310	C→T	G321G (GGG→GGA)	basS ←	sensor histidine kinase BasS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,310	0	C	T	93.2%	134.9 / ‑3.7	45	G321G (GGG→GGA)	basS	sensor histidine kinase BasS
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (38/3);  total (38/7)
Fisher's exact test for biased strand distribution p-value = 1.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑CCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCG  >  NC_000913/4332066‑4332547                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              cTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑tccc                                                                                                                                                                                                                                            >  2:456147/1‑248 (MQ=255)    aaTCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGAcc                                                                                                                                                                                                                                       >  1:34450/1‑250 (MQ=255)     aTCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCa                                                                                                                                                                                                                                      >  1:646130/1‑250 (MQ=255)         ctGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAAt                                                                                                                                                                                                                                    <  1:5920/248‑1 (MQ=255)           gACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTCGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATcc                                                                                                                                                                                                                                  <  2:676733/248‑1 (MQ=255)           gACGCGCATACTCTCCGCCAGGTGAACGGAGGAGAGTGAAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGGTCAGCGGGCTGGTGGTCAGCAGCTTTCTTGATATCTGGGTTGCCACGTACTGATCTTTCTTCAGCCGGACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTGAA‑TCCCAGACCAATCCCGc                                                                                                                                                                                                                               <  2:90985/251‑1 (MQ=255)             cgcgCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGcc                                                                                                                                                                                                                              >  1:33068/1‑250 (MQ=255)                gcATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGTGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAGAACTGCCCGTGACGCAACT‑GTGTAATGCGGCTGACACTACTTAA‑TCCCAGACCAATCCCGCCATa                                                                                                                                                                                                                           >  1:382236/1‑250 (MQ=255)                gcATACTCTCCACCAGGTTAACGGAGGTGAGTGCAATGAAACACCGTGTTTATGAAAGTTTCACTTCCGTGTTCAGCGTGCCGGTGGTCAGCACCTTTCTTTATATCTGGTTTCCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGCAGTCTCTTGCCGGTTTTGCCGGAACAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCAGCCATa                                                                                                                                                                                                                           >  2:738160/1‑250 (MQ=255)                 cATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAAc                                                                                                                                                                                                                         >  1:399277/1‑251 (MQ=255)                             cAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATAc                                                                                                                                                                                                             >  1:482394/1‑251 (MQ=255)                                  tAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTAc                                                                                                                                                                                                         >  1:628815/1‑250 (MQ=255)                                   aaCGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACa                                                                                                                                                                                                        >  1:225958/1‑250 (MQ=255)                                    aCGGAGGAGAGTGCAATGAAAACCCGTGTTTATGACAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACaa                                                                                                                                                                                                       >  1:202634/1‑250 (MQ=255)                                                                tttATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTActc                                                                                                                                                                          >  2:294534/1‑251 (MQ=255)                                                                       aaGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAACAACTGCCCGTGATGCAACT‑GTGTCATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCCCATTTACTCtcctca                                                                                                                                                                    >  2:497391/1‑250 (MQ=255)                                                                       aaGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTAc                                                                                                                                                                            >  2:100273/1‑242 (MQ=255)                                                                          tttAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATAc                                                                                                                                                                >  2:662331/1‑251 (MQ=255)                                                                          tttAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATAc                                                                                                                                                                >  2:480312/1‑251 (MQ=255)                                                                          tttAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGCGCCGGAAGCCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATAc                                                                                                                                                                >  1:379000/1‑251 (MQ=255)                                                                                            cGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCa                                                                                                                                               >  1:164247/1‑250 (MQ=255)                                                                                             gTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAAc                                                                                                                                             >  1:466090/1‑251 (MQ=255)                                                                                                        cagcagCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCCGACGTACCAAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGa                                                                                                                                                                                                                                         >  1:411475/1‑148 (MQ=255)                                                                                                        cagcagCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCTGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GGGTAATGCGGCTGACAATACTTAA‑GCCCAGa                                                                                                                                                                                                                                         <  2:411475/148‑1 (MQ=255)                                                                                                                  ctttATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGACCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCAtc                                                                                                                        >  2:254670/1‑251 (MQ=255)                                                                                                                      atatCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCAtcttct                                                                                                                    >  2:284245/1‑251 (MQ=255)                                                                                                                          cggttgtcccACGTACTTCTCTTCCTTCACCCACACCCGGGCCCGGGTGCGGGAAGTCTCTGGCGGTTTTTGAAAAAAAACCCCCCCGTAATCCAATT‑GTGAAAGGCGGCTGAAAATACTAAA‑CCCCAGCCCAACCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCa                                                                                                                <  2:276402/243‑1 (MQ=255)                                                                                                                          cTGGTATGCCACGTACTGATCTTTCTTCAGCCGTGCCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAACAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTAATCAATACCTGGTCCTTCATCTTAAACTGCCATGACCGCTCCGTCATCTTCTTGcc                                                                                                                >  2:709162/1‑250 (MQ=255)                                                                                                                                  ccACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCCGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATc                                                                                                        >  2:709196/1‑251 (MQ=255)                                                                                                                                   cACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATc                                                                                                        >  1:685770/1‑250 (MQ=255)                                                                                                                                                  attcAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGCCCAATCCCGCCATAACGGCTGTCCATACGAAAAAACGCTTTACTCAACGCCCCGCATTTACTCTCATCAATACCTGGTCATTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATAATAATGTTGCTGCCtg                                                                                        >  2:214087/2‑250 (MQ=255)                                                                                                                                                                    cggcgTGCCGAAAGTCTCTTGCCGTTTGTGAAGAAAACACGCCCCGGGA‑GCAATTTGTGTAAGGCGGCTAAAAATACTAAA‑CCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATggg                                                                      <  1:456147/249‑1 (MQ=255)                                                                                                                                                                               aaGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCACGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTAcc                                                           >  2:452894/1‑251 (MQ=255)                                                                                                                                                                                   ctctTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGt                                                       >  1:681099/1‑251 (MQ=255)                                                                                                                                                                                                 gCAGGAAAACCTGCCCGTGATGCAACT‑GTGTAATGCGACTGACAATACTTAC‑TCCCACACCAATCCCACCATAACGGCTGTCCCTACGTACAACCGCTTTACTCCACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCACATACCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGCTGGGCGTTTTCTTCCAGGGTTCGCACTAACATc                                         >  2:736821/1‑251 (MQ=255)                                                                                                                                                                                                     gAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGACCTTCATCTTCAACTGCCATGACCGCCCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCgc                                      >  2:493368/1‑250 (MQ=255)                                                                                                                                                                                                            tGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGt                              >  1:170219/1‑251 (MQ=255)                                                                                                                                                                                                               ccGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCCATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGc                           >  2:278431/1‑251 (MQ=255)                                                                                                                                                                                                               ccGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCCTCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGc                           >  1:514681/1‑251 (MQ=255)                                                                                                                                                                                                                cGTGATGCAACT‑GTGTAATGCGGCTGCCAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGc                           >  2:644828/1‑250 (MQ=255)                                                                                                                                                                                                                           t‑gtgTAATGAGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGAAATAACATCCGCAGCAGGGTTGCATCACCCTGacc               >  2:619453/1‑249 (MQ=255)                                                                                                                                                                                                                                 aaTGCGGCTGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTg           >  2:581729/1‑250 (MQ=255)                                                                                                                                                                                                                                       gCTGACAAT‑CTTAATCCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcgg    <  1:254670/251‑1 (MQ=255)                                                                                                                                                                                                                                         ttaCAATCCTTAA‑TCCCAGACCAATCCCGCCATACCGGCTGTCAATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTATCATCACTACCTGAACCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGGTTTCTACCAGGTTTCGCAATAACATCCGCCGCCGGGTTGCATCACCCTGAACGGTGATGTcggcg  >  2:653807/3‑251 (MQ=255)                                                                                                                                                                                                                                         tGACACTACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcggcg  >  2:304649/1‑251 (MQ=255)                                                                                                                                                                                                                                         tGACAATACTTAA‑TCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGCTGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGCTGTcggcg  >  2:411824/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              CTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACT‑GTGTAATGCGGCTGACAATACTTAA‑CCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCG  >  NC_000913/4332066‑4332547

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