Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,161,248	G→T	G61V (GGC→GTC)	fabR →	DNA‑binding transcriptional repressor FabR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,248	0	G	T	95.5%	137.1 / ‑4.8	45	G61V (GGC→GTC)	fabR	DNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (27/15);  total (29/16)
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.81e-01

AAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCC  >  NC_000913/4161003‑4161492                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aaaGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                                                                                                                                                                                 >  2:124075/1‑251 (MQ=255)  aaGAGTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGGTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGc                                                                                                                                                                                                                                                >  2:193531/1‑251 (MQ=255)        tGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAc                                                                                                                                                                                                                                          >  1:26965/1‑251 (MQ=255)               aaTATGTTTTATTGCGTTACTGTTCTTTCAAAATAGTGGAGCAATCGAGTATGTTCATTCTCTGATATATTGCCAGCAGTACTTTGGGAAAGGTTACAGACATCGTGTGGGGCGTAAGAGCGAAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                                                                                                                                                                                   <  1:354827/251‑1 (MQ=255)                      tttATTGCTTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAATTGGCGCTTGAAGCGGGCATTGCTCCCCCCTCTTTTTATCGGc                                                                                                                                                                                                                            >  2:281347/1‑251 (MQ=255)                            gCGTTACCGTTCATCCACAATACTGGAGAAATCCAGTATGTTCATTCTCTGGTATAGTCCCAGCAGTACTTTTGGAAAGGATTCAGACATGGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                                                                                                                                                                                      <  2:4248/251‑1 (MQ=255)                                    gTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                                                                                                                                                                                                               >  1:196404/1‑250 (MQ=255)                                            acaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                                                                                                                                                                      >  2:107988/1‑251 (MQ=255)                                                  ctgtagccctcccgTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAACCGCATTTAGCCAATTAAGTGCCGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGACATCTCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt                                                                                                                                                                                                >  1:30305/14‑251 (MQ=255)                                                    ggAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTg                                                                                                                                                                                               >  1:76847/1‑250 (MQ=255)                                                           tCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                                                                                                                                                                                       >  1:317784/1‑251 (MQ=255)                                                                     tctttCTCTGGTATAGTGCCAGCAGTCTTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGagag                                                                                                                                                                              <  1:251060/247‑1 (MQ=255)                                                                     tCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                                                                                                                                                                                >  2:251060/1‑248 (MQ=255)                                                                           tctGGTCTACTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCg                                                                                                                                                                            >  2:53446/1‑246 (MQ=255)                                                                           tctGGTATAGTGCCAGCAGTACTTTGGGAAAGGATTCAGACATCGTGAGGGGCTAAAGAGCCCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                                                                                                                                                                        <  1:371071/250‑1 (MQ=255)                                                                             tGGTATAGTGCCAGCATTACTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTaa                                                                                                                                                                      <  1:81723/250‑1 (MQ=255)                                                                                tataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCt                                                                                                                                                                  >  2:84889/1‑251 (MQ=255)                                                                                           cagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTcc                                                                                                                                                       >  2:90297/1‑250 (MQ=255)                                                                                                                 gACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTAt                                                                                                                                 >  1:264718/1‑251 (MQ=255)                                                                                                                    aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTACTGCTACGCCAACTCATGCGCCAGGCGCGTCCGCGTATCGc                                                                                                                              >  2:288847/1‑251 (MQ=255)                                                                                                                           tGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGc                                                                                                                       >  1:28312/1‑251 (MQ=255)                                                                                                                            gggCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCg                                                                                                                      >  1:208067/1‑251 (MQ=255)                                                                                                                            gggCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCg                                                                                                                      >  1:208052/1‑251 (MQ=255)                                                                                                                                              aaaaGAAAAAAAGCGCCGTTCGCTGGTGGAAGCCGCATTTACCCAATTAAGTGCTGAACGCCGCTTCGCAAGCGTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGATCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGCCGAACTGGGTCTGACCATGGTTGATGAAAGCGGTTTAATGCTACGCCAACTAATGCACCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGGGATCCGCAccc                                                                                                    >  1:240743/1‑250 (MQ=255)                                                                                                                                                        ccccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCcac                                                                                          <  1:107988/250‑1 (MQ=255)                                                                                                                                                                  cGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGGGATCCGCACCTCGGTCTCCACATTTATgg                                                                                  >  1:70876/1‑249 (MQ=255)                                                                                                                                                                   gCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGt                                                                               <  1:359990/251‑1 (MQ=255)                                                                                                                                                                          gAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCTGGAAGCGGTCATTGCTCCCACCTCTGTTTATCGGCATTTCCGCGACGGAGACGAAATGGGGCTGACCATGGTTGATGAGGGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCCCACCTCGGTCCCCACAATTATGGAGttaagcgg                                                                        >  1:206347/1‑245 (MQ=255)                                                                                                                                                                           aaGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGGTTAATGCTACGCCAGCTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGGTCATCGGt                                                                       <  1:231425/251‑1 (MQ=255)                                                                                                                                                                               cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCATGCGTTTGCGTGAAGTGGAGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGGGACGTAGGCGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGGGATCCGCACCTCGGTCCCCACATTTATGGAGTTCATCTGtaat                                                                    >  1:189225/1‑250 (MQ=255)                                                                                                                                                                                        ccAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTGTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACg                                                          <  2:98196/251‑1 (MQ=255)                                                                                                                                                                                             tAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTc                                                                                                                                                        >  2:364651/1‑152 (MQ=255)                                                                                                                                                                                             tAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTc                                                                                                                                                        <  1:364651/152‑1 (MQ=255)                                                                                                                                                                                              aaGTGCTGAACGCAGCTTCGCCAGGCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                                     <  2:171582/250‑1 (MQ=255)                                                                                                                                                                                                    tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGTACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGAt                                                                                                            >  1:91732/1‑189 (MQ=255)                                                                                                                                                                                                    tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGTACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGAt                                                                                                            <  2:91732/189‑1 (MQ=255)                                                                                                                                                                                                      aaCGCAGCTTCGCCAGCCTGCGTTTGCGTGAAGTGACGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttattt                                            >  2:38720/1‑250 (MQ=255)                                                                                                                                                                                                       aCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                           <  2:287576/251‑1 (MQ=255)                                                                                                                                                                                                        cGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCCCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTattatt                                          >  2:369342/1‑251 (MQ=255)                                                                                                                                                                                                                        tGAGTTGGCGGGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCGTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAGGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAAAGCCTTCCGGTTATTATTGCGGGAACGCTCCg                            <  1:358433/249‑1 (MQ=255)                                                                                                                                                                                                                               gCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCg                   >  1:355672/1‑251 (MQ=255)                                                                                                                                                                                                                                 gTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCt                 >  1:13889/1‑251 (MQ=255)                                                                                                                                                                                                                                               gAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCAGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGGGTGATCCGCACCTCGGGCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTgccgc   >  1:233980/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaGCGGTCATTTCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTgccgcc  <  1:298664/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaGCGGTCATTGCTCCCACCTATTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGTCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATTGAGTTCATCGGTAATAATCCTAACGCCGTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTgccgcc  <  1:298674/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     AAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCC  >  NC_000913/4161003‑4161492

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