Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,310,770	A→G	H32H (CAT→CAC)	yciI ←	protein YciI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,310,770	0	A	G	94.4%	222.8 / ‑1.5	73	H32H (CAT→CAC)	yciI	protein YciI
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (36/32);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GAACACGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTG  >  NC_000913/1310525‑1311006                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              gAACACGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAACTGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGCTCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CACACGACCTTCATCGTGcgg                                                                                                                                                                                                                                          >  2:32259/1‑249 (MQ=255)  aacccgtttGCAAAAACCTTGCCCTTAAATTTTGAATCCCTATCCAGACCACTTTCTTAAATGGTTTCACTGAACCGTGTTCATAGACCCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGt                                                                                                                                                                                                                                         <  2:367648/244‑1 (MQ=255)        gTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGt                                                                                                                                                                                                                                   >  1:65508/1‑251 (MQ=255)                aaCCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGc                                                                                                                                                                                                                           >  1:913016/1‑251 (MQ=255)                  ccGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCt                                                                                                                                                                                                                          >  1:862360/1‑250 (MQ=255)                     tGCCCTTAAATATTGACTCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAt                                                                                                                                                                                                                      >  2:709714/1‑251 (MQ=255)                             aaTATTGAATCTCTATTCAGGACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACAACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCgg                                                                                                                                                                                                               >  2:668970/1‑250 (MQ=255)                                atTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACg                                                                                                                                                                                                            >  1:334068/1‑250 (MQ=255)                                    caTCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGGCAACCGCTAGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGAGCCGGACGAACgg                                                                                                                                                                                                       >  2:533173/2‑251 (MQ=255)                                       ctctATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATAGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCTGTAACTGTAAACGTGCTAAATGTGCCGGACGAACggca                                                                                                                                                                                                     >  1:69406/1‑248 (MQ=255)                                        tctATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTACCGGACGAACGGaa                                                                                                                                                                                                     >  2:243220/1‑249 (MQ=255)                                            ttCAGAACACTATCTTAAATGGTTTCACTGAAACGTGTTCATAAACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTTTCTACCGCTGGCATTGGACCAGCTGTCAC‑CAAACTACCTTCATCATGCAGTAACTGTAAACGTGCTACATGTGCCGGACAAACGGAAAGGCgg                                                                                                                                                                                               >  2:979429/1‑250 (MQ=255)                                                 aaCACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTc                                                                                                                                                                                          >  1:669716/1‑251 (MQ=255)                                                 aaCACTTTCTTAAAGGGTTTCACTGAAACGTGTTCATAGACCCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTc                                                                                                                                                                                          <  1:294754/251‑1 (MQ=255)                                                       ttCTTAAATGGTTTCACTGAAACGTGTTAATAGACCCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAgg                                                                                                                                                                                      <  2:771400/249‑1 (MQ=255)                                                        tCTTAAATGGTTCCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCt                                                                                                                                                                                    >  1:339723/1‑250 (MQ=255)                                                          ttAAATGGTTTCACTGACACGTGTTCATAGACTCCTTCCTCTACGTACGGGTCAGCATCCTCCCAGGCCTGAGCTGCTTCCAGCGATTCACATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGTATCGTTACTGTCTACCGCTGGCATTGGACCCGCTGTCCA‑CCAACGACCTTCATCTTGCAGTAACTGTAACCGTCCTCAATGTTCCGGCCGAACGGaacagcgcttttctagtctttc                                                                                                                                                                                 >  2:284353/1‑231 (MQ=255)                                                           tAAATGGTTTCACTGAAACGTGTTCATAGCCTCCTGCCGCTACGTACGGGTCAGCATCGGCACAGGCCTGAGCTGCTTCCAGCGATTCAACTTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAC‑CAAACGAACTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGACACTGCGTTTTTCGAGGCTATc                                                                                                                                                                                 >  2:481332/1‑250 (MQ=255)                                                              aTGGTTTCACAGAAACGAGATCATAGAATCCTGCCGCTACGAACGGGTCAGCATCGGCCCAGGACTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAAACCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAACCGACCTGAATCGTGCGGTAACTGAAAACGTGCTACATGTGCCGGCCGAACGGAAAGGCGTTTTTCGAGGTTATCaact                                                                                                                                                                             >  2:549330/1‑248 (MQ=255)                                                                  tttCACTGAAACGTGTTCATAGACTCATGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTAtgaacacggt                                                                                                                                                                         >  1:363365/1‑242 (MQ=255)                                                                         gAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAACCGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGc                                                                                                                                                                  >  2:825039/1‑251 (MQ=255)                                                                           aaCGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCCGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTACATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCgg                                                                                                                                                                >  1:498907/1‑250 (MQ=255)                                                                           aaCGTGTTAATAGACTCATGCCGCTACGTACGGGGCAGCAACGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATATCGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTGTCGAGGCTAtgaacacggtttagtgaaa                                                                                                                                                                <  2:363365/251‑19 (MQ=255)                                                                                gtttaTAGACCCCTGCCGCTACGTACGGGGCAGCATCGGCCCAGGCCTTAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAAt                                                                                                                                                            <  2:550450/246‑1 (MQ=255)                                                                                gtTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAAt                                                                                                                                                            >  1:120701/1‑250 (MQ=255)                                                                                 ttCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATa                                                                                                                                                           <  1:675313/250‑1 (MQ=255)                                                                                  tCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATAGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAAc                                                                                                                                                         >  1:702044/1‑251 (MQ=255)                                                                                       gACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGGAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATac                                                                                                                                                     >  1:549381/1‑250 (MQ=255)                                                                                                      gtacgGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTc                                                                                                                                    >  1:890722/1‑251 (MQ=255)                                                                                                         cggttCAGCATCGGCCCAGGCCTGATTTGCTTCCAGCGATTCAAATTCAGTAATAACGGTTGCGCCAGTAAATCCCGCAGCCCCTGGATCGTTCCGGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGTCCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GTCGCTACTCCtt                                                                                                                                 <  2:549381/247‑1 (MQ=255)                                                                                                              cAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTTAGCGTAAATAACATACAACACG‑GGCACAACGCCTTGTTc                                                                                                                             <  2:490183/250‑1 (MQ=255)                                                                                                               aGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCACATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCCGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGCGGCTTTCAGCTTTGTATTGAGCGTACATAACATACAACCCG‑GGCACCACTCATTGTCCgg                                                                                                                           >  2:751968/1‑251 (MQ=255)                                                                                                                 cATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAAC‑CGAGGAACAACTCCTTGTTCggg                                                                                                                          >  1:546528/1‑250 (MQ=255)                                                                                                                           ggCCTTAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTaaa                                                                                                               <  2:494136/251‑1 (MQ=255)                                                                                                                           ggCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTaaa                                                                                                               <  1:238861/251‑1 (MQ=255)                                                                                                                                gAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACg                                                                                                          >  2:380753/1‑251 (MQ=255)                                                                                                                                 aGCTGCTTCCAGCGCTTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGt                                                                                                         >  2:941152/1‑251 (MQ=255)                                                                                                                                    tgctTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTAt                                                                                                      >  1:980775/1‑251 (MQ=255)                                                                                                                                     gctTCCAGCGATTCAAATTCAGCAATAACGGTTTAGCCAGTAAATCCCTCAGCCCCTGGATCGTTACTGTCTATCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTAtt                                                                                                     <  1:583748/251‑1 (MQ=255)                                                                                                                                        tCCAGCTATTCATATTCAGCACTAACGCTTGTGCAGGTCAAACCCGCAGCCCCAGGATCGTTACTTGCTTGCGCTGGCATGTGCGCTGCT‑TCAAACAAACGACCTTCATCGTGCAGAAACTGTAAACGTGATAAATGTGCCGGCCGAACGGAATGGCGTTTTTCGAGGTGATCAGCTTTGTCTTGAGCGTAAATAACATACAACATG‑GGCACAACTCCTTGTTCGGGAAAGTTGTACAGCACGTTAtt                                                                                                     <  2:771411/248‑1 (MQ=255)                                                                                                                                        tCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAAGCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTAtt                                                                                                     <  2:661653/248‑1 (MQ=255)                                                                                                                                            gCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAgg                                                                                              <  2:273694/251‑1 (MQ=255)                                                                                                                                             cGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAggg                                                                                             >  1:737241/1‑251 (MQ=255)                                                                                                                                                  cAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGa                                                                                          >  1:264816/1‑249 (MQ=255)                                                                                                                                                    aaTTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGa                                                                                          >  1:237636/1‑247 (MQ=255)                                                                                                                                                        cAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGc                                                                                  <  1:941152/251‑1 (MQ=255)                                                                                                                                                          gCAATAACGGGTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCaa                                                                                <  2:890722/251‑1 (MQ=255)                                                                                                                                                            aaTAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACg                                                                              <  1:623068/251‑1 (MQ=255)                                                                                                                                                             aTAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCGGCAACgg                                                                             <  1:8346/251‑1 (MQ=255)                                                                                                                                                              aaaCGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGa                                                                            <  2:196375/250‑1 (MQ=255)                                                                                                                                                               aaCGGGTGAGCCAGTAAATCCCGAAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCAGCGTGAAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGATGTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGGTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGaa                                                                           <  2:769549/251‑1 (MQ=255)                                                                                                                                                                     tGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAgatg                                                                      >  2:313801/1‑250 (MQ=255)                                                                                                                                                                      gAGCCAGTAAATCCCGCAGACCCTGGATCGTTACTGGCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTGGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAgatgat                                                                    <  2:906472/251‑1 (MQ=255)                                                                                                                                                                       agtcAGTAAAGCCCGATGCCCCGGGCTCGTTACTGTCTACCGCGGGCAGTGGACCATCTGTCAA‑CAAAAGACTGTCGTCGTGCAGTAACTGTAAACGTGCTAAATGTTCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCATGAGCGTAAATAACATACAACTCG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGAACTGCAACGGAAAgatgatt                                                                   <  2:24155/248‑2 (MQ=255)                                                                                                                                                                                       aGCCCCGGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCAt                                                   <  2:528484/251‑1 (MQ=255)                                                                                                                                                                                        gCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATg                                                  <  2:65508/251‑1 (MQ=255)                                                                                                                                                                                          cccTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCa                                                <  1:179074/251‑1 (MQ=255)                                                                                                                                                                                           ccTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCAt                                               >  2:492589/1‑251 (MQ=255)                                                                                                                                                                                              ggATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATa                                              <  1:350686/249‑1 (MQ=255)                                                                                                                                                                                              ggATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCAGCGTGCAGTAACTGTAAACGTGCTAAAGGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCGTTGTTAAGGCCATGCATaa                                             <  2:980775/250‑1 (MQ=255)                                                                                                                                                                                              ggATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCGTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACGCG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATaa                                             <  1:121707/250‑1 (MQ=255)                                                                                                                                                                                                aTCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAACGGGCGACGATATGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGt                                          >  2:906308/1‑251 (MQ=255)                                                                                                                                                                                                aTCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGt                                          >  1:709774/1‑251 (MQ=255)                                                                                                                                                                                                     tACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAgg                                      <  1:71098/250‑1 (MQ=255)                                                                                                                                                                                                      aCTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAggg                                     <  1:957051/250‑1 (MQ=255)                                                                                                                                                                                                        tGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATACCATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTATTTGTTAAGGCCATGCATAAAGTAAGGGTa                                   >  1:895713/1‑250 (MQ=255)                                                                                                                                                                                                         gTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAAt                                 <  1:648208/251‑1 (MQ=255)                                                                                                                                                                                                                 cTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGAAAGGGAAATTACGCCa                          >  1:660156/1‑250 (MQ=255)                                                                                                                                                                                                                  tGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTATTTACGCCaa                         >  2:712572/1‑250 (MQ=255)                                                                                                                                                                                                                               gCTGTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTc            >  1:493957/1‑250 (MQ=255)                                                                                                                                                                                                                                  gTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGTCGAAGATCTGCAACGGAAAGATGATGTCGTTGGTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCAct         <  1:982840/250‑1 (MQ=255)                                                                                                                                                                                                                                  gTCAA‑CAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCActg        <  1:930274/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaCGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCActgatcctg  >  2:563193/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              GAACACGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAA‑CAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACG‑GGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTG  >  NC_000913/1310525‑1311006

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