Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,003,276	A→C	S111A (TCT→GCT)	fliC ←	flagellar filament structural protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,003,276	0	A	C	95.5%	234.7 / 0.9	66	S111A (TCT→GCT)	fliC	flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (39/24);  total (42/24)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

AGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGAC  >  NC_000913/2003031‑2003511                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGAc                                                                                                                                                                                                                                         <  1:292986/250‑1 (MQ=255)    aacaacaGCAGTTACGTGAACACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATAATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCATCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCACTTTCATCCAGCCGGGATTTATTTTCGTCCTGGATAGCAGACatat                                                                                                                                                                                                                                     >  2:516341/7‑248 (MQ=255)      cAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATc                                                                                                                                                                                                                                    >  2:46515/1‑250 (MQ=255)      cAAAAGCAGTTACTAGAGCACTAGAGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATc                                                                                                                                                                                                                                    >  1:777021/1‑250 (MQ=255)                                 aaacgtgatagTTAATTTTAACGCTAAACCCATCAAGGCCAAGAGGTTTAGCATCAATCTTCGTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCACCCGGGATTTGCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTgtcc                                                                                                                                                                                                        <  2:21906/241‑4 (MQ=255)                                      gTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAc                                                                                                                                                                                                   >  1:413601/1‑251 (MQ=255)                                        aTCGTTATTTTTAACGCTAAATCCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAcc                                                                                                                                                                                                  <  1:265498/250‑1 (MQ=255)                                        aTCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAcc                                                                                                                                                                                                  <  1:869230/250‑1 (MQ=255)                                        aTCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCg                                                                                                                                                                                                 >  1:66538/1‑251 (MQ=255)                                        aTCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCAGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCg                                                                                                                                                                                                 >  1:424227/1‑251 (MQ=255)                                           gttttttttAACGCTAAAAACATCAAGACCAATAGTTTTAGCATCAATCTTTTTTATATCGATAGTGATAGTCTGGTTTTCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATGTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCa                                                                                                                                                                                              <  2:538953/247‑1 (MQ=255)                                                    aaCGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGcc                                                                                                                                                                                     >  1:367974/1‑250 (MQ=255)                                                     aCGCTAAAACCTCCAAGGCCAAGAGTTTTAGCAACAATCTGCTTCAGATCGATAGGGATAGTCTGTTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCGTTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTgcacgcac                                                                                                                                                                                    <  2:370068/251‑1 (MQ=255)                                                      cGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGcacc                                                                                                                                                                                   >  1:393865/1‑250 (MQ=255)                                                          aaaaCCATCAAGGCCAAGAGTTTTAGCATAAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGAATCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGAACACGCt                                                                                                                                                                               >  1:516061/1‑251 (MQ=255)                                                                  cAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAgttg                                                                                                                                                                       >  1:114907/1‑251 (MQ=255)                                                                              ttttAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTACGTTGTTGTTGATTTCg                                                                                                                                                           >  1:6500/1‑251 (MQ=255)                                                                                ttAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCgg                                                                                                                                                          >  1:603922/1‑250 (MQ=255)                                                                                     aTCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACATGGATTTTCCTGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTAGGTCTGACCAGATACCCGGTCAATTTCATCCAGACGGGCTTTAATTTCGTCCTGGATAGCAGACAGAACAGACTCCGAGTTAGAACCGGTAGTGGCCTGTACCGTCAGTTCACGCAACCGCTGTAAGTTGTTGATGATTTCggcccgcg                                                                                                                                                    >  1:549569/1‑245 (MQ=255)                                                                                          tatgcatCAGATCGCTGGTGATAGACTGGTTATCATTTGCGCCAACCTGTATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGACGTTGAACTGGGTCTGCCCAGATCCGCGGACAATTTCATCCAGACGGGATTTAATTTCGACCTGGATAGCAGAACGCTCAGATTCAGAGTTAGTACCGGTAGTGGACTGTTACTTCAGTTCACGCACACGCCGGAATATGTTGTTGATTACGGACAACGAGCCt                                                                                                                                               >  2:513667/7‑251 (MQ=255)                                                                                            gcattCAGATCGAACGTGATGTTCGGGTTATCGTTTGTGCAAACCGGGATTTGCATTGAGCCTTTTTTGGCCATCCTTTTCCCGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAAGTGCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTATTGGCCTGGACCGTCAGTTCACGCATACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTc                                                                                                                                             <  2:433472/248‑1 (MQ=255)                                                                                                tcataTCGATTGTGATATGCTGGTTATCATTTGCGCCAACCTGGATTGTCATGGTGCCATTTGTTGCCAGCACGTTCACGCCGTTGAACTGGGGCTGACCAGATACGCGGTCAATTTCATCCAGACTGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTTggt                                                                                                                                          <  2:739182/246‑1 (MQ=255)                                                                                                tCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGCCCTGGATAGCAGACCGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTATACCGTCAGTTCACGCACCCGCTGTAAGTTGTTGTTGATTTCGGACAGAGCGCCTTCggtg                                                                                                                                         >  2:855121/1‑251 (MQ=255)                                                                                                 cAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCggtgg                                                                                                                                        >  1:824694/1‑251 (MQ=255)                                                                                                    aTCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTTTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCt                                                                                                                                     >  1:298350/1‑251 (MQ=255)                                                                                                      cGATAGTGAGAGTCTGGTTATCATTTGCGCCAACCGGGATTTTCATGGAGCCATTTTTTGCCAGCTCGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTgc                                                                                                                                   <  2:706288/251‑1 (MQ=255)                                                                                                           gtgatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGCGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAAc                                                                                                                              <  2:57217/251‑1 (MQ=255)                                                                                                            tgatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACg                                                                                                                             >  1:610620/1‑251 (MQ=255)                                                                                                             gatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACgg                                                                                                                            >  1:653922/1‑251 (MQ=255)                                                                                                             gatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACg                                                                                                                             <  1:5836/250‑1 (MQ=255)                                                                                                                agtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGGCTGCGCAACgg                                                                                                                            >  1:764829/1‑248 (MQ=255)                                                                                                                        tATCATTTACGAAAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTGAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTCCCGGTAGTGGCCTGTACAGTCAGTTCACTCACCCGCTGTACGTTGTTGTTGGTTTCGGCCAGCGCGCATTCGGTGGTATGCGCCCCGGAGATACcgttg                                                                                                                 >  2:72035/1‑249 (MQ=255)                                                                                                                           cATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGtt                                                                                                               <  1:654556/250‑1 (MQ=255)                                                                                                                           cATTTGCGCCAACCGGGATTTGCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTTAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTg                                                                                                              <  2:440394/251‑1 (MQ=255)                                                                                                                            aTTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTgg                                                                                                             >  1:124193/1‑251 (MQ=255)                                                                                                                              ttGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGCCGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACCCGCTGTAAGTTGGTGTTGGTTTCGGACAGCGCGCCTTCGGGGGTATGCGCAACGGCAATACCGTCGTTGGCg                                                                                                           >  2:65969/1‑251 (MQ=255)                                                                                                                                 cgcCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTCCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGtt                                                                                                         >  1:416277/1‑250 (MQ=255)                                                                                                                                   ccAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGCTAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCCGTTCACGCACACGATGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCCACGGAGATACCGTCTTTGGCGTTACg                                                                                                      >  1:428067/1‑251 (MQ=255)                                                                                                                                     aaCCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACggg                                                                                                    >  1:54791/1‑251 (MQ=255)                                                                                                                                     aaCCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGGTGGCGTTACggg                                                                                                    >  1:191472/1‑251 (MQ=255)                                                                                                                                       ccgggATTTTCATGGAGCAATTTTTTGCCGGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGgcc                                                                                                  <  1:149370/248‑1 (MQ=255)                                                                                                                                                    ggAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGc                                                                                     <  1:18820/251‑1 (MQ=255)                                                                                                                                                                      cACGTTCACGCCGTTGAACTGGGTCTGTCCAGATACGCGGTTAATTTCATCCATACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGa                                                                   <  2:774435/251‑1 (MQ=255)                                                                                                                                                                          ttCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGGGGTCTGCGCAACGGCGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGGGGTGAAACg                                                               >  1:650835/1‑251 (MQ=255)                                                                                                                                                                               gCCGTTGAACAGGGACTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAACTTCGTCCTGGATCGAAGACAGCTCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTACACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGCGGGCTGCGCAACGGAGGTACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCCT‑GATGTTCGAAGTGCAACGGTTAGc                                                         >  2:395720/1‑251 (MQ=255)                                                                                                                                                                                   ttGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAAt                                                      >  1:665897/1‑251 (MQ=255)                                                                                                                                                                                       aCTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCgac                                                  <  1:121194/251‑3 (MQ=255)                                                                                                                                                                                             tCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGGGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAAc                                                                >  2:192873/1‑231 (MQ=255)                                                                                                                                                                                             tCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAAc                                                                <  1:192873/231‑1 (MQ=255)                                                                                                                                                                                              cTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCg                                           >  1:198879/1‑251 (MQ=255)                                                                                                                                                                                                gACCAGATAGGCGGTCAATTTCATCCAGACGGGAGTGAATTTCGGCCTGGATAGCAGACAGATCAGACTCAGAGTGAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTGCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCt                                         <  2:653922/251‑1 (MQ=255)                                                                                                                                                                                                gACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTGAGAGGTGAAACGGTTAGCAATCGACTGACCCGCt                                         <  1:654170/251‑1 (MQ=255)                                                                                                                                                                                                     gATACGCGGTCAATTTCATCCAGACGGGATTTACTTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTATGCGCCACGGAGATCCCGTCGTTGGCTTTACGGGACGCCTGAGTCAGGCCTTTAATGTTAACGGTGAAACGGTTAGCCATCGACTGACCCGCTTCGTc                                    >  2:348408/1‑251 (MQ=255)                                                                                                                                                                                                           cgGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCtt                              <  2:22058/251‑1 (MQ=255)                                                                                                                                                                                                            ggTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTc                             >  1:301526/1‑251 (MQ=255)                                                                                                                                                                                                            ggTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCAAACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGGGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAAAGTTAGAGGTGAAGCGGTTAGCAATCGAATGACCCGCTGCGGAATCCTTc                             >  2:281138/1‑251 (MQ=255)                                                                                                                                                                                                                   ttCATCCAGACGGGCTTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGGGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCgcgcgtt                      >  2:345601/1‑248 (MQ=255)                                                                                                                                                                                                                        ccAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAg                                                          <  1:725995/210‑1 (MQ=255)                                                                                                                                                                                                                                gATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGAAAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAAGGTTAGAGGTGGAACGGTTATCAATCGACTGACCCGCTGAGTCATCCTTCGCGCTGTTAATACGCAAGcc         >  2:227167/1‑251 (MQ=255)                                                                                                                                                                                                                                 aTTTAATTTCGTCCTGGATAGCAGACAGATCAGTCTCAGTGTTAGTACCGGTAGTGGCCTGTACCGTTAGTTCACGCACCCGCTGTAAGTTGTTGTTGATTTCGGATAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCa        <  2:665852/251‑1 (MQ=255)                                                                                                                                                                                                                                    tAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGATGTTGTTGATTTCGGACAGCACGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGGGTTACGGGCCGCCTGAGTCAGGCCTTTTATGTTAGAGGAGATACGGTTAGCAATCGACTGACCCGCTGCGGCATCCTTAGCGATGTTAATACGCAACCCaga      >  1:105918/1‑250 (MQ=255)                                                                                                                                                                                                                                       tttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCCACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCagaaaac  >  1:659104/1‑248 (MQ=255)                                                                                                                                                                                                                                       tttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCa        >  1:443997/1‑245 (MQ=255)                                                                                                                                                                                                                                       tgtcggccTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCa        <  2:443997/239‑1 (MQ=255)                                                                                                                                                                                                                                        ttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAagac  >  1:569991/1‑250 (MQ=255)                                                                                                                                                                                                                                        ttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGACAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGCAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAgacgac  >  1:791681/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            AGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGAC  >  NC_000913/2003031‑2003511

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