Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,951,823	T→C	intergenic (+1266/‑378)	ilvX → / → ilvM	uncharacterized protein IlvX/acetolactate synthase II subunit IlvM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,951,823	0	T	C	98.5%	198.7 / NA	67	intergenic (+1266/‑378)	ilvX/ilvM	uncharacterized protein IlvX/acetolactate synthase II subunit IlvM
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/38);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTT  >  NC_000913/3951577‑3952069                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         aCATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCt                                                                                                                                                                                                                                                      <  1:296083/251‑1 (MQ=255)  cctgccTGGCGTTACGACCAGCCCGGTGACGCTATCTACGCGCCGTAGTTGTTAAAACAACTGTCAGTACGTAAACCTTCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATTTGGTCTGCGCAACACATCGCCCCCACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGTGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑g                                                                                                                                                                                                                                                    <  2:177726/247‑1 (MQ=255)         ggCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGGTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTcc                                                                                                                                                                                                                                             <  2:288172/251‑1 (MQ=255)                 gACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACgg                                                                                                                                                                                                                                     <  1:28529/251‑1 (MQ=255)                  cccATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGc                                                                                                                                                                                                                                    <  2:149572/250‑1 (MQ=255)                  aCCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCGGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGc                                                                                                                                                                                                                                    <  2:69403/251‑1 (MQ=255)                           gTGACGCTATCTACGCTCCGTTGTTGTTAAAACCACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCCCATCTCCCCCACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTGGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCatg                                                                                                                                                                                                                           <  1:470692/251‑1 (MQ=255)                                 cTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGTCT‑GTATCTCCCGTGACGGCTCTTTCctgatgca                                                                                                                                                                                                                      >  2:63307/1‑248 (MQ=255)                                           gcCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGGTGGCGCGCAAGTCGCGCGACCAAACGATAACGGCGGCT‑GTATCTACGGGGACGGCTCTTTCATGATGAATGAGGAagagg                                                                                                                                                                                                           >  2:52726/1‑250 (MQ=255)                                                   tgttAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCAc                                                                                                                                                                                                    <  1:99127/250‑1 (MQ=255)                                                    ggtAAAACTACTGTCAGATTGTAACCCTGCGGATTGCGTCGTGACGACATATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACCCTGGCCCGGAGAATTTCATCCCCGCCAGCGGCTTAGTCGCCATGGGGGTTGGGTGACCGGCGTTGGGTGGCGCGCAAGGCGCGCGACCAAACTATATCGGTGTCG‑GTATCTCCGGTGACGGCTCTTTCGTGATGAATGTGCAAGAGCTGGGCAcc                                                                                                                                                                                                   <  1:188275/248‑1 (MQ=255)                                                        aaaCAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGGGCACAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGGACCGTaa                                                                                                                                                                                               >  1:202389/1‑250 (MQ=255)                                                         aacaacTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCCCACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTcc                                                                                                                                                                                                                                             <  2:451648/203‑1 (MQ=255)                                                         aacaacTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTCGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGGAAGTCGCGCGACCAAACGACACCGGCGTCC‑GTATCTcc                                                                                                                                                                                                                                             >  1:451648/1‑203 (MQ=255)                                                             acTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGAGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCTGCGGTTGGCGCGCAAGGAGCGCGACCAAACCATACCGTCGTCT‑GTCTCTCCGGTCACGGCCCTTTCCTGATGAATGGTCTAGAGCTggggcccggataacgcc                                                                                                                                                                                         >  2:86215/1‑237 (MQ=255)                                                                gtcagATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTTGCGCGCCAGACGCGCGACCCAACCACACCCACGTCTTGT‑TATCCGGCGCCGGCTCTTTCATGATGACTGTGCAAGTGCTGGGccacgcacaaagtccgc                                                                                                                                                                                      >  2:94338/5‑235 (MQ=255)                                                                gtcagATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGCGCACCGTAAAACG‑CAAGc                                                                                                                                                                                     <  2:409858/247‑1 (MQ=255)                                                                       cGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCCAGTCGCGCGACCCAACGATACCGGCGGCT‑GTATCTCCGGGGACGGCTCTTTCAATATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑ccaacaggtacc                                                                                                                                                                              >  2:467456/1‑240 (MQ=255)                                                                             ccTGGGGATTGCTGCGTGACCCCAGAGGTGGGGCAGCGCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCTGCTTAGGCACCATTGGTTGTGGTTTACCGGCGGCGGTTGGCGCGCAATTCGCGCGACCAAACGATACCGTCGTCT‑GTAGCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCGGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGat                                                                                                                                                                        <  2:78858/251‑2 (MQ=255)                                                                                  ggATTGCGTCGTGACCACAGATGCGGGGCAGCACCAGATGCGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAAAGATACCGTCGTCT‑GTTTCCCCGGTGACAGCCCTTTCATGATGAATGTGCAACAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCCTTGTAAAt                                                                                                                                                                     >  1:175486/1‑249 (MQ=255)                                                                                      tGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCt                                                                                                                                                                <  1:114142/250‑1 (MQ=255)                                                                                       gCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTAACCGTTGAAAATCGTCtt                                                                                                                                                               <  2:401999/250‑1 (MQ=255)                                                                                         gtcgtGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACt                                                                                                                                                            >  1:195108/1‑251 (MQ=255)                                                                                         gtcgtGAACACCGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCCCAGCGCCCACTCTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTGTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCGTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAAGCGTCTTACt                                                                                                                                                            <  2:266724/251‑1 (MQ=255)                                                                                                ccACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAAc                                                                                                                                                     <  1:252089/251‑1 (MQ=255)                                                                                                    aGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAAc                                                                                                                                                 <  1:381738/251‑1 (MQ=255)                                                                                                      aggtgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCGCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACgg                                                                                                                                               <  2:202389/249‑1 (MQ=255)                                                                                                       tgtgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCACGCCAATCGCGCCACCAAAAGATACCGTCGTCT‑GTATCTCCGGTCACCGCTCTTTCATGATGAATGGGCAAGAGCTGGGCAACGTAAAACGACAA‑CAGTTATCATTGAAAATCGACTTACTCCATACCCAACggg                                                                                                                                              >  2:346111/1‑250 (MQ=255)                                                                                                          ggggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAg                                                                                                                                           <  1:63307/251‑1 (MQ=255)                                                                                                           gggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAgg                                                                                                                                          <  1:148946/251‑1 (MQ=255)                                                                                                            ggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACCGTTatgg                                                                                                                                         >  2:149849/1‑248 (MQ=255)                                                                                                              cagcaCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGAt                                                                                                                                       <  2:335708/251‑1 (MQ=255)                                                                                                                  aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGtt                                                                                                                                   <  1:238240/251‑1 (MQ=255)                                                                                                                  aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAAAGATAACGTCGTCT‑GTATCTCCGGGGACGGCTCTTTCATGATTAATGTGCAAGAGCTGGGCCCCGTAAAACG‑CAAGCAGTTACCGGTGAAAATCGTCCTACTCGATAAACAACGGGTAGGGAttggt                                                                                                                                   >  2:34521/1‑247 (MQ=255)                                                                                                                      gATGTGGGCTGCGCAGCACATCGCCCCCACTNGCCCGGAAAATTTCATCACCTCCAGCGGATTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTCCTCGATAGCCAACGGTTAGGGATGGTTCGAc                                                                                                                               <  2:227896/251‑1 (MQ=255)                                                                                                                       gtgtgGGCTGCGCAGCACAGCGCCCACACGCNCCCGTAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGATTTTGGTGTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGGCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCTAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTAGCGTTGAAAATCGTCTTACTCCATAACCAACGGTTATGGATGGTTCGACa                                                                                                                              <  2:234096/250‑1 (MQ=255)                                                                                                                                    agctcaTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTg                                                                                                                 <  1:241764/247‑1 (MQ=255)                                                                                                                                        caTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAATCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCCGTTAACGTTGAAAATCGTCTTACTCGAAAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttt                                                                                                             >  1:96026/1‑251 (MQ=255)                                                                                                                                         aTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGAACAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTGCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGttttt                                                                                                            <  2:459278/251‑1 (MQ=255)                                                                                                                                           cGCCCACACTCGCCCGGAAAATTTCATCACNCCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttttc                                                                                                           <  2:224874/250‑2 (MQ=255)                                                                                                                                                     cGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGa                                                                                                  <  2:195108/249‑1 (MQ=255)                                                                                                                                                        ccGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCACGTCGCGCGAGCCAACGCTAACGTCGTCT‑GTATCTCCGGGGAAGGCTCCTTAATGATGAATGTGCAAGAGCCGGGCACCGTAAAACG‑CAAGCAGCTACCGTTGAAAATCGTCTTACTCGATAACCAACGGGTAGGGGTGGGTCGACAANGGCAGAAACTGGTTTTTCAGGaaaggtct                                                                                               >  2:234853/1‑243 (MQ=255)                                                                                                                                                         cGGAACATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGCTACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTTATGTTGAATGTGCAAGAGCTTGGCACCGTAACACG‑CCAGCCGTTACCGTTGAAAATCGTATTAATCGATCAACAACGGTTCGGGATGGTTTGACCAANGCCGCAACTGTTTTTTCCAGAACGATaata                                                                                             >  2:222677/1‑247 (MQ=255)                                                                                                                                                         cGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGc                                                                                            >  1:54549/1‑251 (MQ=255)                                                                                                                                                           gAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGGCGCGCGACCAAAAGATACCGTCGTCT‑GTATTTCCCGGGACGGCTCTTTCATGTTGAATGTGCAAGAGCTGGGCACCGTAAAAAG‑CAAGCAGTTACCGTTGAAAAACGTCTTACTCGATAAACAACGGTTTGGGATGGGTTGTCAATGGCCGCAACTTTTTTTCCAGGAACGATATAGcc                                                                                           >  2:126472/1‑249 (MQ=255)                                                                                                                                                              aaTTTCAGCACCTCCAGCGGCTTAGGCACCATGGGGTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCTACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAac                                                                                       <  1:234853/251‑1 (MQ=255)                                                                                                                                                                      caCCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTa                                                                               <  1:248868/251‑1 (MQ=255)                                                                                                                                                                        ccTCCAGCGGCTTAGGCACCATGGGTTTTGGTTGAACGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCATCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTAcc                                                                             <  1:349136/251‑2 (MQ=255)                                                                                                                                                                              gCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTGGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGAGGGTTCGCCAATGGCAGCAACTGTTTTGCCAGGAACGATATAGCGAAACCACCCTTACCGATAAc                                                                       <  2:172151/251‑1 (MQ=255)                                                                                                                                                                                 gcttAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAAcccc                                                                    <  2:176540/249‑1 (MQ=255)                                                                                                                                                                                          ccATGGTTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGAGGAATGTGCAAGAGCTGGGCACCGTAAACCG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTc                                                           <  2:277833/251‑1 (MQ=255)                                                                                                                                                                                            aTGGGTTTTGGTTTCCCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCAt                                                         <  2:79566/251‑1 (MQ=255)                                                                                                                                                                                            aTGGGGTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATCCCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTGCCTCAt                                                         <  2:189470/251‑1 (MQ=255)                                                                                                                                                                                              gggTTTTGGTTTACCGGCGGCGGTTTGCGCGCAAGTCCCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACC‑CAAGCAGTTACCGGTGAAAATCGTCTTAATCGATAACCAACGGGTAGGGCTGGTTCGACAATGGCAGCAAATGTTTTTCCAGGAACGAAATAACGGAACCACCCTTACCGCTAAACCCGATTTCCTCaagt                                                       >  2:100892/1‑248 (MQ=255)                                                                                                                                                                                                gTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGtt                                                      >  1:376604/1‑250 (MQ=255)                                                                                                                                                                                                       tttACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CCAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGCTATAGCGAAAACAACCTTACNCATAACCCCGATTTCCTCCTGtttaccaccc                                              >  2:226047/1‑243 (MQ=255)                                                                                                                                                                                                                   ggTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATcc                                  <  1:40625/251‑1 (MQ=255)                                                                                                                                                                                                                       ggCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTgg                              <  1:31058/251‑1 (MQ=255)                                                                                                                                                                                                                       ggCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTgg                              >  1:150998/1‑251 (MQ=255)                                                                                                                                                                                                                          gcgcaAGTCGAGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCCGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGcc                            >  2:197122/4‑250 (MQ=255)                                                                                                                                                                                                                                gTCGCGCGACCAAACGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACAt                     <  2:115677/251‑1 (MQ=255)                                                                                                                                                                                                                                     gcgacaaaaCGATACCGTCGTCT‑GTATCTCCGGTAACGGCTCTTTCATGAAGAATGTGCAAGAGCTGGGCACCGTAACACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCCGCAACTGTTTTTCCAGGAACGCTATAGCGAAACCACCCTTACCTATAACCCCGATTTCCTCATGTTCGCCAGCGCCTTCGGCATCCCTGGCCCACACATCcccc                >  2:373882/8‑251 (MQ=255)                                                                                                                                                                                                                                         ccaaaCGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACACTGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCgta             >  1:207574/4‑248 (MQ=255)                                                                                                                                                                                                                                          caaaCGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCCGCGCCTTCGGCATCCCTGGCCAACACATCACCCgtaa            >  2:358652/3‑247 (MQ=255)                                                                                                                                                                                                                                            aaCGATACCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCgtaa            >  2:325845/1‑245 (MQ=255)                                                                                                                                                                                                                                                  aCCGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGt   >  1:137642/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ccGTCGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCCCCCGTAAAGCCCAgggt  >  2:69005/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         ACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTT  >  NC_000913/3951577‑3952069

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