Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,038	A→G	A216A (GCA→GCG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,038	0	A	G	98.4%	184.4 / ‑4.3	63	A216A (GCA→GCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (27/34);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

GTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAG  >  NC_000913/22797‑23282                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     gTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGAGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATc                                                                                                                                                                                                                                               >  2:448187/1‑249 (MQ=255)                ttATCCGGCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGc                                                                                                                                                                                                                               <  2:356107/250‑1 (MQ=255)                  aTCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCCTGGACTTCAAAACTGAAGCCACCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACCCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAACGGACGTTGAGTCTTCCGACAAAACTTCTCCGTCCATCTACGTTGCTTTT‑CAGTCGGTCGATCAGGATGCACTGAAAGCaa                                                                                                                                                                                                                             >  2:243926/1‑250 (MQ=255)                  aTCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCaa                                                                                                                                                                                                                             >  1:425430/1‑250 (MQ=255)                          gggCGTGCTGGGCGACTGGTTGCACTCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTg                                                                                                                                                                                                                       <  2:292874/248‑1 (MQ=255)                              gTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGACGCCAACATCATCCGCGCGCTGGGCAAACTCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGAGCGTTGACTGCCGTTCTGCGCTGGCGGCAACGGACGTTGAGTCTTACGACAACACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACCGAAAGCAACAtttacctta                                                                                                                                                                                                                  >  2:320997/1‑243 (MQ=255)                              gTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTa                                                                                                                                                                                                                  >  1:229796/1‑249 (MQ=255)                                  tGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCGAAACTGAAGCCACCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAACACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCCGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAAc                                                                                                                                                                                                            >  2:255105/1‑251 (MQ=255)                                  tGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAAc                                                                                                                                                                                                            >  1:120092/1‑251 (MQ=255)                                   gggcgtatGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGGGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCaa                                                                                                                                                                                                             <  2:348562/242‑1 (MQ=255)                                   gggCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAAc                                                                                                                                                                                                            <  1:303790/250‑1 (MQ=255)                                                 aCCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGATGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTTGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAATCAAAATTTGCCGTAAGCAACGTTAACGGCCCAAtc                                                                                                                                                                                             <  1:369689/251‑1 (MQ=255)                                                   ccGTACCTGACCATGGACTTCAAAACGGAAGCCAACTTCAGCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCGGCACAAAGGCGCGAAGCCGGTACACGGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTGATCGGCCCAAtcgc                                                                                                                                                                                           <  2:231977/251‑3 (MQ=255)                                                     gTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGt                                                                                                                                                                                                          >  1:55836/1‑234 (MQ=255)                                                     gTACCTGACCATGGACTTCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGt                                                                                                                                                                                                          <  2:55836/234‑1 (MQ=255)                                                              cATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGAAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCt                                                                                                                                                                                >  1:281383/1‑251 (MQ=255)                                                                 ggaattCAAAACTGAAGCCAACATCAACCGCGCGCTGGCCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGACGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGAGCTGGa                                                                                                                                                                             <  2:269656/247‑1 (MQ=255)                                                                  gACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGTCAAAAGCATCGGCAACGGTCACCTGCCCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATGTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGac                                                                                                                                                                            <  2:500961/251‑1 (MQ=255)                                                                     ttCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                                                                                                                         >  2:217186/1‑251 (MQ=255)                                                                      tCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACGTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                                                                                                                         <  2:393554/250‑1 (MQ=255)                                                                             tGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGt                                                                                                                                                                 >  1:8002/1‑251 (MQ=255)                                                                              gAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTg                                                                                                                                                                <  1:65444/251‑1 (MQ=255)                                                                                aGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                               >  2:251561/1‑250 (MQ=255)                                                                                             cgcgcgCTGGGCAAACTCACCGGTAACGGTCACCTGTGCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAACCTGCTCCGTCCATCGACGTTG‑TTTCTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCTCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcctgcgaa                                                                                                                                                 <  1:20772/251‑4 (MQ=255)                                                                                                 cgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTAAGCCGTGGACTCTGCCTGCGAACcgc                                                                                                                                             >  1:496499/1‑251 (MQ=255)                                                                                                     gggCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCC‑AATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAAt                                                                                                                                        <  2:237385/251‑1 (MQ=255)                                                                                                      ggCAAAATCATCGGCATCTGGCCCCTGCACAAAGGCGCGATGCCGGTACACTGGTGCGGTGACTGCCGTTCTGCGCTGGCGGAAGCTGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGGCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACTTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAATCGCGCAAt                                                                                                                                        <  2:400878/251‑1 (MQ=255)                                                                                                              catcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTggggat                                                                                                                                                                                  >  2:29508/1‑197 (MQ=255)                                                                                                              catcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtgat                                                                                                                                                                                  <  1:29508/201‑4 (MQ=255)                                                                                                                              cTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACt                                                                                                                 >  1:482778/1‑250 (MQ=255)                                                                                                                                tcacaAAGGCGCGAAGCCGGTACACTGGTGCGTTGTCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTATCCGTCCATCGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTATCGGCCCAAGCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCGTTTGCACTGGATTGCGACTAt                                                                                                               <  2:503381/249‑1 (MQ=255)                                                                                                                                    aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                          >  2:115229/1‑251 (MQ=255)                                                                                                                                    aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                          >  1:379488/1‑251 (MQ=255)                                                                                                                                     aaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGCCTATGCGCTg                                                                                                         >  2:125876/1‑251 (MQ=255)                                                                                                                                                 cggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACAGAAAGCAAAATTTGCCGTCAGAAAAGTTAACGGCCCAATCTCGCTGGTGATCTGTACCACTACGCCGTGGCATCTGCCTGCAAACCGCGCAATCTCTATTCCACCTGATTTCGACAATGCTCTGGAGCAGTTCGAc                                                                                             >  2:274235/6‑251 (MQ=255)                                                                                                                                                  ggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                            <  1:265486/247‑1 (MQ=255)                                                                                                                                                  ggtacACTGGTGCGTTGACTGCCGTTCGGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                            <  2:149249/247‑1 (MQ=255)                                                                                                                                                       aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGATGGTGCAGATCGACGGTc                                                                                         >  1:368787/1‑249 (MQ=255)                                                                                                                                                                   aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                           <  2:86781/251‑1 (MQ=255)                                                                                                                                                                   aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGAGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                           >  2:306621/1‑251 (MQ=255)                                                                                                                                                                           tCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAg                                                                    <  1:251561/250‑1 (MQ=255)                                                                                                                                                                              gcgcGGGCTGAAGCTGAAGTTTAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTTAAAGCAAAATTTGCCTTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTTTATTGCATCTGATTTCGACTATGCGCTGGTGCAGATCGTCGTTCAGGTCGTGATTCTGGCGAAAGATCt                                                                <  1:360000/251‑1 (MQ=255)                                                                                                                                                                                        aagcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAACATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGc                                                      >  1:230879/1‑251 (MQ=255)                                                                                                                                                                                             gaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGATTCTGGCAGAAGATCTGTTTGaagacgcaat                                                 >  2:450379/1‑243 (MQ=255)                                                                                                                                                                                               agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                               >  1:117650/1‑251 (MQ=255)                                                                                                                                                                                                    attactACGACAAATCTTTTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGTATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTTGTGATCTGGACCACTACGCCGTGGACTCTGCCTTCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                          <  1:54215/246‑1 (MQ=255)                                                                                                                                                                                                    aggaTTACGACAAAACTTCTCCGTCCAACGTCGTTGCTGTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCCCTGGTGAGCTGGACCACTACGCCGTGGACTCTGCCGGCGAACCGCGCAATCTCTATTGCACCTGATTTCGAATATGCGCTGGTGCAGATCGACGGTCAGGTCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTACTGCAGCGt                                          <  2:351610/248‑1 (MQ=255)                                                                                                                                                                                                    aGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                          <  2:425430/251‑1 (MQ=255)                                                                                                                                                                                                                aaaCTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTggcc                              >  1:228682/1‑248 (MQ=255)                                                                                                                                                                                                                        ccggccATCGACGTTTCTTTT‑CAGGCGGACGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGACCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCGGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGGATCTGCGTGACCGATTac                        <  2:225698/245‑1 (MQ=255)                                                                                                                                                                                                                         cGTCCATCGACGTTGCTTTT‑CAGGCGGTCGGTCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctg                                                                                                                                                         >  1:277680/1‑119 (MQ=255)                                                                                                                                                                                                                         cGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc                     <  1:78738/251‑1 (MQ=255)                                                                                                                                                                                                                         cGTCCAGCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctg                                                                                                                                                         <  2:277680/119‑1 (MQ=255)                                                                                                                                                                                                                          gtccctctaCGTTGCTTTT‑CAGGCGGTCGACCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCTAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCa                    <  2:281383/243‑1 (MQ=255)                                                                                                                                                                                                                            ccattggcGTTGCTTTT‑CAGGCGGTCGATTAGGATGCACAGAAAGCAAAATTTTCCGTAAGCAACGTTAACGGCCCAATCTATCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCTCCGGAGTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCAt                   <  2:278228/243‑1 (MQ=255)                                                                                                                                                                                                                             cATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTTACCGATTACAcccttt                 >  2:6420/1‑247 (MQ=255)                                                                                                                                                                                                                              agcGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                <  1:267733/249‑1 (MQ=255)                                                                                                                                                                                                                                cGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                <  1:356569/249‑1 (MQ=255)                                                                                                                                                                                                                                cGACGTTGCGTTT‑CAGTCGGTCGATCAGGATGCACTGAAAGCAAAAGTTTCCGTAAGAAACGTTAACGGCCCAATCTCGCGGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATtctca              <  2:460447/251‑2 (MQ=255)                                                                                                                                                                                                                                  aCGTTGCTTGT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCGGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagc            <  1:74963/251‑4 (MQ=255)                                                                                                                                                                                                                                       gCTTTT‑CAGGCGGTCGATCAGGATGCCCTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACgg        <  2:32924/250‑1 (MQ=255)                                                                                                                                                                                                                                            ttcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTaaa    >  1:445739/3‑250 (MQ=255)                                                                                                                                                                                                                                              cAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAg  <  1:243926/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     GTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAG  >  NC_000913/22797‑23282

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