Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,955,304	C→T	R609R (CGC→CGT)	ilvD →	dihydroxy‑acid dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,955,304	0	C	T	94.0%	194.5 / 1.5	69	R609R (CGC→CGT)	ilvD	dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (27/36);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGCGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCG  >  NC_000913/3955059‑3955543                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 gTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGGGATa                                                                                                                                                                                                                                              >  2:94946/1‑250 (MQ=255)      aCCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCg                                                                                                                                                                                                                                        >  1:62169/1‑250 (MQ=255)             gcggcaagcggcGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGACCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGAGCGATA‑AATCGAAACTggg                                                                                                                                                                                                                                <  1:143656/251‑1 (MQ=255)                   agcggcGGCAGCAGTTGCCTGAATGAAGATGGCGACCTGATCGCTACCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGACAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AAACGAAACTGGGGGGTta                                                                                                                                                                                                                          <  2:132223/251‑1 (MQ=255)                        cggtagcaTGGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATGCAGCTACAGGTAAGCGATGCCGAACTGGCGGCGAGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTCCGTGCTTATGCCTGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGGGTGATA‑AATCGAAACTTGGGGGTTaataat                                                                                                                                                                                                                     <  1:136968/247‑1 (MQ=255)                           cagcaTTTGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCGCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGc                                                                                                                                                                                                                  <  1:468387/251‑1 (MQ=255)                            agctttGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGTCTTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTGGCCCTGCGTGCTTATGCCAGGCGGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAAGGGCt                                                                                                                                                                                                                 <  2:307512/247‑1 (MQ=255)                               aTTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGGGATA‑AATCGCACCTGGGGGGTTAATAATGGCTGAc                                                                                                                                                                                                              >  2:32612/1‑251 (MQ=255)                                            aaGATGGAGACATGATCGCTATCTAAATCCCGAACCGTGGATTTCAGTTACAGGAAAGCGATGCCGACCTTGCGGCGCGTCGTGAAGCGCAAGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGGGAACGTAAGGGCTCCTTTGCCCTGCGTGCTTATGCCCGCCTGGCAACAAGCGCCGAAAAAGACGCGGTGCGC‑ACACATTCGAAAATAGAGAGGTAATAATGGATGACTCGCAAccccagt                                                                                                                                                                                                 >  2:153792/1‑248 (MQ=255)                                             aGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGAGGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTc                                                                                                                                                                                                <  1:105588/251‑1 (MQ=255)                                               aTGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCg                                                                                                                                                                                              >  1:459198/1‑251 (MQ=255)                                                 ggTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCgg                                                                                                                                                                                             <  1:274665/250‑1 (MQ=255)                                                  gTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGAGTCGAAACCCCTGTCCgg                                                                                                                                                                                             >  2:422734/1‑249 (MQ=255)                                                  gTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGGTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTg                                                                                                                                                                                           <  2:79738/251‑1 (MQ=255)                                                      ccTGATCGCTGTCGACATCCCGAACCGTGGCATGCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGACTCCTTTGCCCTTCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTcc                                                                                                                                                                                       <  2:419138/251‑1 (MQ=255)                                                      ccTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTcc                                                                                                                                                                                       <  2:419144/251‑1 (MQ=255)                                                       cTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCg                                                                                                                                                                                      >  1:306853/1‑251 (MQ=255)                                                           tCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGa                                                                                                                                                                                    >  1:184007/1‑249 (MQ=255)                                                            cGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGGGCTCCGGAAgg                                                                                                                                                                                 >  2:73054/1‑251 (MQ=255)                                                                aTCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGcc                                                                                                                                                                             >  1:48170/1‑251 (MQ=255)                                                                 tCGACATCCCGAACCGTGGCATTCAGTTACAGGCAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGACCGTCAGGTCTCCTTTACCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGCGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGAATAGAAACCCCTGTCCGGTGCTCCGGagggggacg                                                                                                                                                                            >  2:219971/1‑243 (MQ=255)                                                                               cGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGCCGTGAAGCGCAGGACGCTCGAGGTCACAAAGCCTGGACGCCGAAAAACCGTGAACGTCAGGTCGCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCACTCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAGCCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCa                                                                                                                                                              <  1:404800/251‑1 (MQ=255)                                                                                  ggCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGt                                                                                                                                                            <  1:387102/250‑1 (MQ=255)                                                                                   gCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGtgc                                                                                                                                                          >  1:311398/1‑251 (MQ=255)                                                                                     aTTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGtgctg                                                                                                                                                        <  1:92325/251‑1 (MQ=255)                                                                                     aTTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGtgctg                                                                                                                                                        <  1:209916/251‑1 (MQ=255)                                                                                            tACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGGGCTCCGGAAGGGGCCGAATATTTAAGAGCAGTGCTgcgcgcgc                                                                                                                                                 >  2:93458/1‑251 (MQ=255)                                                                                             aCAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAACGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAACGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGTAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGcc                                                                                                                                                >  2:458112/1‑251 (MQ=255)                                                                                                    aGCGAAGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTAc                                                                                                                                         <  2:371802/251‑1 (MQ=255)                                                                                                     gCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACg                                                                                                                                        >  1:318913/1‑251 (MQ=255)                                                                                                       gATGCCGAACTGGCGGCGCGTCGTGACGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACg                                                                                                                                        >  1:308778/1‑249 (MQ=255)                                                                                                             gAACTGTCTGCGCGGCGTGAAGCGTAGGACGCTCGAGGTGACAAAGCCTGGACGCAGAAAAAGCGCGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGTGTTGCGTGATC‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGTTGCCGAATATTTAATAGCAGTGCTGCGCGCGCCGGTTTACGAGgcggcg                                                                                                                                <  1:136944/251‑1 (MQ=255)                                                                                                                 gtggggcgaggcgtGAAGCGCAGGACGCTCGAGGTGACCAAGCCTGGACGCCGACAAATCGTGAACGGCAGGGCTCCTGGGCCCTGCGTGCTTATGCCAGCATGGCAACCAGCGCTGACAAATGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCACCCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATGTAAGAGCAGTGAAGCGCGCGCCGGTGTACGAGGCGGCGCAgg                                                                                                                            <  2:52011/240‑1 (MQ=255)                                                                                                                   gcggcgCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGtt                                                                                                                          <  1:416077/251‑1 (MQ=255)                                                                                                                     ggcgCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTAc                                                                                                                        <  1:35455/251‑1 (MQ=255)                                                                                                                     ggccgtgcgtGAAGCGCAGGACGTTCGAGGTTGCAAACCCTTGGCGCCTAAAAACCGTGAACGTCAAGCCTTCATTGCCCTGCGTGCTTATGCCAGACTGGCAACCAGCGCAGGCAAAGGCGCGGTGCGTGATA‑AAGCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCTCCGGTTTACGAGGCGGCGCAGGTTAc                                                                                                                        <  2:223027/244‑1 (MQ=255)                                                                                                                         cgtcgtGAAGCGCAGGACGCTCGGGGTGCCAATGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGCAACTGGGGAGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCAGGAATGAGCCGTATATTTTAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAggctaagacc                                                                                                                    >  2:38023/1‑243 (MQ=255)                                                                                                                            cgtGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCt                                                                                                                  <  1:217345/250‑1 (MQ=255)                                                                                                                               gAAGCGCAGGCCGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTACGGTGCTCCGGAAGGAGCCGAATCTTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCGGGTTAGGCCGCTAcc                                                                                                               >  2:222470/1‑249 (MQ=255)                                                                                                                                      aGGACGCTCGGGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATgg                                                                                                       <  2:444926/251‑1 (MQ=255)                                                                                                                                         aCGCTCGAGGTGACAAAGCCTGGTCGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAATAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGaaa                                                                                                    <  1:64127/251‑1 (MQ=255)                                                                                                                                          cGCTCGAGGTGACAAAGCATGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGACGACAAAGGCGCGGTGCGTGATC‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGGGCCGCATATTTCAGAGCAGTTCTGCGCGCGCCGGGTTACGAGGCGTCGCCGGTTACGCCGAAACAAAAACTGGaaa                                                                                                    >  2:210152/1‑250 (MQ=255)                                                                                                                                           gCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGATGACTCGCAACACCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTAAGCCGCTACAAAAAATGGaaaa                                                                                                   >  2:3078/1‑250 (MQ=255)                                                                                                                                             tCGAGGTGAAAAAGCCTGGACGCCGAAACATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGAAAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAAc                                                                                                 <  2:76453/250‑1 (MQ=255)                                                                                                                                                   tGACAAAGCCTGGACTCAGAAAAATCGTGAACGTCTGGTCTCCTTTGCCCTGCGTGCTAATGCAAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACCGGGGGGTTAATAATGGATGAATCGCAACCCCTGACCGGTGCTCCGGAAGGAGCCGAATATTTAAGAGCAGTGGTGCGCGCGCCGGTTTACGCGGCGGCGCCGGTTAAGTCGCTACcaaaacagggaaaacaggcggc                                                                                          >  2:124156/1‑234 (MQ=255)                                                                                                                                                      cAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTcgcg                                                                                       <  1:93458/251‑1 (MQ=255)                                                                                                                                                        aaGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTc                                                                                     <  1:377056/251‑1 (MQ=255)                                                                                                                                                              ggACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGa                                                                                 <  1:323943/249‑1 (MQ=255)                                                                                                                                                                  gCCGAAGAATCGTGAACGTCAGGTCTCCTTTGCCCTTCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGt                                                                           <  1:214535/251‑1 (MQ=255)                                                                                                                                                                        aaaTCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCt                                                                     >  1:142969/1‑251 (MQ=255)                                                                                                                                                                               gAACGTCAGGTCTCCTTTGACCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGATGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAg                                                              <  2:349479/251‑1 (MQ=255)                                                                                                                                                                                aaCGTCAGGTCTCCTTTGCACTGCGTGCTTATGCCGGCCAGGCAACCAGAGCCGCCAAAGGCGCGGTGCGTGATA‑CATCGAAACTGGGGGGTTAATAATGGCTGACTCGCACCCCCTGTCCGGTGCTCCGGAAGGTGCCGCATATTTCAGAGCGGTGCTGCGCGCGCCGGATTACGCGGCGGCGCAGGCTACGCCGCTACAAAAAATAGAAAACCTGTCGTCGCGGCTTGTTAACGTCATTCTGGTGGAgc                                                             >  2:458339/1‑251 (MQ=255)                                                                                                                                                                                   gTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGCCAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCTCGCGCCGGTTTACGAGGCGGCGCCGGTTACGCCGCTACAAAAAATGGAACAACTGTCGTCTCGTCTTGATCACATCATTCTTGGGAAGcgcg                                                          >  2:69816/1‑251 (MQ=255)                                                                                                                                                                                     cAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCCACCCCTGTCCGGTGCTCCGGAGGGTGCCGAATATTTAGGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAGATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGACGCGCGa                                                         >  2:195700/1‑250 (MQ=255)                                                                                                                                                                                     cAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGa                                                         >  1:298881/1‑250 (MQ=255)                                                                                                                                                                                         tctcCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGa                                                         >  1:269379/1‑246 (MQ=255)                                                                                                                                                                                           cccTTTTCCCTGCGTGCTTAGGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGGGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGGGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTTAAGCGCGTAGATCgc                                                  <  2:454990/250‑1 (MQ=255)                                                                                                                                                                                            cctgtGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCg                                                   <  1:458112/245‑1 (MQ=255)                                                                                                                                                                                            ccTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCg                                                   <  1:94946/249‑1 (MQ=255)                                                                                                                                                                                                 gCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccagcca                                            <  2:166450/251‑1 (MQ=255)                                                                                                                                                                                                 gCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccagcc                                             >  1:229270/1‑250 (MQ=255)                                                                                                                                                                                                       ggTGCTTATGCCAGCCGGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACGCGCAACCCCTGTCCGGGGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGcac                                      <  2:62169/250‑1 (MQ=255)                                                                                                                                                                                                         tGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGCATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCACCAAACAATGGAAAAACTGTCGTCGCTTCTTGTAAACGTCATTATGGAGAAggggtcagatcgcccaccatagacc                                      >  2:9293/1‑225 (MQ=255)                                                                                                                                                                                                         ggCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGAGTTATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGGCCGGTGCTCCGGAAGGTGCCGAAGATTTAAGAGCAGGGCTGCGCGCGCCGGTTTACGATGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAg                                    <  2:308778/250‑1 (MQ=255)                                                                                                                                                                                                          gCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGc                                   >  1:346613/1‑251 (MQ=255)                                                                                                                                                                                                                 ccAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGCGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGc                            >  1:139245/1‑251 (MQ=255)                                                                                                                                                                                                                        ggCAACCAGCGCCGACAAAGGCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgcgc                       <  2:163059/249‑1 (MQ=255)                                                                                                                                                                                                                                         aaGGGGCGGTGGGTGATA‑AATCGAAACTGGGGGGTTAATAATGTCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCGGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGGCGCGTCTGGATAACGTCATTATGGTGAAGCGCGAAGATCGCCAGCCAGGGCACAGCGTTAAGCGGCGCGGCGCATACGCCATGATgg    <  2:44920/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggCGCGGTGCGTGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGACCTTCTGGTGAAGCGCGAAGATCACCAGACAGTGCACCGCTTTCAGCTGAGCGGCGCATACGCCATGATGGCg  >  1:153010/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 GTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGCGATA‑AATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCG  >  NC_000913/3955059‑3955543

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