Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	77,117	G→A	R61R (CGC→CGT)	sgrR ←	DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	77,117	0	G	A	98.0%	147.9 / ‑6.6	50	R61R (CGC→CGT)	sgrR	DNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/24);  total (26/24)
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.49e-01

TTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCAATAATTCC  >  NC_000913/76871‑77362                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ttCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGcc                                                                                                                                                                                                                                                    >  2:59370/1‑251 (MQ=255)              aaTGCGCTGCCAGGTAGCAGATTACGCCACGGACGATAGTATAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTc                                                                                                                                                                                                                                      >  1:337410/1‑251 (MQ=255)                   gcTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACg                                                                                                                                                                                                                                 >  2:196715/1‑251 (MQ=255)                     tGCCAGGTAGCAGATTCCGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGGCAGCCGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTc                                                                                                                                                                                                                               >  2:341847/1‑251 (MQ=255)                     tGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCACAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTc                                                                                                                                                                                                                               >  2:196622/1‑251 (MQ=255)                        cAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTcagc                                                                                                                                                                                                                            >  1:293557/1‑251 (MQ=255)                              gCAGATTACGCAACGGAAGATAGGAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTGGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCcg                                                                                                                                                                                                                      <  2:94464/251‑1 (MQ=255)                                     aCGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCt                                                                                                                                                                                                               <  1:319276/251‑1 (MQ=255)                                       gCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGc                                                                                                                                                                                                             >  1:291/1‑251 (MQ=255)                                               cGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGtt                                                                                                                                                                                                     <  1:64232/251‑1 (MQ=255)                                                              cgcAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGTCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCat                                                                                                                                                                                      <  1:68235/251‑1 (MQ=255)                                                               gcAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTGTAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCata                                                                                                                                                                                     <  1:96128/251‑1 (MQ=255)                                                               gcAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCat                                                                                                                                                                                      >  1:364899/1‑250 (MQ=255)                                                                 aGGATGTGCCGCCCCTGTCGGAAGCTGCGTCCCAGATGAGCAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCAATAGTTGGTCGAGACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGGATAGAGAAATGTCAGACGCGACCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCTGGCATGGGGTTGAGCAGGGTggcgatatg                                                                                                                                                                                   >  2:140402/1‑243 (MQ=255)                                                                 aGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATg                                                                                                                                                                                   >  2:100199/1‑251 (MQ=255)                                                                      gtgCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATgacgac                                                                                                                                                                              <  1:53573/251‑1 (MQ=255)                                                                          cgccCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGAcgcgg                                                                                                                                                                          >  1:386939/1‑250 (MQ=255)                                                                              ccTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGagc                                                                                                                                                                        >  1:2480/1‑249 (MQ=255)                                                                                    ggAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATCCGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAAt                                                                                                                                                                <  1:1994/251‑1 (MQ=255)                                                                                       aGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGCAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCT‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGATCAATacc                                                                                                                                                             >  2:458973/1‑249 (MQ=255)                                                                                                ccAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGCATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCCCGTCAGCCAGCCGCGATCCTGCATGGCGTTGCGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATACCGCTTCCAGt                                                                                                                                                    >  2:12092/1‑251 (MQ=255)                                                                                                ccAGATGAGAAACCAGCATGTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGt                                                                                                                                                    <  1:331804/251‑1 (MQ=255)                                                                                                  aGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTc                                                                                                                                                  >  1:103125/1‑251 (MQ=255)                                                                                                  aGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTc                                                                                                                                                  <  1:100199/251‑1 (MQ=255)                                                                                                        gaAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTg                                                                                                                                              >  1:428851/1‑249 (MQ=255)                                                                                                               gCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTtgtg                                                                                                                                     <  1:341847/251‑1 (MQ=255)                                                                                                                 aTTTGCCGCCCAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCGTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTc                                                                                                                                   <  2:428851/251‑1 (MQ=255)                                                                                                                      ccGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCg                                                                                                                              <  1:361037/251‑1 (MQ=255)                                                                                                                      ccGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCg                                                                                                                              >  1:373509/1‑251 (MQ=255)                                                                                                                                 cttttccGCCAACCAACTGCTCTAGTTGGTCGTTCCGTTCCTGCCTCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTAccc                                                                                                                      <  2:386939/242‑1 (MQ=255)                                                                                                                                   ttGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCgc                                                                                                                  >  1:145219/1‑250 (MQ=255)                                                                                                                                             ccaacTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCcac                                                                                                       <  1:59370/251‑1 (MQ=255)                                                                                                                                                   gCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCg                                                                                                  <  1:391418/250‑1 (MQ=255)                                                                                                                                                     actaGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGAt                                                                                               >  1:133971/4‑251 (MQ=255)                                                                                                                                                        aGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGa                                                                                             >  1:407462/1‑250 (MQ=255)                                                                                                                                                                         cTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGATGTTGCAGaag                                                                             >  1:61217/1‑247 (MQ=255)                                                                                                                                                                           gcccttccgagcgttCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAAGGTCAGACGCGAGCGGTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGGTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCa                                                                         <  2:171294/238‑1 (MQ=255)                                                                                                                                                                           gCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGc                                                                          >  1:323471/1‑250 (MQ=255)                                                                                                                                                                                        ttCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCa                                                                         >  1:171294/1‑238 (MQ=255)                                                                                                                                                                                         tCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGAGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGCCGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCTCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCAAAAAAgg                                                           >  1:26014/1‑251 (MQ=255)                                                                                                                                                                                         tCCGCCCGCGGTTGCTGAAGCCCCAGCCCGGTATAGAGGAATGGCCGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCCCGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATaaa                                                              <  2:358411/248‑1 (MQ=255)                                                                                                                                                                                              ccgcTGTTGCGGACGCGCAAGCTCTGTATATAGGAATGGTAGACGCTATCTTTCACCACGCCCTGTCT‑CCGCCTCCCTCGTCAGCCAGCCGCGATCCTGCATGGTTTTGATCGGGGTGCGCATATGACGACGCTAGCAGCGCAATAACCCTGCCAGTTCGTGGAGCGTTGTGTCCTGCGATTGACCCCAGTAGCATTGCGGCAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGTACAAAAGGGGa                                                        <  2:123501/249‑1 (MQ=255)                                                                                                                                                                                                         gAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GCCTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCaaa                                           <  2:323409/251‑1 (MQ=255)                                                                                                                                                                                                              gcCAGCACGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTGCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACa                                      <  2:364899/251‑1 (MQ=255)                                                                                                                                                                                                                cAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGc                                    <  1:196715/251‑1 (MQ=255)                                                                                                                                                                                                                                aaTGTCAGACGCGAGCGTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCta                    <  1:440658/251‑1 (MQ=255)                                                                                                                                                                                                                                       gacccgAGCGTTTACCGCGCCC‑GCCTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGCTGAACTGTTGTTGCAGACGAGCAGATGGCATAAACGGGGAAATCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTa               >  2:184150/5‑249 (MQ=255)                                                                                                                                                                                                                                             ggCGTTTACCACGCAC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCaata       <  2:407462/250‑1 (MQ=255)                                                                                                                                                                                                                                                gTTTACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCAATAAtt    >  1:59448/1‑251 (MQ=255)                                                                                                                                                                                                                                                   tACCACGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCAATAATTcc  <  1:231564/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       TTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCC‑GACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCAATAATTCC  >  NC_000913/76871‑77362

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