Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	22,468	C→G	P26P (CCC→CCG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	22,468	0	C	G	95.1%	182.1 / 3.1	61	P26P (CCC→CCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (31/27);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-02

AGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTG  >  NC_000913/22225‑22702                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           agtagaGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGc                                                                                                                                                                                                                                     >  2:435486/1‑251 (MQ=255)   tagtGCTGCGTAAAAAAATACGCCATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGAATATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTg                                                                                                                                                                                                                                   >  1:402199/1‑251 (MQ=255)   tagtGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCACGCGCGAACCGGGAATGCTg                                                                                                                                                                                                                                   >  1:297583/1‑251 (MQ=255)              aaaaaaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGa                                                                                                                                                                                                                         <  2:181992/250‑1 (MQ=255)                 aaaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtg                                                                                                                                                                                                                      >  1:342798/1‑250 (MQ=255)                  aaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgat                                                                                                                                                                                                                    <  2:676968/251‑1 (MQ=255)                         cAATGCGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGACTTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCACGCGCGAACCGGGAATGCTGGCGCGTTGGACTGCTGTTGATc                                                                                                                                                                                                             >  2:307024/1‑251 (MQ=255)                                             cGCTGGACGAACTGAACGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGATTAAGCCTGTAATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGATAGACTATAAATCAACCCTGCATTTGCCGGAAACAGAGTTCCAGCTGCTATGCGATCACGCCAAGCGCGAACCGGGAATGCTGGCACGTTGGCATGATGATAATCTTTACGACATCATCCATGcc                                                                                                                                                                                         >  1:273948/1‑250 (MQ=255)                                               cTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGc                                                                                                                                                                                       >  1:522496/1‑251 (MQ=255)                                                tGGACGAACTGAACGCGAAGATTGCGCGTGATGATTTCACCGCCCGCGAATTTTTTGGGATAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGTAACACGGCACCGAGACTCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGACACAGGGTTCCCGCTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGCTGATCTGTACGGCATCATCCGAGCggat                                                                                                                                                                                      >  2:276331/1‑249 (MQ=255)                                                 ggACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGc                                                                                                                                                                                       <  1:669245/249‑1 (MQ=255)                                                       aCTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTGTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACGCGGAACCGAGAATCGGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTGCCCGAGGCGTGGCGATCGGGCCGAGGGCGAACCGGGAATGCTGGCGCGTGGTACGGATGATGATCTGTACGGCATCATCCGTGGGGCTAAAAAAg                                                                                                                                                                               <  1:389366/251‑1 (MQ=255)                                                          gAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAACCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAACAAAGGc                                                                                                                                                                             >  1:240636/1‑250 (MQ=255)                                                          gAAAGCGCAGATTGCGCGTGATGAATTAACCGCCAGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGCTGCGTGGCGATCTCGCCAAGCGAGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGGACGGCATCATCCGTGCGGCTaaaaaggga                                                                                                                                                                             >  2:116897/1‑246 (MQ=255)                                                             aGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGTCGATCTCGCCAAGCGCGACCCGGGAATGCTGGCGCGTTGGGATGATGATGATCTGTACGGCATCATCCGGGCGGCTaacaacgccaaa                                                                                                                                                                          >  2:459193/1‑240 (MQ=255)                                                             aGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCaaaa                                                                                                                                                                         <  2:342798/251‑1 (MQ=255)                                                               cgcAGATTGCGCGTGATGAATTAACCGCCAGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGACTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGTTGATGATCTGTACGGCATCATACGAGCGGCTAAAAGAGGCaaaaa                                                                                                                                                                        >  2:347187/1‑250 (MQ=255)                                                                    aTTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCtt                                                                                                                                                                  >  1:572193/1‑251 (MQ=255)                                                                       gcgcGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAAGCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCa                                                                                                                                                                <  1:633930/250‑1 (MQ=255)                                                                          cgtaatgaATTAACCGCCCGCGAATTTTTTGGGCTACCACCACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGACTTTGCCGGACCCAGGGTTCCCGCTCCGTGCCGATCTCGCCAAGCTCGAACCGGGAATGCTGTCGCGTTGGACTGATGCTGTTCTGTACGGCAGCATCCGTGCGGCTAAAAACGGCAAATAACCCTTCATTCt                                                                                                                                                            >  2:25876/5‑251 (MQ=255)                                                                           gTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATATCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCACACAAACCTTCATTCTg                                                                                                                                                           >  1:311752/1‑251 (MQ=255)                                                                                    tAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAAACGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGCCTGATGATGATCTGTACGACATCATCCGTGCGGCTAAAAAAGGCAAAAAACCCTTCATTCTGCATGATGGc                                                                                                                                                  >  2:666613/1‑251 (MQ=255)                                                                                             gcgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCACGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTAt                                                                                                                                         >  2:120142/1‑251 (MQ=255)                                                                                              cgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCCTCATCCGTGCGGCTAAACAGGGCAAAAAAACCTTCATTCt                                                                                                                                                            >  2:544895/1‑231 (MQ=255)                                                                                              cgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCt                                                                                                                                                            <  1:544895/231‑1 (MQ=255)                                                                                                     tttGGGCTTACAAAACCGGCTTAAGCCCTTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCGGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATg                                                                                                                                  <  2:673128/250‑1 (MQ=255)                                                                                                         gactaaacaaccctgCTTAAGCCTGTTATGTAATCCAACCGAAACACGGAACCGAGAATCAGATGAGTGACTATAACTCAACCATGAATTTGCCGGAAACAGGCTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGTTGATCTGTACGGCTTCATCCGTACGGCTAAAAACGGCAAAAAAACCTTCATACTGCATGATGGCCCTCCTTATGCGACTGgaagc                                                                                                                             >  2:65221/15‑251 (MQ=255)                                                                                                          gCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAAACGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTACCTTCAGCATGATGGTCCTCCTTATGAGAATGgaagc                                                                                                                             >  2:203341/1‑250 (MQ=255)                                                                                                             aaCAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACAGAGAGTCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCCCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTACAAACGGCAGAAAAACCTTCATTCTGCATGATGACCCTCCTTATGCGACAGGCAGCAtt                                                                                                                          >  2:338140/1‑250 (MQ=255)                                                                                                              aCAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTc                                                                                                                         >  1:105271/1‑250 (MQ=255)                                                                                                              aCAAAACCGGCTTAAGCCTGTTATGTAATCAAACAGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                                                                                                                                                                                                       >  2:172475/1‑156 (MQ=255)                                                                                                              aCAAAACCGGCTTAAGCCTGTTATGTAATCAAACAGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                                                                                                                                                                                                       <  1:172475/156‑1 (MQ=255)                                                                                                                aaaaCCGGCTTAAGCCTGTTATGTAATCAAAACGAAACCCGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCat                                                                                                                       >  1:183049/1‑250 (MQ=255)                                                                                                                  aaCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATAtt                                                                                                                    <  1:228810/251‑1 (MQ=255)                                                                                                                   aCCGGCTTAAGCGTGTTATGAAATCAAACCGAAACACGGACCCGAGAATCTGCGGAGTGACTATAAATCAACCCTGAATTTGACAGAAACAGGGTTCCCGATGCGTGGCGATCTCGCAAAGCGCGAACCGGGAATGCTTGAGCATTGGACTGATGATGATCTGTACGGAATCATCCGTGCGGCTAAAAAAGGACACAAAACCTTCATTCTGCATCATGCACCACCTTAGACGAGTGGCAGAATTCaacttg                                                                                                                   >  2:73482/1‑246 (MQ=255)                                                                                                                     cGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGt                                                                                                                 <  1:276331/251‑1 (MQ=255)                                                                                                                          tAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACt                                                                                                             <  1:363361/250‑1 (MQ=255)                                                                                                                           aaGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTc                                                                                                            <  1:335030/250‑1 (MQ=255)                                                                                                                                   tATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGATTGACTATAAAGCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACa                                                                                                    <  2:294274/250‑1 (MQ=255)                                                                                                                                       tAATCAACCCGAAACACGGAACCGTGAATCTGATGAGTGTCTATAAATCAACCCTGAATTTGCCGGAATCTGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGCATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGAt                                                                                                <  2:145571/250‑1 (MQ=255)                                                                                                                                           cAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGa                                                                                           <  1:22285/251‑1 (MQ=255)                                                                                                                                            aaaCCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGaa                                                                                          <  2:102501/251‑1 (MQ=255)                                                                                                                                                cGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTCAACAAGATTCTGAAAGAc                                                                                      <  1:298835/251‑1 (MQ=255)                                                                                                                                                      aCGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAAGTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATGCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATc                                                                                <  2:315689/251‑1 (MQ=255)                                                                                                                                                               agaATCTGATGAGTGACTATAAATCACTCCTGAAGTGGCCGGAAACAGGGTTCCCGATGAGTGGCGAACTCGCCAATCGCGAACCGGGAATGCTGGCGCGATAGACTGATGATGATCTGTACTGCATCTTCCTTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGTCCCTCCTTATGCGAATGGCAGCATTCATATTTGTCACTCGATTAACAAGATTCTGAAAGACATTATCGTGAAGTcc                                                                       <  2:188533/251‑1 (MQ=255)                                                                                                                                                                       atgaGTGCCTATAACTCAACCCTGAAGTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAACGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCTGCTAAAAAAGGCAAAAAAACCTTCATTCTGCCTGATGTCCCTCCTTATCCGAATGTCAGCATTCATATTGGTCACTCAGTTACAAAGATTCTGCAAGACATTATCGTAAAGTCCAACGGGCt                                                               >  2:61258/1‑251 (MQ=255)                                                                                                                                                                         gaGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCtt                                                              <  1:92727/250‑1 (MQ=255)                                                                                                                                                                           gTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTcc                                                           >  1:249693/1‑251 (MQ=255)                                                                                                                                                                             gACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCgg                                                         <  1:256757/251‑1 (MQ=255)                                                                                                                                                                                      tCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACTGATGATGAGCTGTACGGCATCCTCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGACGCATTCATATTGGTAACTCGGTTAACAAGATTATGAAAGACATTCTCGTGAAGGCCAAAGGGATTTCCGGTTATGACTc                                                >  2:208830/1‑251 (MQ=255)                                                                                                                                                                                                    gCCGGAAACAGGGTTCCCGCTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTCCGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATg                                                                                                                                  >  2:626314/1‑155 (MQ=255)                                                                                                                                                                                                    gCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATg                                                                                                                                  <  1:626314/155‑1 (MQ=255)                                                                                                                                                                                                     ccGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAATCGGGAATGCTGGTGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTGTCCGGTTATGACTCGCCGTATGTGCctgg                                 <  1:573513/251‑1 (MQ=255)                                                                                                                                                                                                          aaCAGGGTTCCCGATGCGTGGCTATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACTGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCAGTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTggg                            <  2:240636/251‑1 (MQ=255)                                                                                                                                                                                                                   cccGATTCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGTACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACACGGTTAACAAGATTCCGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACGCGCCGTATGTGCCTGGCTGGGACTGCCAc                    <  1:535043/250‑1 (MQ=255)                                                                                                                                                                                                                     cGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGt                 >  1:634958/1‑251 (MQ=255)                                                                                                                                                                                                                        tGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGAAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGACTg              >  2:596832/1‑251 (MQ=255)                                                                                                                                                                                                                           gttgcgatcgtgCCAAGCGCGAACCGGAAATGCTGGCGCATTGGACTGATGATGATCTGTACGGCATCATCCGTGCTGCTAAAAAAGTCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTGATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATGCTTAAAGACATTATCGAGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGGATGTGCCTGGCTTGGACGGCTACGGGCTGCCg           <  2:315767/240‑1 (MQ=255)                                                                                                                                                                                                                             ggCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGAt         >  1:65609/1‑251 (MQ=255)                                                                                                                                                                                                                               cGATCGCGCCAAGCGCGAACCGGGAATGCTGGCGCGGTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATGATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCAACGGTCTGCCGATCg       <  2:530647/251‑1 (MQ=255)                                                                                                                                                                                                                                     cGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGTATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTg  <  2:192234/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           AGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTG  >  NC_000913/22225‑22702

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