Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,145	G→T	intergenic (+113/‑489)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,145	0	G	T	95.6%	127.9 / ‑4.1	45	intergenic (+113/‑489)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (14/29);  total (15/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 = 9.82e-01

CACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCT  >  NC_000913/88904‑89390                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      cACCGCGATGTCGCAGAGCTGGTGCTGGAGATTGTCCTGGTAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAc                                                                                                                                                                                                                                              <  2:343393/251‑1 (MQ=255)    cgcgATGTCGCAGGGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAcccg                                                                                                                                                                                                                                           <  2:81791/251‑1 (MQ=255)       gATGTCGCAGAGTGTGGGCTTGAGATTGTCCTGGCAAGCATGGACGACCCGCGTAAGCCAAAACCTGGTTTAACTCGCATTAAACGTAAGCTCGATCGCGGCGGCGTGCTCAGCCGGAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATGTCCCTATGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTAGTGATAAGTTCTTTTAAAcccgccc                                                                                                                                                                                                                                        <  2:827764/251‑1 (MQ=255)                      gtgCTGGAGATTGTCCTGGCAATCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAcc                                                                                                                                                                                                                         <  2:263692/251‑1 (MQ=255)                             agaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAATCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGaa                                                                                                                                                                                                                  <  1:16643/251‑1 (MQ=255)                                    ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTcac                                                                                                                                                                                                            <  1:699800/250‑1 (MQ=255)                                       ggCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaat                                                                                                                                                                                                         <  1:862926/250‑1 (MQ=255)                                        gCAAGCCAGGACACACCGCGTAAGCAAAAACCTTGTTTAACGCGCATTAAACGTAATCTCTAACGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCAAACAAAAATACGCGCCATGTGAATTTCCCTCTTGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTTATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTTAAATTCACaataa                                                                                                                                                                                                       <  2:375728/251‑1 (MQ=255)                                         cAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATACGACTATTCCTGACTATTATTGCTAAGTTCTTTTAAACCCGCCCGTTAATTAACTCCCCAGCTGAAATTCACattaat                                                                                                                                                                                                      >  2:183955/1‑251 (MQ=255)                                          aaGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAAtt                                                                                                                                                                                                     <  1:29695/251‑1 (MQ=255)                                                  aCGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCTTAGAGTTCGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGata                                                                                                                                                                                             <  2:848709/251‑1 (MQ=255)                                                     aaCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAAGTCACAATAATTAAGTGATATCg                                                                                                                                                                                          <  2:289382/251‑1 (MQ=255)                                                       ccGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAc                                                                                                                                                                                        <  1:315910/251‑1 (MQ=255)                                                        cgcgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGACTTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTTAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTAACCAGCTGAAATTCCCAATAATTAAGTGATATCGaa                                                                                                                                                                                        >  2:22558/1‑249 (MQ=255)                                                            tAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAgcgc                                                                                                                                                                                   >  1:84650/1‑251 (MQ=255)                                                                   aaaCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                                                                                            <  1:850668/251‑1 (MQ=255)                                                                   aaaCCTGGTTTAACACGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                                                                                            <  2:766889/251‑1 (MQ=255)                                                                      ccTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCTCGGCGTGCTCAGCCGTAGCTAAGCCGAGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGATATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATAGCGACATCGCGTTTTTGCat                                                                                                                                                                         <  2:780517/251‑1 (MQ=255)                                                                         ggTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCatta                                                                                                                                                                       <  1:332828/250‑1 (MQ=255)                                                                              aaCGGGCATTTAACGTAAACTCCATCGCCGCGGCGTGCTTAGCCGTAGCTAAGCCGCGAAAAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTGTTTGCATTATTTTGt                                                                                                                                                                 <  2:346025/251‑1 (MQ=255)                                                                                           cGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATgcag                                                                                                                                                        >  2:853206/1‑245 (MQ=255)                                                                                           cGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgcagagat                                                                                                                                                    <  1:453537/251‑7 (MQ=255)                                                                                                 cgctATCGCCGCGGCGTGCTCAGCCGTAGCGAAGCCGCGAACAAAAAGACGCGCCAGGTTAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAAtt                                                                                                                                               <  2:667460/248‑1 (MQ=255)                                                                                                   ctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGCGTACGGGAATGGACATCAATATACTTAAAGGAAATAAGCCTATTCCTGACTATTATTGTTACGTTCTTTTACACCCGCCCGTTAATTAACTCACCAGCTGAAATTTACAATACTTATGTGATATCGCCAGCGCGTTTTTGCATTATTTTGTTACATGCCAAGATGAATTGcc                                                                                                                                            >  2:865544/1‑251 (MQ=255)                                                                                                           gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATc                                                                                                                                                                                           <  1:246892/196‑1 (MQ=255)                                                                                                           gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATc                                                                                                                                                                                           >  2:246892/1‑196 (MQ=255)                                                                                                             ggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTaaca                                                                                                                                   >  1:169844/1‑250 (MQ=255)                                                                                                             ggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGAGAATTTCCATCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGCTACGTTCTTTTAAACCCGCCCGTTAATTAAATCACCAGCTGTAATTCACAATAATTAAGTGATATCGAAAGCTCGTTTTTGCATTATTTTGTTACATACAGAGATGAATTgacgcttaaact                                                                                                                                   >  2:118913/1‑239 (MQ=255)                                                                                                                     cAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATGTCCCTCTGGCGCGTAGAGTACGTGAATGGTCATCAATATCCTTAAAGTAAATAAGACTATTCCTGTCTATTATGGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGGGATATCGGCAGCGCGTTTTTGCAGTATGTTGTTACGTGCAGAGATGAATTGCCGATTGAACAAGCACtt                                                                                                                            <  2:846675/249‑1 (MQ=255)                                                                                                                                     gcgAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCa                                                                                                           >  1:64366/1‑250 (MQ=255)                                                                                                                                       gAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAc                                                                                                        >  2:193337/1‑251 (MQ=255)                                                                                                                                                 aCGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCa                                                                                               <  1:865540/250‑1 (MQ=255)                                                                                                                                                         ggTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGcc                                                                                      <  1:786963/251‑1 (MQ=255)                                                                                                                                                                  cccTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGAc                                                                              <  1:603214/250‑1 (MQ=255)                                                                                                                                                                             gTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACt                                                                  <  1:818029/251‑1 (MQ=255)                                                                                                                                                                                                  ccatatCCTTAAAGTAAATAAGACTAGTCCTGACTATTATTGATATGTTCTTTTAAACCCGCCCGTTAAGTAACTCACCAGATGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCCCTGTTCTTTGCTTTTGCGCTAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTc                                             <  2:226494/249‑1 (MQ=255)                                                                                                                                                                                                    atatCCTTCAAGTAAATAAGACTATTCCTGCCTATTATTGATAAGTTCTTTTAAACCCGACCGTTAATTAAATCACCAGCTGCAATTAACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGAGATGAATAACCGATTTAACAAGCACTTTTCTTTCCTTTTGCGCAAACCCGCTAGCATCAAGCGCCACACAGACGCAACATGGACTGTTAACCGGGAAAGatatttgct                                           >  2:80093/1‑246 (MQ=255)                                                                                                                                                                                                      atccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCTCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTa                                          <  1:690524/247‑1 (MQ=255)                                                                                                                                                                                                                  aaTAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcgc                             >  1:871222/1‑251 (MQ=255)                                                                                                                                                                                                                             tCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTACCAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACATACGCAACAGGGCCTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcgcaggaacgga                    >  2:873789/1‑240 (MQ=255)                                                                                                                                                                                                                             tCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcgcg                            <  1:873789/241‑1 (MQ=255)                                                                                                                                                                                                                                tGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGac                >  1:714039/1‑250 (MQ=255)                                                                                                                                                                                                                                       tattGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTtat         <  1:193337/250‑1 (MQ=255)                                                                                                                                                                                                                                        attGATAAATTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACTCTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTgtat       <  2:397238/251‑1 (MQ=255)                                                                                                                                                                                                                                             tAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTACCCAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGCCACTTTATATTTGCt  >  2:594036/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      CACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCT  >  NC_000913/88904‑89390

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