Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	97.6%	124.4 / NA	42	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (25/16);  total (26/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.99e-01

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACC  >  NC_000913/3485959‑3486450                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTCCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCtt                                                                                                                                                                                                                                                   <  2:23144/251‑1 (MQ=255)    gTAATACATTGATGTACGGCATGTATACAAAGGACGTCTCATTACGGTGCAGTGCAGTTGATAGCCCCTTCCCAGGTCGCGGGAAGCATATTACGGCAATCCAGAGACAGCGGCGTTATCTGGCGCGGGAGAAAGCTTATAACAGAGGATAACCCCGCATGGTGCTTGGCAACCCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAtt                                                                                                                                                                                                                                                <  1:303176/251‑1 (MQ=255)       aTACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAAAAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAATGGCAAGCTTATTCAc                                                                                                                                                                                                                                             >  2:7426/1‑251 (MQ=255)              gATGTACTGCATGTATGCAAAGGACGTCACAGTACCGTGCAGTACAGTTGATAGCCCCTCCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAACCAGACCCGACTCTCGAATGGTCCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGt                                                                                                                                                                                                                                       <  2:113492/250‑1 (MQ=255)                  tACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTATCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaa                                                                                                                                                                                                                                   >  1:283964/1‑250 (MQ=255)                          gTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAACCCCTTCCCAGGTAGAGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGCAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGAATGGTTCCTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGaa                                                                                                                                                                                                                           >  1:201519/1‑250 (MQ=255)                            aTGCAACGGTCGACACATTACCGTTCAGAACAGTTGATAGTCCCTTCCAAGGGACCGGGAAGCCTATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGATCTGGAGAAAGCTTATAACCAAGGATAACCGCGCATGGTGATAGGCAAACCGCAACCAGACCAGCCTCTCGAATGGTTCTTGTCTCATTGCCACATGCATACGTACCCATCCAAAAGCAAGCTTATTCAACAGTGTGAATAAGCGGCAACg                                                                                                                                                                                                                        >  2:348524/1‑251 (MQ=255)                                               aCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCTTTAAAg                                                                                                                                                                                                     <  1:75530/251‑1 (MQ=255)                                                  gTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGct                                                                                                                                                                                                  <  1:228394/251‑1 (MQ=255)                                                           aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                         <  1:352833/251‑1 (MQ=255)                                                             ttGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTg                                                                                                                                                                                       >  1:25903/1‑251 (MQ=255)                                                              tGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTg                                                                                                                                                                                       <  1:171141/250‑1 (MQ=255)                                                                      ccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACCTTCATACGTACCCATCCAAGAGAAAGCTTATTCACCAGCGTGAAAAAGCGGAAACGCTGTACTACAGCGTTCAAGGCTCTGTGGCAGTGCGGATCaa                                                                                                                                                                               >  2:238735/1‑250 (MQ=255)                                                                      ccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaaa                                                                                                                                                                              <  1:340585/251‑1 (MQ=255)                                                                                     ggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGACTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAACGGCTCTGTGGCAGGGCTGATCTACGACGAagcgtttca                                                                                                                                                                >  2:96876/1‑243 (MQ=255)                                                                                       aaGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaa                                                                                                                                                                >  2:108089/1‑248 (MQ=255)                                                                                                 cggcaagcaatagaCAGCGTCTTTATCTGGCACTGGAGAAAGCTAAAAACAGAGGATAACCGCGAATGGTGCTTGGAAAACCGCAAACAGCCCCGACTCTCGAAGGGTTCTTGTCTCATTGCCACATGCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGACTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATccgc                                                                                                                                                   <  2:423121/240‑3 (MQ=255)                                                                                                   gCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCtctc                                                                                                                                                 <  1:200041/251‑1 (MQ=255)                                                                                                       tCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTAt                                                                                                                                             >  1:47155/1‑251 (MQ=255)                                                                                                        ccAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATc                                                                                                                                            <  1:31424/251‑1 (MQ=255)                                                                                                        ccAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATc                                                                                                                                            >  2:110991/1‑251 (MQ=255)                                                                                                              aCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                                                                                                      >  1:438939/1‑251 (MQ=255)                                                                                                                 gcggcgATATCTGGCTCTGGAGACAGCTTATAACAGAGGATAACCGCGCATGGCGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAgg                                                                                                                                   >  2:186927/1‑251 (MQ=255)                                                                                                                    gcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGCCGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTg                                                                                                                                >  2:189584/1‑251 (MQ=255)                                                                                                                              gCTCTGGAGAAAGCTTATAACAGCGGATAAAAGTGCACGGTGCTTGGCACAACGCAACCCGACCCGACTCTCGAATGGATCTTGTCTCATTGCCACATTAATAAGTACCCATCCAAGAGCAAGCTTCTTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGCGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGCTTTTATTgg                                                                                                                      >  1:181241/1‑251 (MQ=255)                                                                                                                                                                 cATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTa                                                                                   <  1:352058/251‑1 (MQ=255)                                                                                                                                                                      tGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCCGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAgcgcc                                                                              >  2:428335/1‑250 (MQ=255)                                                                                                                                                                               acccGCAAAAAGACCCGACTCCTGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCATGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGTCAGGGCTGATCAAAGACGTAGAGGGTAAAGAAATTATCCTCTCCTATCTGAATCATGGTGATTTTATGGTCCAACTGGGCCTTTTTTAAGAGTGTCAGGAACGTAGGGCATGGGTACGt                                                                     <  2:444953/249‑1 (MQ=255)                                                                                                                                                                               aaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGCAACGCTGTACTACATCGTTAAAGGCTCTGTGTCAGTGCTGAGCAAAGACGAAGAGGGTAAAGAAATTATCCTCTTCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCGGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGt                                                                     <  2:241586/251‑1 (MQ=255)                                                                                                                                                                                       aCAGCCCCGCCTCTCGAATGGTTCTTGGCTAATCGCCACATTAATAAGTACCCATCCCAGAGCACGCTTCTTAACCAGGGTGAAAAAGCTAAAACGCTGTACTACATCGTTAACGACTCTGCGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGCAGAGGGCCTGGAACGTAGCGCATGGGTACATGCGaacac                                                             >  2:81286/1‑248 (MQ=255)                                                                                                                                                                                             ccGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGAGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtg                                                       >  1:323740/1‑251 (MQ=255)                                                                                                                                                                                                   ctcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAgg                                                 >  1:294390/1‑250 (MQ=255)                                                                                                                                                                                                        atggatattGTCGCATTGCCACATTCATAAGTACCCATCCAAGGGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGCACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCACAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGCTTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCCAGGGTACGTGCGAAAACCGCCTGTGAAGTggatg                                            >  2:328463/8‑248 (MQ=255)                                                                                                                                                                                                            ttCTTGTCTCATTGCCACATTCATAAGTACCCCTCCAAGAGCAAGCTTATTCACCAGGGTGAAAACGCGGAACCGCTGTACCACATCGTTAAAGGACCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTCACGAAATGATCCTCTCCTATCTGCATAAGGGTGATTTTACTGGCGAACTGGGCCTGTTTGAAGAGGCACAGGACCGTAGCGCATGGGTACGTGCGAAAACCGACTGTGAAGTGGCTGAAAt                                        >  1:27778/1‑251 (MQ=255)                                                                                                                                                                                                                            acaTTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGAGTGAAAAAGAGGAACAGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGAGGCTGAAATTTCGTACAAAAACttt                        >  1:379083/1‑251 (MQ=255)                                                                                                                                                                                                                                 cATAAGTACCCATCCTAGCGCAAGCTTATTCACCAGGGTGAACAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGATGATCAAAGACGAAGAGGGTAAAGAAATGATCCTATCCTATCTGAATCAGGGGGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGAGAAGTGGCGGACATTTCGTACAAAAAATTTcacc                    >  2:100018/1‑247 (MQ=255)                                                                                                                                                                                                                                  aTAAGTCCCCACCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCaa                  <  1:291783/251‑1 (MQ=255)                                                                                                                                                                                                                                   tAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTa                                                                                                                          >  2:266215/1‑146 (MQ=255)                                                                                                                                                                                                                                   tAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTa                                                                                                                          <  1:266215/146‑1 (MQ=255)                                                                                                                                                                                                                                   tAAGTACCCATCCAAGAGCAAGATTATTCACCAGGGTGAAACAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGCGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGAGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAt                 >  2:198982/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaGTACACATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCGGAAATTTCGTACAAAAAATTTCGCCAAtt                >  2:425590/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gagCAAGCTTATTCACCAGGGGGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCTCCAATTGATTCAGGTAAAcc  >  2:35920/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACC  >  NC_000913/3485959‑3486450

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