Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,787,100	C→T	7.3%	G5S (GGC→AGC)	ppsA ←	phosphoenolpyruvate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,787,100	0	C	T	7.3%	156.7 / 4.4	55	G5S (GGC→AGC)	ppsA	phosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (2/2);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

ATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTCA  >  NC_000913/1786829‑1787355                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                aTCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCAtt                                                                                                                                                                                                                                                              >  1:63305/1‑275 (MQ=255)        ggCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGAc                                                                                                                                                                                                                                                       <  2:63304/275‑1 (MQ=255)          cGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACa                                                                                                                                                                                                                                                      <  1:848855/274‑1 (MQ=255)            gatgtgctccgccGCTTTCTCAGGCTGAGTAACATCGTCAAGATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTACAGAAACTGGGTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAACAACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGCGCCATTGTTGGACATCg                                                                                                                                                                                                                                                   <  2:1286517/264‑1 (MQ=255)              tttGCGGGCCCGCTTTCGCAAGCGTGGTAAAATCGTCAATGTCCGTTTTATCCAGCAGTTCATAAATGCGCTGGGTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGATGCATTGTTGACCCTAACCCTGGCTACATCATTCATGCCGAGTTGGTTATACCAAAGTACCAGCGGTGACGAGCCATTGTTGGACATCGaa                                                                                                                                                                                                                                                 <  2:422517/275‑1 (MQ=255)              tttGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAAGATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGaa                                                                                                                                                                                                                                                 <  1:375493/275‑1 (MQ=255)                                                 tCAATATTCGTTTTATCCAGCAGTTAATACAAGCGCTGTTTTACGCCACATTGGTACAGAAACTGGTTAAACGCGACGGCGGTTGTGGCGAAACCATTCGGAACGGAAACTCCCATTCCGGAAAGATTAGCAATCATGTCACACAGGGAGGCATTTTTGCCCTCAACCCTGTCTACTACATTCATGCCGAGTTGGTTATACCAACACACCAGCGGTGCCGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATCAATGAACGGTTTGCGaaaa                                                                                                                                                                                                              >  2:235851/1‑274 (MQ=255)                                                         cggtttATCCAGCAGGTCATAAATGGGCGGGTTTACGCCGCTTGGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCAGTCGGAACTGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  2:1499616/272‑1 (MQ=255)                                                         cGTTTTGTCCAGCAGGTCATAAATGCGCTGGGTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  2:559028/275‑1 (MQ=255)                                                         cGTTTTATCCAGCAGTTAATAAATGCGCTGGTTTACGCCGCTTTGGTCCATAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGATGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  1:1497189/275‑1 (MQ=255)                                                         cGTTGTATCCAGCAGTTAATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGTAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATTCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  1:1564283/275‑1 (MQ=255)                                                                     cagTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACTTCGAACAATCCTTTTGTGCTAAATGAACGGTTTGAGAAACACTTTTCTGCGCATTTATCt                                                                                                                                                                                          >  1:265310/1‑275 (MQ=255)                                                                        ttCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAACGCACCAGCGGTGACGAGCCATTGTTGGACAt                                                                                                                                                                                                                                                     >  2:1136306/1‑213 (MQ=255)                                                                            tAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGaa                                                                                                                                                                                                                                                 <  1:1136307/213‑1 (MQ=255)                                                                                  cgcTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAAAGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATAATCTCACCCAGGGAGGCATTTTTGCCCCCAAACCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACCCATTTTTGGTCATTTATCTTTGCATATTTnan                                                                                                                                                                             >  2:1173173/1‑274 (MQ=255)                                                                                       gTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAt                                                                                                                                                                        >  1:701845/1‑275 (MQ=255)                                                                                                 cTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACAt                                                                                                                                                                                                          >  1:607104/1‑231 (MQ=255)                                                                                                  tttGGCCCAGAAGCTGGATAAACAAGTGGGCGGTTATGGCAAAACCATTCGGAAGGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTAGATCATTCATGCCGAGTTGGTTATACAAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTGCTGCGCTTTTATCTTTGCATATTTAACCGGATGAAAAAAACg                                                                                                                                                              <  2:265309/274‑1 (MQ=255)                                                                                                     ggTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  2:607103/231‑1 (MQ=255)                                                                                                            aaaCTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCAATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGTTGAAAAAAACGGTAAAACGTTc                                                                                                                                                   >  2:750800/1‑275 (MQ=255)                                                                                                                  gTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGc                                                                                                                                               >  2:646828/1‑273 (MQ=255)                                                                                                                            tCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTATGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATAttt                                                                                                                                     >  1:279536/1‑273 (MQ=255)                                                                                                                            tCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACAAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTATTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATGTTTACAGCAAATTTAttttt                                                                                                                                   >  2:482115/1‑275 (MQ=255)                                                                                                                               gcggTTGTGGCGAAACCATTCGGAACTGAAACACCCATTCCGGAAAGCTTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCTAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATTAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTATCCGGATGAAACAAACGGTGATTCGTTCAAGCAAATATATTTTTTTa                                                                                                                                <  1:235852/275‑1 (MQ=255)                                                                                                                               gcggTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTa                                                                                                                                <  1:646829/275‑1 (MQ=255)                                                                                                                               gcagTTTTGGCCAATCCATTCGGAACGGAAACACCCATGCCGGAAAGATTATTAATCAGTTCACCCATGGGGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCGGCGGTGACGGGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACAGTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTa                                                                                                                                <  2:175627/272‑1 (MQ=255)                                                                                                                                         cGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGa                                                                                                                      <  1:1101775/275‑1 (MQ=255)                                                                                                                                           aaaCCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACt                                                                                                                    <  2:701844/275‑1 (MQ=255)                                                                                                                                                                    ttCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGc                                                                                            >  2:466894/1‑274 (MQ=255)                                                                                                                                                                     agcggacggcttAGGAATCATTTGCGCCAGGGAGGCATTTTTGACCCCAACCCTGGCTACAATATCCATGCCGCGTTGCTGATGCCCAAGCAACAGTACTGACGGGCCCTGGTTGGTGATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATATTTGCATATGTAACCGGTTGAAAAAAACGGTGAATCGTTCAAGCAAATGTATTTTTTTACTTTTTAAGCCTGATCGCAGCGTTGCGCAAATCTGct                                                                                           <  2:380422/264‑1 (MQ=255)                                                                                                                                                                                 aGTAATCATTTCACCCAGTGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAAtt                                                                              <  1:884454/275‑1 (MQ=255)                                                                                                                                                                                        aTTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATaaa                                                                                                                                                                                                                               >  1:586345/1‑123 (MQ=255)                                                                                                                                                                                        aTTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCACCGGTGATGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCACGCAAATATATTTTTTTACTTTTTAAGACTGATGCCAGCGTTGCGCAAATc                                                                                               >  1:380423/1‑251 (MQ=255)                                                                                                                                                                                            cACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGaa                                                                                                                                                                                                                           <  2:586344/123‑1 (MQ=255)                                                                                                                                                                                               ccAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGTACAATCCTTTTGTGATAAATGAACTGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGAATGAAAAAAACGGTGAAGCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  1:876956/275‑1 (MQ=255)                                                                                                                                                                                               ccAGGGAGGAATTTTTGCCACCAACCCTGTCTAAATCATTCAGGCCGAGTTGGGTATACCAAAGCAGCAGCCGTGACGAGCCATTGTTGGACATCTAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCGGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  2:279535/275‑1 (MQ=255)                                                                                                                                                                                                 aGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCg                                                              >  2:1571093/1‑275 (MQ=255)                                                                                                                                                                                                                     aaCCCTGTTTACATAAGTCATGCCGAATTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATaaaa                                          <  2:352844/275‑1 (MQ=255)                                                                                                                                                                                                                            tCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATaaa                                                                    >  1:453265/1‑242 (MQ=255)                                                                                                                                                                                                                            tCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACCTTTCTGCGCCTTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATTTGCTCTCTTCGACAATTTCCACATaaa                                                                    >  1:1226203/1‑242 (MQ=255)                                                                                                                                                                                                                               aCATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCAt                                              >  1:352845/1‑261 (MQ=255)                                                                                                                                                                                                                                catcatTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  2:453264/242‑1 (MQ=255)                                                                                                                                                                                                                                   cattcatGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTTCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGttt                            <  1:750801/275‑1 (MQ=255)                                                                                                                                                                                                                                        atGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATa                       <  1:482116/275‑1 (MQ=255)                                                                                                                                                                                                                                          gCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGAAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATTTGCTCTCTTCGACAATTTCCACATACACCATCGGTATAATCGAAACCATAAATATGAAATGCTGtttt                           >  2:628927/1‑269 (MQ=255)                                                                                                                                                                                                                                           ccGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATaaaa                    >  2:296448/1‑275 (MQ=255)                                                                                                                                                                                                                                              aGTTGGTTAAACCAAGTCACCAGCGGTGACTAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGAATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATa                       <  1:628928/269‑1 (MQ=255)                                                                                                                                                                                                                                                ttGGTTATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATccc                                                                                                              >  2:810147/1‑180 (MQ=255)                                                                                                                                                                                                                                                   gTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAAGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCAAGCGTTGCGCaaa                                                                                                 >  1:1087982/1‑190 (MQ=255)                                                                                                                                                                                                                                                    ttATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCg                                                                                                          <  1:810148/180‑1 (MQ=255)                                                                                                                                                                                                                                                       taCCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTg                                                                                             <  2:1087981/190‑1 (MQ=255)                                                                                                                                                                                                                                                             aGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAgg      >  2:46796/1‑271 (MQ=255)                                                                                                                                                                                                                                                             aGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTc   >  1:858670/1‑274 (MQ=255)                                                                                                                                                                                                                                                                 ccAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAGGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTCa  <  1:46796/271‑1 (MQ=255)                                                                                                                                                                                                                                                                             aGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAg       >  2:1470441/1‑254 (MQ=255)                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                ATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTCA  >  NC_000913/1786829‑1787355

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