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	97.9%	171.8 / NA	48	R609R (CGC→CGT)	ilvD	dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/26);  minor base G (1/0);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGCGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCC  >  NC_000913/3955071‑3955547                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    gcggcaagcggcGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGGGATAAATCGAAACTggg                                                                                                                                                                                                                                    >  2:112819/1‑251 (MQ=255)  cggcaagcggcGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTgggg                                                                                                                                                                                                                                   >  1:26663/1‑251 (MQ=255)   ggcaagcggcGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTggggg                                                                                                                                                                                                                                  >  1:131446/1‑251 (MQ=255)     caagcggcGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGt                                                                                                                                                                                                                                <  2:79287/251‑1 (MQ=255)        gcggcgGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTtaa                                                                                                                                                                                                                             <  1:110893/251‑1 (MQ=255)            cggcAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTaataat                                                                                                                                                                                                                         <  2:103373/251‑1 (MQ=255)             ggcAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATg                                                                                                                                                                                                                        >  1:76238/1‑251 (MQ=255)               cagcaTTGGCCTGCTTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAACGCCTGGACGCCGAACAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAACGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATACTGGc                                                                                                                                                                                                                      >  2:331364/1‑251 (MQ=255)                      ggCCTGATTGAAGATGGTGACCTGATCGCTATCGCCATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCGGGACGCCGAAAAATCGTGAACTTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTTACTCg                                                                                                                                                                                                               <  2:293465/251‑1 (MQ=255)                      ggCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTc                                                                                                                                                                                                                <  2:170591/250‑1 (MQ=255)                               gAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTg                                                                                                                                                                                                      <  1:70230/251‑1 (MQ=255)                                       tGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGc                                                                                                                                                                                              >  1:360314/1‑251 (MQ=255)                                            tGATCGCTATCGAAAGCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCgg                                                                                                                                                                                         >  1:227365/1‑251 (MQ=255)                                            tGATCGCTATAGACATCCCGAACTGTGGCATTCAGCTACAGGTAAGCGATGCCGACCTGGCGGCGCGTCCTGACGCGCAGGACGCTCGAGGTGACAACGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGAGCTTATGCCCGCCTGGCACCCAGCGCCGACAAAGTCGCGGCGCGTGCTAAATCGAAAATGGGGGGTTAATAATGGCTGACTCGCATCCCCTGTCCGGCGCTCCgg                                                                                                                                                                                         >  2:238692/1‑251 (MQ=255)                                             gATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGa                                                                                                                                                                                        <  1:305001/251‑1 (MQ=255)                                                              cgcaCCGTGGCACTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGCGACAAAGCCTGGACGACGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGCACGGTGCTCCGGACGGAGCCGAATATTTa                                                                                                                                                                        >  2:168523/4‑250 (MQ=255)                                                                aaCCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCCCCGGAAGGTGCCGAATATTTa                                                                                                                                                                        >  1:185681/1‑248 (MQ=255)                                                                    gTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAg                                                                                                                                                                 <  2:113430/251‑1 (MQ=255)                                                                           tCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGtgctgc                                                                                                                                                           <  1:110904/250‑1 (MQ=255)                                                                            cAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTgcg                                                                                                                                                          <  2:322711/250‑1 (MQ=255)                                                                               ttACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGGGCCGAATATTTAAGAGCAGTGCTgcgcgcg                                                                                                                                                      >  2:66009/1‑251 (MQ=255)                                                                                       aaGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGttt                                                                                                                                               >  2:304882/1‑250 (MQ=255)                                                                                       aaGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTa                                                                                                                                              <  2:26663/251‑1 (MQ=255)                                                                                              gCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAggcg                                                                                                                                       >  1:107948/1‑251 (MQ=255)                                                                                                      ggcggcGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGATGCGGCGCAg                                                                                                                                 <  1:277516/249‑1 (MQ=255)                                                                                                            gcgTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTAcgcc                                                                                                                         <  2:116782/251‑1 (MQ=255)                                                                                                                cgtGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTa                                                                                                                     <  1:71570/251‑1 (MQ=255)                                                                                                                  tGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTAc                                                                                                                    >  1:2828/1‑250 (MQ=255)                                                                                                                     aGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACaaaa                                                                                                                <  1:114714/251‑1 (MQ=255)                                                                                                                                      gTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTgtcgt                                                                                               <  1:119335/251‑1 (MQ=255)                                                                                                                                           aaaGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGt                                                                                          >  2:122487/1‑251 (MQ=255)                                                                                                                                            aaGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTc                                                                                         <  1:87776/251‑1 (MQ=255)                                                                                                                                              gCCTGGACTCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACGGTCGTCGCGTCtt                                                                                       <  2:300003/251‑1 (MQ=255)                                                                                                                                                cTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTg                                                                                      >  2:216031/1‑250 (MQ=255)                                                                                                                                                 tGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGAt                                                                                    >  1:344406/1‑251 (MQ=255)                                                                                                                                                    aCGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAAc                                                                                 <  1:124241/251‑1 (MQ=255)                                                                                                                                                          aaaaaTCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACTAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCAtt                                                                           <  1:23615/251‑1 (MQ=255)                                                                                                                                                          aaaaaTCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCAtt                                                                           <  2:70672/251‑1 (MQ=255)                                                                                                                                                           aaaaTCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTc                                                                          <  2:2828/251‑1 (MQ=255)                                                                                                                                                                           ggTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGa                                                          <  2:353129/251‑1 (MQ=255)                                                                                                                                                                            gTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGa                                                             >  1:290717/1‑247 (MQ=255)                                                                                                                                                                                    tGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGAGCgccagc                                                  <  1:317400/250‑1 (MQ=255)                                                                                                                                                                                          gCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGGCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGca                                           <  1:355976/251‑1 (MQ=255)                                                                                                                                                                                           cGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGCAGCGCGCAGATCGCCAGCCAGTGca                                           >  2:33017/1‑250 (MQ=255)                                                                                                                                                                                           cGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGcac                                          >  1:250681/1‑251 (MQ=255)                                                                                                                                                                                                              cAACCAGCGCCGACAAAGGCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcggcg                       <  1:135524/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCGCGGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCg      >  1:150241/1‑251 (MQ=255)                                                                                                                                                                                                                                   cgGTGCGTGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGcc  >  2:336587/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    GCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCGTGGCATTCAGTTACAGGTAAGCGATGCCGAACTGGCGGCGCGTCGTGAAGCGCAGGACGCTCGAGGTGACAAAGCCTGGACGCCGAAAAATCGTGAACGTCAGGTCTCCTTTGCCCTGCGTGCTTATGCCAGCCTGGCAACCAGCGCCGACAAAGGCGCGGTGCGCGATAAATCGAAACTGGGGGGTTAATAATGGCTGACTCGCAACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCC  >  NC_000913/3955071‑3955547

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