Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,905,532	T→C	intergenic (‑114/‑25)	queE ← / → yqcG	7‑carboxy‑7‑deazaguanine synthase, queosine biosynthesis/membrane stress resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,905,532	0	T	C	100.0%	154.5 / NA	43	intergenic (‑114/‑25)	queE/yqcG	7‑carboxy‑7‑deazaguanine synthase, queosine biosynthesis/membrane stress resistance protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/21);  total (22/21)

CTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGAT  >  NC_000913/2905398‑2905646                                                                                                                                       |                                                                                                                   ctcGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTcaata                                                                                                                >  2:1237127/1‑139 (MQ=255)      tAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                           <  2:489913/139‑1 (MQ=255)                      aGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGtatta                                                                                           >  2:581929/1‑139 (MQ=255)                           tctcTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCa                                                                                      >  2:224520/1‑139 (MQ=255)                             tctGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCaga                                                                                    >  1:1098944/1‑139 (MQ=255)                                  gAAGTGGATAATTGTTAATTATTGCAGATCTTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                    <  2:111776/118‑1 (MQ=255)                                  gAAGTGGATAATTGTTAATTATTGCAGATCTTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                    >  1:111776/1‑118 (MQ=255)                                  gAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                    >  2:779389/1‑118 (MQ=255)                                  gAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                    <  1:779389/118‑1 (MQ=255)                                           aaTTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaa                                                                      <  2:1098944/139‑1 (MQ=255)                                                ttAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                     >  1:162461/1‑135 (MQ=255)                                                ttAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                     <  2:162461/135‑1 (MQ=255)                                                ttAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGa                                                                 >  2:727643/1‑139 (MQ=255)                                                  aaTTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGAtt                                                               >  2:396157/1‑139 (MQ=255)                                                      attGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAAt                                                           <  1:595201/139‑1 (MQ=255)                                                      attGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAAt                                                           <  2:138865/139‑1 (MQ=255)                                                        tgaaGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                           <  2:1064729/86‑1 (MQ=255)                                                        tGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                           >  1:1064729/1‑89 (MQ=255)                                                               ccTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGAt                                                                <  2:391120/125‑1 (MQ=255)                                                               ccTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGAt                                                                >  1:391120/1‑125 (MQ=255)                                                               ccTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCaaa                                                  >  2:1465122/1‑139 (MQ=255)                                                                  gCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaaa                                                                                                          <  1:1333769/80‑1 (MQ=255)                                                                  gCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaaa                                                                                                          >  2:1333769/1‑80 (MQ=255)                                                                         caaTCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCa                                        <  1:1465122/139‑1 (MQ=255)                                                                            tCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGGGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCa                                                    >  2:484871/1‑124 (MQ=255)                                                                            tCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGGGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCa                                                    <  1:484871/124‑1 (MQ=255)                                                                                      ttaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                     <  2:1123707/97‑1 (MQ=255)                                                                                      ttaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                     >  1:1123707/1‑97 (MQ=255)                                                                                      ttaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTAtatttat                           >  1:957804/1‑139 (MQ=255)                                                                                             tCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCt                    >  1:724815/1‑139 (MQ=255)                                                                                             tCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCt                    >  2:264853/1‑139 (MQ=255)                                                                                                gTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTa                 <  1:581929/139‑1 (MQ=255)                                                                                                   aCAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCa                                                         <  1:757668/96‑1 (MQ=255)                                                                                                   aCAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCa                                                         >  2:757668/1‑96 (MQ=255)                                                                                                          ggTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTa       <  1:396157/139‑1 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                           >  1:756558/1‑38 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                           <  2:756558/38‑1 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGAt                                                                <  2:298384/81‑1 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGAt                                                                >  1:298384/1‑81 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTaa      >  2:324817/1‑139 (MQ=255)                                                                                                             ttAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATg    <  2:724815/139‑1 (MQ=255)                                                                                                               aaaGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGAt  <  1:264853/139‑1 (MQ=255)                                                                                                                     aaaTGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                   >  1:98757/1‑36 (MQ=255)                                                                                                                     aaaTGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                   <  2:98757/36‑1 (MQ=255)                                                                                                                                       |                                                                                                                   CTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGAT  >  NC_000913/2905398‑2905646

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