Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,857,607	G→A	R172H (CGC→CAC)	mutS →	DNA mismatch repair protein MutS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,857,607	0	G	A	98.0%	163.5 / NA	50	R172H (CGC→CAC)	mutS	DNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/23);  minor base C (1/0);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACGCACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTG  >  NC_000913/2857364‑2857848                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  ccATCTGCGAACAAAGTGGCGATCCAGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCTTCAGGACAACCTGCGGGCGGCTATCTGGCAGGAAAGCAAAGTTTTCGGCTGCGCGCCGCTGGATATCAGGTCCGGTCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAGCGATGTCGGCAGAACTGCACCACACTAAt                                                                                                                                                                                                                                            <  2:205789/251‑1 (MQ=255)          aaCAAAATGGCGATACGGCGACCCGCAAAGGACCGGTTGAGCGCAAAGTTGTGCGGATCGTTACGCCAGGCACCATCAGCGATGAAACCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGCCAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAATTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCCGCACTGCAACACACTAATCCTGCggc                                                                                                                                                                                                                                    >  2:92886/1‑249 (MQ=255)            cAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAAct                                                                                                                                                                                                                                 <  1:228178/251‑1 (MQ=255)                  ggCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGt                                                                                                                                                                                                                            >  1:20379/1‑250 (MQ=255)                         cGGCGACTAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCCGCGCTGAAGCCCTGTTGCAGGAGCGTCCGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGCCGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAACCGATGGCGGCAGAACTGCCACACCCTCATCCTGCGGAACTGCTGTATGccga                                                                                                                                                                                                                     >  2:131272/1‑247 (MQ=255)                               ccAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACACCCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGAtttt                                                                                                                                                                                                              <  1:100132/251‑1 (MQ=255)                                      aGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGaa                                                                                                                                                                                                        >  1:124/1‑250 (MQ=255)                                        gTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGCTATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTTCAACACACTAATCCTTCGGAACTGCTGTATGCAGAAGATTTTGCTGaactg                                                                                                                                                                                                     >  2:48753/1‑248 (MQ=255)                                               tGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAAAGTCGTTa                                                                                                                                                                                               <  1:229523/250‑1 (MQ=255)                                                            gggCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTGGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATAGCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGAAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTc                                                                                                                                                                                  <  2:124/248‑1 (MQ=255)                                                             tgCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTcgc                                                                                                                                                                                >  1:315330/1‑251 (MQ=255)                                                                gTATCGTTACGCCAGGCACCATCCGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCACAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGc                                                                                                                                                                             >  2:315074/1‑251 (MQ=255)                                                                gTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGCCAGCAAAGGGTTCGGCTACGCGGCGCTGGATATCAGTTCCGGGCGCTTTTTCCTGAGAGAACCGGCTGCCCGCGAAACGATGGAGGCAGCACTGCAACACACTAATCCTGCGGACCTGCTGAATGCCGAAGATTTTGCTGCAATGTCGTTACTTGAAGGCCGTCGAGGc                                                                                                                                                                             >  2:234123/1‑251 (MQ=255)                                                                gTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCTTGCGGAACTGCTGTATGCAGAAGATTTTGCTGCAATGTCGTTAATTGAAGGCCGTCGCGGc                                                                                                                                                                             >  2:80669/1‑251 (MQ=255)                                                                 tATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTTCTGAAATGTCGTTAATTGAAGGCCGTCGCGGcc                                                                                                                                                                            >  2:240131/1‑251 (MQ=255)                                                                              ggCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTTAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGccc                                                                                                                                                               <  1:291888/251‑1 (MQ=255)                                                                                         cGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGGCGCTGGAGATCAGTGCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTTTGGGAGt                                                                                                                                                    <  1:269834/251‑1 (MQ=255)                                                                                                      ttGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGGATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGgcac                                                                                                                                       >  2:314561/1‑251 (MQ=255)                                                                                                              gCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGcc                                                                                                                               <  1:226039/251‑1 (MQ=255)                                                                                                                          ccTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATc                                                                                                                   >  2:297296/1‑251 (MQ=255)                                                                                                                             gctgGCGGCTATCTGGCAGGTCAGCAAAGGTTTCGGCTACGCGACGCTTGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACGGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAAGTGAAGGCCGTCGCGGCCTGCGCCGTCGTCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGc                                                                                                                <  2:55413/251‑1 (MQ=255)                                                                                                                             gctgGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCATGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGACTCTGc                                                                                                                >  1:234936/1‑251 (MQ=255)                                                                                                                                    gcTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCTGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGTCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTgg                                                                                                         <  2:45449/251‑1 (MQ=255)                                                                                                                                    gcTATCTGGCAGGACAGCAAAGGGTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAACGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACAAACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCAGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGCATCTGCAATTTgg                                                                                                         >  2:283025/1‑251 (MQ=255)                                                                                                                                           ggCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCCGTTGGATCTGCAATTTGGGACCcgc                                                                                                  >  2:235001/1‑251 (MQ=255)                                                                                                                                               ggACAGTAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATc                                                                                              >  1:90954/1‑251 (MQ=255)                                                                                                                                                gACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCt                                                                                             <  1:315074/251‑1 (MQ=255)                                                                                                                                                     cAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCg                                                                                        >  1:275085/1‑251 (MQ=255)                                                                                                                                                          gTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTTTGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAAGCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGCCCCGCGATCTGGTCGGtttt                                                                                   <  2:274071/251‑1 (MQ=255)                                                                                                                                                          gTTTCGGCTACGCGACGCTGGATAGCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAAGCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGtttt                                                                                   <  2:106699/251‑1 (MQ=255)                                                                                                                                                              cGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTGTTGGCg                                                                               <  2:275085/251‑1 (MQ=255)                                                                                                                                                                gCTACGCGGCGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTGTGGCGTc                                                                             <  1:241339/251‑1 (MQ=255)                                                                                                                                                                 cTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCg                                                                            >  2:120860/1‑251 (MQ=255)                                                                                                                                                                   aCGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTTCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCATTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCga                                                                           <  1:283025/250‑1 (MQ=255)                                                                                                                                                                   aCGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgag                                                                          <  1:8178/251‑1 (MQ=255)                                                                                                                                                                        aCGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAAcgc                                                                     <  1:48753/251‑1 (MQ=255)                                                                                                                                                                         cGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTG‑AATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAAcgcgc                                                                   >  1:119080/1‑251 (MQ=255)                                                                                                                                                                          gCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGACGGCGGCAGAACTGAAACACACTAATCCTGCGGACCTGCTGTATGCAGACGATTTTGCTGAAATGTAGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACACGCGATCTGGTCGGTTTTGGCGTCGAGAAcgcgc                                                                   >  1:258953/1‑251 (MQ=255)                                                                                                                                                                            tGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATTGCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGAc                                                                                                      <  1:32478/214‑1 (MQ=255)                                                                                                                                                                            tGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGAc                                                                                                      >  2:32478/1‑214 (MQ=255)                                                                                                                                                                                     gttacGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCCACCCACTAATCCTGCGGAACTGCTGAATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGGGGGAGTTTGAACTCGCCACGGCGCGCCAGCAGTTGAAACTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTAGAGAACGCGCAGCGCGGACtt                                                        >  2:71327/5‑251 (MQ=255)                                                                                                                                                                                     gTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACtt                                                        >  1:210801/1‑251 (MQ=255)                                                                                                                                                                                         cGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAATACTCTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTTCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTtgtg                                                    <  1:12520/251‑1 (MQ=255)                                                                                                                                                                                          gggCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGtgc                                                   >  1:20001/1‑251 (MQ=255)                                                                                                                                                                                                            gAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCa                                 <  2:315330/251‑1 (MQ=255)                                                                                                                                                                                                                         gcgAAACGATGGCGGCAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGAt                    <  1:131272/251‑1 (MQ=255)                                                                                                                                                                                                                               cgagggcggcAGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTTAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTGTGGCGTCGAGAACGCGCCGCGGGGACTTTGTGCTGCCGTTTGTCTGTTGCAGTATGCGAAAGATACCCaa              <  1:311043/247‑1 (MQ=255)                                                                                                                                                                                                                                 aTGGCGGCAGAACTGCAACACACTACTCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCGGTATGCGAAAGATACCCCACg            >  2:150881/1‑251 (MQ=255)                                                                                                                                                                                                                                         aGAACTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGActc    >  1:228892/1‑251 (MQ=255)                                                                                                                                                                                                                                           aaCTGCAACACACTAATCCTGCGGAACTGCTGTATGCAGCAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTg  >  1:199722/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  CCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCGCGAAACGATGGCGGCAGAACTGCAACGCACTAATCCTGCGGAACTGCTGTATGCAGAAGATTTTGCTGAAATGTCGTTAATTGAAGGCCGTCGCGGCCTGCGCCGTCGCCCGCTGTGGGAGTTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCAACGTACGACTCTG  >  NC_000913/2857364‑2857848

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