Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,951,811	G→A	intergenic (+1254/‑390)	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,811	0	G	A	93.9%	205.4 / NA	66	intergenic (+1254/‑390)	ilvX/ilvM	uncharacterized protein IlvX/acetolactate synthase II subunit IlvM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/41);  minor base C (4/0);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACAT  >  NC_000913/3951568‑3952050                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  gCGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGAt                                                                                                                                                                                                                                             >  1:191923/1‑250 (MQ=255)  cGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATAc                                                                                                                                                                                                                                           <  1:213212/251‑1 (MQ=255)   gTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCCGTTGGCGCGCAAGTCGCGCGACCAAACGATAcc                                                                                                                                                                                                                                          >  1:315695/1‑251 (MQ=255)     gatgaACATGCCTGTCGGTACGACCAGCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCTGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGt                                                                                                                                                                                                                                        <  2:14752/251‑1 (MQ=255)          aCATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCt                                                                                                                                                                                                                                   <  1:296083/251‑1 (MQ=255)          aCATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATAACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGGTGGCGCGCAAGTCGCGCGACCCAAccgcaccgtcgtct                                                                                                                                                                                                                                   >  2:20118/1‑238 (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)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  GCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCT‑GTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACG‑CAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACAT  >  NC_000913/3951568‑3952050

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: ‑