Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,464	A→T	E54V (GAG→GTG)	fimC →	type 1 fimbriae periplasmic chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,464	0	A	T	95.5%	134.3 / ‑3.9	45	E54V (GAG→GTG)	fimC	type 1 fimbriae periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (18/24);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTC  >  NC_000913/4544219‑4544686                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                tGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAAc                                                                                                                                                                                                                            >  1:5455/1‑251 (MQ=255)      gcATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAAGTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGcc                                                                                                                                                                                                                       <  1:335600/251‑1 (MQ=255)       cATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTCGGAGCGACTCGCGTAATTTATCCGGCGGGGCAAAAACAAAGGCACCTTGccc                                                                                                                                                                                                                      >  2:329/1‑250 (MQ=255)          cGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTg                                                                                                                                                                                                                    <  1:329/250‑1 (MQ=255)          cGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGa                                                                                                                                                                                                                   <  1:49369/251‑1 (MQ=255)               aTGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGCTGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCACAAACAAGTGCAACTTGCCGTGACAa‑at                                                                                                                                                                                                             >  1:441753/1‑251 (MQ=255)                 gCCCAGGCCTGGGTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATATCAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAa‑ata                                                                                                                                                                                                            <  2:223483/250‑1 (MQ=255)                   ccAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGAGAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAatg                                                                                                                                                                                                         <  2:35960/251‑1 (MQ=255)                          tGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGCAGCGGGAGTGGCCTTAGGTGCTACTCGCGTAATTTATCCGGCAGGGAAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGaaa                                                                                                                                                                                                   >  2:392717/1‑250 (MQ=255)                                 ctTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTa                                                                                                                                                                                             <  1:116932/249‑1 (MQ=255)                                    tAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTa                                                                                                                                                                                         <  2:97061/250‑1 (MQ=255)                                      aCCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTAGCCTGATGTTCATGGCAAGGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGAGGAAAATAGTACCTAttt                                                                                                                                                                                      <  2:272800/251‑1 (MQ=255)                                           tCAGTAAATGTTCAGCAGATAATGTATTAACTGGAACAGGATAATGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTAGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTTCCGGACGCGCTGAAGCGAGCGTTGCATTAGCTGCGACTCGCGTAATTTATCCGGATGGGCAAAAACAAGTGCAACTGGACGTGACAA‑ATAATGAAGAAAATAGTACCCATTTAATTc                                                                                                                                                                                 >  2:474579/1‑251 (MQ=255)                                             aGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCaa                                                                                                                                                                               <  2:5455/251‑1 (MQ=255)                                                 aTTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACAGTCTGCTTGCTGGCAGGTATCCTGATGTTCAGGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCAt                                                                                                                                                                           <  1:483041/251‑1 (MQ=255)                                                                aTGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATc                                                                                                                                                                             <  1:328034/234‑1 (MQ=255)                                                                aTGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATGTAATTCAATc                                                                                                                                                                             >  2:328034/1‑234 (MQ=255)                                                                    gATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATgg                                                                                                                                                        >  1:360962/1‑251 (MQ=255)                                                                       aacaggaacaggAAAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGATTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTACTTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGgtgt                                                                                                                                                     >  2:60924/1‑251 (MQ=255)                                                                                     aGTGAATAATAAAAACGTCAAAGTAAGGAAAACGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGtt                                                                                                                                       <  2:38086/251‑1 (MQ=255)                                                                                          gTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATGTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATc                                                                                                                                  <  1:235795/251‑1 (MQ=255)                                                                                               aaaaaCGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGaa                                                                                                                             <  2:360962/251‑2 (MQ=255)                                                                                                aaaaCGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACg                                                                                                                            <  1:187196/251‑1 (MQ=255)                                                                                                 aaaCGTCAATGTAAGGAAATCGCAGGAAATAACATACTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTAGACGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGATGTAAGGGATGGTCGTTTTATCGTGACGc                                                                                                                           >  1:300429/1‑251 (MQ=255)                                                                                                    cGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGcctc                                                                                                                        >  1:167253/1‑251 (MQ=255)                                                                                                                   cgacgctcttccgatcgtCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGAt                                                                                                           <  2:355423/232‑1 (MQ=255)                                                                                                                        aGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTGTGCGATGAAGGGa                                                                                                    <  2:441753/251‑1 (MQ=255)                                                                                                                           aaaTAACATTCTGCTTGCTGGCAGGTAGCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaaaa                                                                                                 <  2:321880/251‑1 (MQ=255)                                                                                                                               aaCATTCTGCTTGATGGCAGGTATCCTGATGTTCATGGCAATGATGTTTGCCGGGCGCGCTGAAGCGGGAGTGGCCTTAGGTGCGAATCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGT‑ACAACATAATGATGAAACTAGTACCTATTTAATTCAATCATGGGGGAAAAATGACGATGGGGTAATAGATGGTCGTTTTATCGTGCCGCCTCCTCTGTTTGCGTTGAAGGGaaacaa                                                                                               >  2:36940/1‑246 (MQ=255)                                                                                                                                    tCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGAt                                                                                                           >  1:355423/1‑232 (MQ=255)                                                                                                                                       gcttgctGGCAGTTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGGTTTATCGTGACGCATCCTCTGTTTTCGATGAAGGGATAAAAATAGAAt                                                                                        <  2:434178/248‑1 (MQ=255)                                                                                                                                        cttgctGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTAGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCt                                                                                    >  1:321218/1‑251 (MQ=255)                                                                                                                                         gtgctGGCAGGGATCCTTATGTACATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCGTAGGTGCGACTCGCGTAATGTATCTGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑AAAATGATGAAAATAGTACCTATTTAATTCAATCCTGGGTGGAAAATGCCGATGGTGGAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGACAATACCtt                                                                                   <  2:267649/250‑1 (MQ=255)                                                                                                                                           gctGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTa                                                                                  <  1:259421/250‑1 (MQ=255)                                                                                                                                            ctGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCttccg                                                                                >  2:377680/1‑248 (MQ=255)                                                                                                                                             tGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAAt                                                                                                                                                                                   >  2:402624/1‑151 (MQ=255)                                                                                                                                             tGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAAt                                                                                                                                                                                   <  1:402624/151‑1 (MQ=255)                                                                                                                                                                 aTGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATaa                                                           <  2:281525/251‑1 (MQ=255)                                                                                                                                                                                 ccGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAc                                           <  1:323984/251‑1 (MQ=255)                                                                                                                                                                                        cgcTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGCCGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAACTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaa                                     >  2:271173/1‑250 (MQ=255)                                                                                                                                                                                                                     cTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGAtt       >  2:187872/1‑251 (MQ=255)                                                                                                                                                                                                                        gcgTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGttt                                                                                                                <  1:1504/143‑1 (MQ=255)                                                                                                                                                                                                                        gcgTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGttt                                                                                                                >  2:1504/1‑143 (MQ=255)                                                                                                                                                                                                                          gTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc  >  1:179054/1‑251 (MQ=255)                                                                                                                                                                                                                           tAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTATGGATGAACGTTAAAGCGATTCCGTc  >  2:327093/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                TGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAA‑ATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTC  >  NC_000913/4544219‑4544686

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