Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	525,835	C→T	L859L (CTG→TTG)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	525,835	0	C	T	95.8%	146.6 / ‑3.0	48	L859L (CTG→TTG)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (15/31);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGAT  >  NC_000913/525541‑526113                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                       ggggAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACACATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAAc                                                                                                                                                                                                                                                                                   <  1:705482/298‑1 (MQ=255) gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTTGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCAAGCACTTGAAc                                                                                                                                                                                                                                                                                   >  1:594247/1‑300 (MQ=255) gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAAc                                                                                                                                                                                                                                                                                   >  1:362992/1‑300 (MQ=255)        aCCCGGAGACGGGGGAACTGCTGTGGCAGTATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTAGGCAGCATGGCAGGCAGTAATGCCGCATACAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAACCTgg                                                                                                                                                                                                                                                                            >  1:154512/1‑300 (MQ=255)        aCCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGCGAGACGGTGCGCAGCTTAGGCAGCATGGCAGGCAGTAATGCCGCCTACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTgg                                                                                                                                                                                                                                                                            >  1:129808/1‑300 (MQ=255)           cGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACNTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGtat                                                                                                                                                                                                                                                                         <  1:565964/300‑1 (MQ=255)           cGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGGCGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTTAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGtat                                                                                                                                                                                                                                                                         <  1:681830/300‑1 (MQ=255)           cGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTTAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGtat                                                                                                                                                                                                                                                                         <  1:681945/300‑1 (MQ=255)           cGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGtat                                                                                                                                                                                                                                                                         <  1:132409/300‑1 (MQ=255)           cGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGtat                                                                                                                                                                                                                                                                         <  2:659729/300‑1 (MQ=255)                aCGGGGGAACTGCTGTGGCAGCATTAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCTGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTTAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGAc                                                                                                                                                                                                                                                                    <  1:786232/300‑1 (MQ=255)                           gctgTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCAAGCACTTGAACAGCCTGGTATATGACCGTGAATACgg                                                                                                                                                                                                                                                         >  2:710844/1‑300 (MQ=255)                           gctgTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCAGAACAGCCTGCCGCAGGTGGAGTGGATGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTAAATAGCCTGGTACATGACCGTGACTaagc                                                                                                                                                                                                                                                         >  2:285808/1‑297 (MQ=21)                                        gagaCGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCAAATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAAc                                                                                                                                                                                                                                            <  1:406350/300‑1 (MQ=255)                                             gAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCTTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGa                                                                                                                                                                                                                                       <  1:344540/300‑1 (MQ=255)                                              aaaCACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCAGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATTAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGAc                                                                                                                                                                                                                                      >  2:593063/1‑300 (MQ=255)                                                   cGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGTTGCACCGTGAGACGGTGCGCAGCTTAGNCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAAAAGCCTGGTATATGACCGTGACTACGGGTGGAGTGCCAACGGCGAACTGGt                                                                                                                                                                                                                                 >  2:194620/1‑300 (MQ=255)                                                   cGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTACCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTACACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGAATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGCACAGCCTGGTATATGGATGTGACTACGGGTGGAGTGACAACGGCGACCtagt                                                                                                                                                                                                                                 >  2:236583/1‑297 (MQ=18)                                                       tACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACNGGCTGCACCGTTAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCATGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTgcgc                                                                                                                                                                                                                             <  2:184975/300‑1 (MQ=255)                                                         cAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCTGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAACATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCAt                                                                                                                                                                                                                           <  2:136207/300‑1 (MQ=255)                                                             gAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGc                                                                                                                                                                                                                       <  1:745510/300‑1 (MQ=255)                                                              aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCg                                                                                                                                                                                                                      >  1:391854/1‑300 (MQ=255)                                                                                              ctgcgagatgtgtataagagacagACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAgc                                                                                                                                                                                      >  1:60645/24‑300 (MQ=255)                                                                                                ccTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAgcgc                                                                                                                                                                                    <  1:427320/300‑1 (MQ=255)                                                                                                                 ggCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGgaga                                                                                                                                                                   <  1:442665/300‑1 (MQ=255)                                                                                                                  acagACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGagag                                                                                                                                                                  <  2:60511/297‑1 (MQ=255)                                                                                                                           tGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATTGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCAc                                                                                                                                                         <  1:161681/300‑1 (MQ=255)                                                                                                                            ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCTGGACGCCGCTGTTCGAGTAGACGCTGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCAcc                                                                                                                                                        <  2:445407/300‑1 (MQ=255)                                                                                                                                   gTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCAc                                                                                                                                                 <  1:727569/300‑1 (MQ=255)                                                                                                                                              ggAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCAAGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGAc                                                                                                                                      <  1:740269/300‑1 (MQ=255)                                                                                                                                              ggAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAGACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGAc                                                                                                                                      <  1:79082/300‑1 (MQ=255)                                                                                                                                               gAATGAAGCTGGGTGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGTAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGGCGCGGAAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGaaa                                                                                                                                     >  1:416378/1‑298 (MQ=255)                                                                                                                                                               gACGCCTCTGTTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAGCGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCACCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGc                                                                                                                     <  2:110818/300‑1 (MQ=255)                                                                                                                                                                                                     gTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACATAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGc                                                                               <  2:359893/300‑1 (MQ=25)                                                                                                                                                                                                             gTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCAcacg                                                                       >  2:443121/1‑297 (MQ=25)                                                                                                                                                                                                                cgcAGCTTCGGCAGCGTGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCTAGTTACGGAGCCAGCACTTGAACAGCCTGGTATATGACCGGGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTACGCACCCTCGCACCAGACCCGGACATCCGCATCCCGTATGCCACGGACCCTGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGAc                                                                    <  2:336849/300‑1 (MQ=11)                                                                                                                                                                                                                cgcAGCTTCGGCAGCATGTCAGTCAGTAATGCCGCATACGAACTGACCAGCACATACACCCACGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGTACCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGAc                                                                    <  2:763808/300‑1 (MQ=255)                                                                                                                                                                                                                           cagcaTGGCAGGCAGTAATTCCGCATACGAACTGACCAGCACATACACCCCCGCATGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTANGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACGGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTcac                                                         <  2:688407/300‑1 (MQ=255)                                                                                                                                                                                                                           cagcaTGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTcac                                                         <  2:145646/300‑1 (MQ=25)                                                                                                                                                                                                                           cagcaTGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTcac                                                         <  1:148716/300‑1 (MQ=25)                                                                                                                                                                                                                                     ggcagTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGTATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCACCCGGCTGCCGGACCCGGAGCTGCACCCGGACAGCACACTCACAGTGTGGCCg                                               >  2:58221/1‑300 (MQ=255)                                                                                                                                                                                                                                                cGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCAt                                    <  1:512181/300‑1 (MQ=25)                                                                                                                                                                                                                                                              ccAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATgcgc                      <  1:153076/300‑1 (MQ=17)                                                                                                                                                                                                                                                                 gCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACt                   <  1:509492/300‑1 (MQ=17)                                                                                                                                                                                                                                                                  cacaTACACCCCCGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTa                  >  1:7327/1‑300 (MQ=17)                                                                                                                                                                                                                                                                         acCCCCGCAGGCCAGTTACAGATCCAACACTTGAACAGCCTGGTATATTACTGTGACTACGGGTGGATTGACAACGGCGACCTAGTGCGAATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCTGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTAc           <  2:278698/300‑1 (MQ=11)                                                                                                                                                                                                                                                                          cccccGCAGGCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGNCTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTAcc          <  2:577155/300‑1 (MQ=17)                                                                                                                                                                                                                                                                                  ggCCAGTTACAGAGCCAGCACTTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTCCCGCCACGAt  >  2:758177/1‑300 (MQ=17)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                       GTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGAT  >  NC_000913/525541‑526113

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