Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	173,427	A→G	16.9%	R656G (AGA→GGA)	fhuB →	iron(3+)‑hydroxamate import ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	173,427	0	A	G	16.9%	130.2 / 2.1	67	R656G (AGA→GGA)	fhuB	iron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (19/33);  new base G (10/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTT  >  NC_000913/173196‑173654                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     aCCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATT‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAAC‑a                                                                                                                                                                                                                              >  1:313819/1‑239 (MQ=255)        cgaTGACTATTGGACCGTTGATTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCg                                                                                                                                                                                                                        <  1:289673/238‑1 (MQ=255)                    gACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGTGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGGGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCat                                                                                                                                                                                                           >  1:490111/1‑239 (MQ=255)                        gTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATaa                                                                                                                                                                                                       <  1:463103/239‑1 (MQ=255)                         ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGCCTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGGTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCCATTTGTTTGGAAAGC‑CAAGCCGTTAATTTTTCCAAATAAc                                                                                                                                                                                                      >  1:434656/1‑239 (MQ=255)                         ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGT‑TGCTGT‑TCACCTTTTTCTG‑TCGC‑CCACATTTTATCTGTTTTTTTTGAAAGC‑AAAACCCTTAATTTTTTCaataaat                                                                                                                                                                                                      >  2:479831/1‑233 (MQ=255)                         ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGT‑CGCGTT‑CAACCCTTTTC‑GCGCGCCCCA‑ATTTTTTCTTTTTTTTGGGAAAGC‑AAACCCGTTAAttttttccaataaa                                                                                                                                                                                                      >  2:894978/1‑230 (MQ=255)                         ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGGGGGCGGATGGGGCTGTTTCCATTCCAAGTCCCGGCGGGGG‑TG‑TGTGCAAACGTGTTT‑G‑CCGCGCCACAATTTTTTCATTTTTTGTGAAAAC‑AGAGCGTCTAATTTTTTCAAATaaa                                                                                                                                                                                                      >  2:879409/1‑238 (MQ=255)                         ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGGGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCGTGC‑GT‑CAAACTTTTTC‑G‑GGCGGCCCCTTTTTTTCTTTTTGTTTGGAAAACAAGGCCCGTTATTTTTTcaaaaaac                                                                                                                                                                                                       >  2:478067/1‑232 (MQ=255)                         ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGGGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAAATCCCGGCGGGGC‑TGCCGT‑CTACCTTTTTC‑G‑GCGCGCCATATTTTTTTTATTTGTGGGGAAACC‑AAAGCCGGTAATTTTTTCAAATaaa                                                                                                                                                                                                      >  2:860280/1‑238 (MQ=255)                           taGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACgg                                                                                                                                                                                                    <  1:490511/238‑1 (MQ=255)                                tttttgGTTTTAGGGATCCGCTTTTTTCGAGGATTATGGGCTTTCAGCGCACGATGCCGCACCTCGTAAGTTCGGCGCTGGTGGGGCGTTTTCTGCTGGTGTGCGCTGACTGGTGTGGGCGGATGGGTCTGTTTCCATTCCAGATCCCGGCGGGGC‑AGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTg                                                                                                                                                                                               <  2:740669/236‑1 (MQ=255)                                   ttgGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGG‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTTTCTATTTGTTGGGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTa                                                                                                                                                                                             >  1:59519/1‑238 (MQ=255)                                   ttgGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGTGTT‑CAACCTTTGTC‑G‑GCGCGCCATATTTTATCCTTTTTTTGTGAAAGC‑CAAGCCGTTAATTTTTCCCTTTAACGGCCCTGTa                                                                                                                                                                                             >  1:858197/1‑238 (MQ=255)                                   ttgGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCCTATTTTATCTATTTGTTGAGAAAGC‑CAAGCCGTTAATTTTTCCTTATAACGGCCCTGTa                                                                                                                                                                                             >  1:191510/1‑238 (MQ=255)                                      gTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGGGAAAGC‑AGAACCGTTAATTTTTCCATATAACGGCCCTGTAAGcc                                                                                                                                                                                         >  1:142926/1‑239 (MQ=255)                                                cccGCGTTTTGCGCGGGTTGTGGTGTTTCGACGGGCGTGGCCACACATCGTAATTTCGGGGCTGGGGGGTGGTTTACCGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑TCGTGCCATACTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATa                                                                                                                                                                               <  2:341025/238‑1 (MQ=255)                                                 ccGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATA‑G‑GCGCGCCATATTTTTTCTCTTTGTTTGGGAAGC‑AGAACCGGTAATTTTTCCCACTAACGGCCCTGTTCGCCCGGAAGGGTaa                                                                                                                                                                              >  1:302705/1‑239 (MQ=255)                                                         tGCGCGGATGGTGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa                                                                                                                                                                      <  1:649317/239‑1 (MQ=255)                                                                gttgtgGGCTTTCGACGGACGATTCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc                                                                                                                                                                <  2:711124/234‑1 (MQ=255)                                                                          ttCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CCACCTTTATC‑G‑GCGCCCCATCTTTTATCTTTTTGGTGGGCAAAC‑AGAGCCGTTAATTTTTCCATATAACCGCCCTTTACGCCCTGCCGGATAAAGCGTTTCAGCCACAAAAGAGGccc                                                                                                                                                     >  2:272047/1‑239 (MQ=255)                                                                                                cGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTTTC‑G‑GCCCGCCATATTTTATCTATTTGTTGGGAAATC‑AAAGGCGGTAATTTTTCCCTATTACGGCCCTGTACGCCTGGCCGGATAAAACTGTTAGCCGACAACCGATGCCCCACGGCGCACCGGCGCGTCTTa                                                                                                                               >  2:441835/1‑239 (MQ=255)                                                                                                cGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCCTATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGGTAAAGCGTTTCCGCGTCAACGGATGCCCGGTGCGACGCTGGCGCGTCTTa                                                                                                                               >  1:376780/1‑239 (MQ=255)                                                                                                cGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATCTTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAACGTTTAACCGACAACGGATGCCCGATCCGCCCCGGGCCCGCCTTa                                                                                                                               >  1:824182/1‑239 (MQ=255)                                                                                                  tAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTTAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGGCAACGGATGCCCGATGCGACGCTGGCGCGTCTTATc                                                                                                                             >  1:426715/1‑239 (MQ=255)                                                                                                       tCGGCGCGGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTTTGGGCGGATGGTGCTTTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGcc                                                                                                                        <  1:269541/239‑1 (MQ=255)                                                                                                        cGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGCCTTTTCGGGCCt                                                                                                                       >  1:445007/1‑239 (MQ=255)                                                                                                             gggcaagggttttttCTGCTTGTTTTCCCCGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGGTCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACaaa                                                                                                                  <  2:220513/228‑1 (MQ=255)                                                                                                                gtgggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGGTAAAGCGTTTTCCCGACAACGGATGCCCGATGCGACCCTTGCGCGTCTTATCGGGCCCACaaaac                                                                                                                >  1:704438/1‑236 (MQ=255)                                                                                                                          tACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTTGAGGGATAAAGCGTTTACCCTACAACGGATGCCCGATTCGACCCGGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAAACg                                                                                                     >  2:147057/1‑239 (MQ=255)                                                                                                                          tACTGCTGGTGTTCGATGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCTGCGGGGG‑TCCCTG‑CAACCTTTTTT‑G‑GCGCGCCCGATTTTATTTATTTGTTCAGAAAAC‑AGAGCCGTTAATTTTTTCATATATACGCCCTCTTGGCCCGGCCGGTTAAACCGTTTACGCGACAAGGGGTGCCCTCTGTGACCCCTGGGCGCCCTTTAGGGCCTACAAACCACCCCCaaaca                                                                                                     >  2:686543/1‑236 (MQ=255)                                                                                                                               ctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATctgtc                                                                                                  >  1:220513/1‑233 (MQ=255)                                                                                                                                    gtTGGCTGAATGGTGGGGGCGGATTGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAt                                                                                           <  1:366401/239‑1 (MQ=255)                                                                                                                                       cGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTTGGAAAGC‑AGAGCCGTTAATTTTTTCATATAACGGCCCTTTACGCCTGGACGGACAAAGCGTTTAACCGACACCGGCTGCCCGATGCGCCGCCTGCGCGGCCTTTCGGGCCCCCAAACCGCCCCGCAACGGACGTCgggaaaa                                                                                        >  2:10195/1‑234 (MQ=255)                                                                                                                                                 gttgGGCCGATGGTGCGTTGTCCCTTCCAGGTCCGGGCGGGGC‑TGCTGT‑TAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGGGCCGTTAATTGTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcg                                                                               <  1:28284/236‑1 (MQ=255)                                                                                                                                                    gggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACCCTGGCGCGTCTTATCGGGCCCACAAACAGCCCCGAATCGTAGGTTGGCTAACGCGTTCCcgccgc                                                                           >  1:153441/1‑239 (MQ=255)                                                                                                                                                           tggttttGTTTCCTTTCCAGACCCGGGGGGGGT‑TCCGGT‑ACCCCTTTATG‑G‑GTGGGCAATATTTTTTCTATTTGTTGAAAAGGC‑AGAGCCGTTAATTGTTCCATAAACCGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                                                                     <  2:824514/232‑1 (MQ=255)                                                                                                                                                           tGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCCAAACCTAGGGCGGATAAGGCGTTCCCGCCGCACCCGAc                                                                    >  1:824514/1‑239 (MQ=255)                                                                                                                                                              tGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAAttt                                                                                                                                                                                                                 <  1:471497/95‑1 (MQ=255)                                                                                                                                                               gCTTTTCCCATTCCAGATCCCGGCGGGGT‑TGCTGG‑CAACCTTTATC‑G‑GCGCGACATATTTTATCTATTGGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa                                                                 <  1:574392/238‑1 (MQ=255)                                                                                                                                                                    ttcctttccaGTTCCGGGGGGGGT‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGc                                                           <  1:11358/239‑1 (MQ=255)                                                                                                                                                                    tcacattccaGAGCCGGGGGGGGC‑TGCTGT‑AAACCTTTAGC‑G‑GCGCGCCATATTTTATCTTTTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCGGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGc                                                           <  1:415997/236‑1 (MQ=255)                                                                                                                                                                        attccaGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑C‑GCGCGCCATATTTTATTCATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATaa                                                                                                                                                                                                       >  2:471497/1‑95 (MQ=255)                                                                                                                                                                                   cggcggGCC‑TGCTGT‑CAACCTTAAAC‑G‑GCGGGCAATATTTTACTTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCGGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCg                                            <  1:479831/239‑1 (MQ=255)                                                                                                                                                                                      gggggC‑TGCTGG‑CAACCTTTATC‑G‑GCGCGCCATTTTTTATTTATTTGTTGAGAAAGC‑AGGGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCa                                          <  1:147057/237‑1 (MQ=255)                                                                                                                                                                                              gctgT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTTTTTGTTGAGAAAGC‑AGAGCCTTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCg                                  <  1:663550/239‑1 (MQ=255)                                                                                                                                                                                              gctgT‑CAACCTTTAAC‑G‑GCGCGTCATTTTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCg                                  <  1:514820/239‑1 (MQ=255)                                                                                                                                                                                              gctgT‑AACCTTTTACC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCg                                  <  1:459742/239‑1 (MQ=255)                                                                                                                                                                                              gcgtt‑CAACTTTTATC‑G‑GCGTGCCATTTTTTTTCTATTTGGTAAGAAACA‑CAAGCCTTTAATTTTTCCATATACAGGCCTGGTACGACTGGAGAGATCAAGCGTTTACGCGACAACGGATGCCCGATGCGACGGTTGGGCGTCTTTTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTCCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCg                                  <  2:199032/235‑1 (MQ=255)                                                                                                                                                                                                      aCCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGCCGCTGGCGCGACTTATCGGGCCTACAACCGGCCCCGAATCGTCGGCCGGATAAGGCGTTCACCCCCCACCCAACACACCATACGGGTCAACACCTTAGCAAAAAccccccccgc                            >  2:103822/1‑231 (MQ=255)                                                                                                                                                                                                               tcggcagcctatatTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTCCGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGa                    <  1:700405/229‑1 (MQ=255)                                                                                                                                                                                                                aaaagagacaaatTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTCCCATAGAACGGCCCTGTACGCCTGGACGGATAAACCTTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGa                    <  1:285641/227‑1 (MQ=255)                                                                                                                                                                                                                   cgcgcCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGcc                                                                                                                                                      <  2:879618/105‑1 (MQ=255)                                                                                                                                                                                                                   cgcgcCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGcc                                                                                                                                                      >  1:879618/1‑105 (MQ=255)                                                                                                                                                                                                                   cgcgcCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGCCAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCCTCGATggg                >  1:551116/1‑238 (MQ=255)                                                                                                                                                                                                                    gcgcCATATTTTATCTATTTGTTGAAAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGgt                <  1:894978/238‑1 (MQ=255)                                                                                                                                                                                                                           tttttAT‑TATTTGTTGGGAAAGG‑AGGCCCGTTAATTTTCCCATATAACGGCCCTGGACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  1:421651/238‑1 (MQ=255)                                                                                                                                                                                                                            tttttttttttttttGGGAAAGC‑AGAGCCGTTAAATTTTCCAAATAACGGCCCGGGCCGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTAAAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  2:704438/230‑1 (MQ=255)                                                                                                                                                                                                                            tttttttttttttttGAGAAAGC‑AGGGCCGTTTATTTTTCCATATACCGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGGGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  2:142926/230‑1 (MQ=255)                                                                                                                                                                                                                            ttttttCTATTTGTTGAGAAAGC‑AGGGCCGGTAATTTTTCCATAAACCGGACCTGTACGTCTGGCAGTATAAAACGTTTCCGCCGCAAGGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  2:490111/239‑1 (MQ=255)                                                                                                                                                                                                                            ttttATTTATTTGTTGAGAAAGC‑AGGGCCTTTAATTTTCCCATAGAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  1:478067/239‑1 (MQ=255)                                                                                                                                                                                                                            ttttATCTTTTTGTGAAGAAAAC‑AGAGCTTTTAATTTTTCCATAAAACGGCCCGGTGCGCGGGGACGGTAAAAGCGTTTACGCGAAAACGGATCCCGGATGGGACGCGGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  2:551116/239‑1 (MQ=255)                                                                                                                                                                                                                            ttttAACTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat       <  1:493018/239‑1 (MQ=255)                                                                                                                                                                                                                              tttttttttttttGAGAAAGG‑AGTGCCGTTATTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTCCGCGACAACGGATGCCCGGTGCGACGCTGGCGGGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGatat     <  2:426715/232‑1 (MQ=255)                                                                                                                                                                                                                                  cTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGGTCACGCCCCATCCGACAAACCATGCTGGGTCACAACTTCCCAAACCCCCGCCCCGCCACAACGAGGGAGTTATGGAAATCtt  >  2:457532/1‑238 (MQ=255)                                                                                                                                                                                                                                      ttgttgAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGa                    >  2:141526/1‑216 (MQ=255)                                                                                                                                                                                                                                      ttgttgAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGa                    <  1:141526/216‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     ACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGC‑TGCTGT‑CAACCTTTATC‑G‑GCGCGCCATATTTTATCTATTTGTTGAGAAAGC‑AGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTT  >  NC_000913/173196‑173654

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