Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,904,716	C→T	intergenic (+1342/+31)	ygcE → / ← queE	putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,904,716	0	C	T	100.0%	94.5 / NA	28	intergenic (+1342/+31)	ygcE/queE	putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/16);  total (12/16)

TTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTT  >  NC_000913/2904589‑2904840                                                                                                                                |                                                                                                                             ttCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCAccc                                                                                                                         <  1:14727/133‑1 (MQ=255)       ccAGCGAAAAATTAATGCCACAGAATTTGTAGAG‑ATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTtata                                                                                                                  <  1:565695/133‑1 (MQ=255)        cAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGTAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTtata                                                                                                                  <  2:343133/133‑1 (MQ=255)        cAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGTAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTtata                                                                                                                  >  1:343133/1‑133 (MQ=255)        cAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGAATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTtata                                                                                                                  <  2:607622/133‑1 (MQ=255)         aGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATaa                                                                                                                 >  2:490075/1‑133 (MQ=255)                 aTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTtata                                                                                                         <  2:342513/133‑1 (MQ=255)                       gCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGtt                                                                                                   <  1:419828/133‑1 (MQ=255)                          acaGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTaa                                                                                                <  1:490075/133‑1 (MQ=255)                               aTTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAgg                                                                                           <  2:599215/133‑1 (MQ=255)                               aTTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAgg                                                                                           <  2:542365/133‑1 (MQ=255)                               aTTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAgg                                                                                           <  2:459297/133‑1 (MQ=255)                                 ttGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCa                                                                                         <  2:53056/133‑1 (MQ=255)                                           tGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGat                                                                               <  1:504525/133‑1 (MQ=255)                                                     tAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATgtgt                                                                     >  2:598032/1‑133 (MQ=255)                                                     tAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATgtgt                                                                     >  1:578753/1‑133 (MQ=255)                                                        cccTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTg                                                                  >  2:553694/1‑133 (MQ=255)                                                             gATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCg                                                             >  1:147947/1‑133 (MQ=255)                                                                     aaaTGGAGGATTTTTTAGAGGCGACAAAAAATATAAATGAGCAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGc                                                     <  2:257051/133‑1 (MQ=255)                                                                           aggagTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTa                                               >  1:466733/1‑133 (MQ=255)                                                                               gTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTAcgcg                                           <  1:103813/133‑1 (MQ=255)                                                                                ttttCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACgcgc                                          >  1:440931/1‑133 (MQ=255)                                                                                           gACAAACAATATAAATGAGTAGAAGTCTGTACTGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGttt                               <  2:242009/133‑1 (MQ=255)                                                                                                              tAGAAGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGAAATGCAGATTTAAACGCACAATCGTGTGGc            >  1:206889/1‑133 (MQ=255)                                                                                                                  aGTCTTTACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGcatca        <  2:578753/133‑1 (MQ=255)                                                                                                                      tttACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGTATCAGTGCACAAAAGTGTGGCATCATCCt    >  2:186867/1‑133 (MQ=255)                                                                                                                        tACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCCCAATCGCGTGGCACCATCCttt  >  1:315583/1‑133 (MQ=255)                                                                                                                        tACAGGTTTCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCttt  >  2:28086/1‑133 (MQ=255)                                                                                                                                |                                                                                                                             TTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAATATTTAGATATTTATGTGTTTGCATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTT  >  NC_000913/2904589‑2904840

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