Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,516,959	C→T	P253P (CCG→CCA)	fecR ←	ferric citrate regulator FecR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,516,959	0	C	T	93.0%	105.0 / ‑1.7	46	P253P (CCG→CCA)	fecR	ferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (18/22);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCC  >  NC_000913/4516724‑4517195                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             agaAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGa                                                                                                                                                                                                                                 >  1:218393/1‑249 (MQ=255)         aaaaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGTTCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCAGCCGGTGCGGGAACTGGTTAGCGGTGCTTA‑CACCTCCCCCAGTGGGTTATCGCTGAAGCTCAggct                                                                                                                                                                                                                       >  1:341145/1‑249 (MQ=255)         aaaaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGAATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGGTTTTCGCTGTAGCTCAgggt                                                                                                                                                                                                                       >  1:402186/1‑249 (MQ=255)         aaaaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTAAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCTCAGCGCCGCACGCCCGTGCGGTAACGGGTGAGGGTGGGTCT‑CACCTCACCCAGTGGATTATAGCTGTAGCTCAgg                                                                                                                                                                                                                         >  2:43867/1‑249 (MQ=255)             aGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTcc                                                                                                                                                                                                                   <  2:56856/251‑1 (MQ=255)                acaacctaaaTAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTccgtc                                                                                                                                                                                                                <  2:175706/244‑4 (MQ=255)                aCAATCGAAATAAGAATGATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTc                                                                                                                                                                                                                <  2:106368/251‑1 (MQ=255)                           aaGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAACGCAGCGCAGCACGCCGTTTCGGGAAAGGGGTAGCGTGGC‑ATACACCTCCCCCCGAGGGTTATCGCTGAAGACCAGGGTGGCCTTCGTCCAGCTTGt                                                                                                                                                                                                     >  2:464096/1‑251 (MQ=255)                           aaGAAATATTTGCCTTAAAGTGGTGAAATGTTTATCCATGACCGCGTAATAGACTGAAGGTTAACGGGAAGCGTGTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAGCGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGTGTTAGCGTGGCTAT‑CACCTCCCCCAGTGGTTTATCGCTGAAGCTCAGGATGTACGTCGGCCCGCTTGt                                                                                                                                                                                                     <  2:221497/251‑1 (MQ=255)                               attattTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCCGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCCTCGTCCAGCTTGTActc                                                                                                                                                                                                 >  1:45396/1‑251 (MQ=255)                                    tttCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGAt                                                                                                                                                                                                                       <  2:51595/224‑1 (MQ=255)                                    tttCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGGTAGCGTGGCTAT‑CACCTCACCCAGGGGTTTATCGCTGAAGCTCAgggt                                                                                                                                                                                                                       >  1:51595/1‑222 (MQ=255)                                       ccTTCCAGTGTTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGTGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCATCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTGT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCGCGTCagc                                                                                                                                                                                         <  2:223075/251‑3 (MQ=255)                                          tACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCGCGTCATcc                                                                                                                                                                                        <  1:474673/249‑1 (MQ=255)                                           aCAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGGGGCTAT‑CACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGc                                                                                                                                                                                     >  1:171639/1‑251 (MQ=255)                                                    aaaTGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACt                                                                                                                                                                            <  1:52542/251‑1 (MQ=255)                                                                 gTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATTGCAGCGCCGCCCGCCCTTGCGGTAACGGGTTAGCCGGGCTAT‑CACCCCACCCAGCGGTTTATCGCTGAAGCTCACGATGTCCTTCGCCCAGCTTGTACTCTCGTCATCCCGCGGtttttatccgccaaactccg                                                                                                                                                               >  2:24332/1‑234 (MQ=255)                                                                        gTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCt                                                                                                                                                        <  1:376740/251‑1 (MQ=255)                                                                           aTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGGAACGGGGTAGCGTGGGTAT‑CACCCCACCCAGTGGGTTTTCGCTGAAGCTCCGGGTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCCGTTTCACTGCGCCCAACTCACGGGCGCTGa                                                                                                                                                      >  2:96847/1‑250 (MQ=255)                                                                                      tttAACGGGAAGCGTTTGCGCGATAACGGTCAGGATCGCATCGGTATTTTTACGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTAGCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCtt                                                                                                                                           <  2:2942/250‑1 (MQ=255)                                                                                          aCGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGGTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGGACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                                                                                                      >  1:112237/1‑251 (MQ=255)                                                                                               aaGCTTTTGCGCGACAACGTTCAGGATCGCATCGGTAATGTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGt                                                                                                                                 <  1:157870/251‑1 (MQ=255)                                                                                                   ggttGCGCGATAACGTTCAGGATCGCAACGGTAGTTTTCAGCGGGAACGGCCCGCTCAGGCGGAGCCCGGCAACGGGGGGATCGCAGCGCATCGAGTCGTTGCGGTAACGGGTTAGCGTGGCTAG‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTTCTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGTACTGTAGGCTTTCACCAGCGTTCAc                                                                                                                             <  2:341203/249‑1 (MQ=255)                                                                                                                 gTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTgcg                                                                                                               >  1:9466/1‑251 (MQ=255)                                                                                                                       gATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggggg                                                                                                         <  1:153022/251‑1 (MQ=255)                                                                                                                          cGCAGCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCAc                                                                                                      <  2:366458/251‑1 (MQ=255)                                                                                                                           gCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCCGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTTTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCCTGTCCTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGACGAACTGCCGGCTTTCACCAGCGGTCCCGATGCGTTTTTTCGCGGGGGCAcc                                                                                                     >  2:434129/1‑250 (MQ=255)                                                                                                                            cATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGGTCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTAT‑CACCTCACCCAGGGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCCAACTCAGAAGCGCTGAACTTCAGACTTTTACCAGCGGTCACGATGGGGTTTTGCGCGGGGGCACTg                                                                                                    >  2:265867/1‑251 (MQ=255)                                                                                                                             aTCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCATCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCACCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTgg                                                                                                   <  1:239895/251‑1 (MQ=255)                                                                                                                                              ggAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGGTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCCCCCAGTGGTTTTTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCCTCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCGGGCTTTCACCCGCGTTCACGCTGCGTTTTTTCTCGGGGGCACTGGCGAGAAGCCCTTcccc                                                                                   >  1:84791/1‑248 (MQ=255)                                                                                                                                              ggAACGGCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCca                                                                                    <  1:504886/249‑1 (MQ=255)                                                                                                                                               gAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAAATGCAGGCTTTCACCAGCGGTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCcac                                                                                   >  1:7596/1‑249 (MQ=255)                                                                                                                                                    tCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGGGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTACTTCGTCCAGCCTGTACTCTCGTCATACAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGAGTTTTGGCGCGGGGGCACTGGCGAGACGCACTTCCACAGCGtgc                                                                            >  1:380066/1‑251 (MQ=255)                                                                                                                                                      ccGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctg                                                                          >  1:136357/1‑251 (MQ=255)                                                                                                                                                           cAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctgctg                                                                       <  2:112237/249‑1 (MQ=255)                                                                                                                                                                  aGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAg                                                              <  2:365743/251‑1 (MQ=255)                                                                                                                                                                       ggCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCg                                                         >  1:46663/1‑251 (MQ=255)                                                                                                                                                                                gggATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTAt                                                <  1:96847/251‑1 (MQ=255)                                                                                                                                                                                     cgcagcgcagCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCATGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCctg                                            <  2:218393/250‑1 (MQ=255)                                                                                                                                                                                       cagcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGCAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCCCCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGCAATTATCctgctg                                         >  2:137885/1‑251 (MQ=255)                                                                                                                                                                                                  gCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTaa                              <  2:378132/251‑1 (MQ=255)                                                                                                                                                                                                         cGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGt                       <  1:22981/251‑1 (MQ=255)                                                                                                                                                                                                                gggTTAGCGTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTaa                 >  1:4032/1‑250 (MQ=255)                                                                                                                                                                                                                        gTGGCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTg          >  1:47931/1‑249 (MQ=255)                                                                                                                                                                                                                          ggCTAT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCt      <  1:248580/251‑1 (MQ=255)                                                                                                                                                                                                                              aT‑CACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGcc  >  1:160393/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             AGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAT‑CACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCC  >  NC_000913/4516724‑4517195

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