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%	191.8 / NA	53	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 (28/25);  total (28/25)

CTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGATAACGA  >  NC_000913/2905398‑2905651                                                                                                                                       |                                                                                                                        ctcGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTcaata                                                                                                                     >  1:321930/1‑139 (MQ=255)       aaTCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaaa                                                                                                               <  1:396005/139‑1 (MQ=255)         tCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             <  2:1164942/139‑1 (MQ=255)         tCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             <  2:321930/139‑1 (MQ=255)           gggTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGa                                                                                                           >  2:322098/1‑139 (MQ=255)             gTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                         <  1:1207679/139‑1 (MQ=255)             gTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                         <  2:1117541/139‑1 (MQ=255)                cTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAAc                                                                                                      >  1:963549/1‑139 (MQ=255)                           tctcTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAAtt                                                                                                                  >  1:1122265/1‑116 (MQ=255)                           tctcTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAAtt                                                                                                                  <  2:1122265/116‑1 (MQ=255)                           tctcTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        <  2:650589/126‑1 (MQ=255)                           tctcTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        >  1:650589/1‑126 (MQ=255)                             tctGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             >  1:698799/1‑119 (MQ=255)                             tctGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             <  2:698799/119‑1 (MQ=255)                             tctGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAg                                                                                                          <  1:1135145/122‑1 (MQ=255)                             tctGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAg                                                                                                          >  2:1135145/1‑122 (MQ=255)                                            aTTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCATGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                          <  2:963549/139‑1 (MQ=255)                                             ttGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAAt                                                                         >  1:1331766/1‑139 (MQ=255)                                                ttAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        >  2:180524/1‑105 (MQ=255)                                                ttAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        <  1:180524/105‑1 (MQ=255)                                                ttAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGa                                                                      >  1:1060824/1‑139 (MQ=255)                                                  aaTTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                                >  1:392266/1‑95 (MQ=255)                                                  aaTTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaa                                                                                                                <  2:392266/95‑1 (MQ=255)                                                                  gCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCATGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACa                                                    >  2:489704/1‑139 (MQ=255)                                                                  gCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACa                                                    >  2:1344379/1‑139 (MQ=255)                                                                       aacaaTCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATa                                                                                 <  2:178582/105‑1 (MQ=255)                                                                       aacaaTCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATa                                                                                 >  1:178582/1‑105 (MQ=255)                                                                                   ttattaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             <  2:280558/65‑1 (MQ=255)                                                                                   ttattaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             >  1:280558/1‑65 (MQ=255)                                                                                    tattaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTAtattt                                  >  2:1148902/1‑139 (MQ=255)                                                                                     attaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             <  2:363137/63‑1 (MQ=255)                                                                                     attaACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             >  1:363137/1‑63 (MQ=255)                                                                                           aTTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        <  2:402124/62‑1 (MQ=255)                                                                                           aTTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        >  1:402124/1‑62 (MQ=255)                                                                                           aTTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACa                                                    <  1:415491/114‑1 (MQ=255)                                                                                           aTTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACa                                                    >  2:415491/1‑114 (MQ=255)                                                                                              cTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             <  2:924266/54‑1 (MQ=255)                                                                                              cTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAgg                                                                                                             >  1:924266/1‑54 (MQ=255)                                                                                              cTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGa                                                                                                           <  2:961357/56‑1 (MQ=255)                                                                                              cTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGa                                                                                                           >  1:961357/1‑56 (MQ=255)                                                                                               tGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                         <  2:292992/57‑1 (MQ=255)                                                                                               tGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGt                                                                                                         >  1:292992/1‑57 (MQ=255)                                                                                                gTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        <  1:457051/57‑1 (MQ=255)                                                                                                gTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTa                                                                                                        >  2:457051/1‑57 (MQ=255)                                                                                                   aCAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaaa                                                                                                               <  2:1094541/47‑1 (MQ=255)                                                                                                   aCAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTaaa                                                                                                               >  1:1094541/1‑47 (MQ=255)                                                                                                   aCAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTAtt                    <  1:147620/138‑1 (MQ=255)                                                                                                   aCAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTAtt                    >  2:147620/1‑138 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAAt                                                                <  1:1139365/86‑1 (MQ=255)                                                                                                           gTTTAAAGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAAt                                                                >  2:1139365/1‑86 (MQ=255)                                                                                                                aaGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGATa      >  1:1020944/1‑139 (MQ=255)                                                                                                                 aGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                          <  1:911648/70‑1 (MQ=255)                                                                                                                 aGAAAAATGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGaaaa                                                                          >  2:911648/1‑70 (MQ=255)                                                                                                                    aaaaTGCGAAACATATCGTCAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGATAACGa  >  2:664348/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                        CTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGATAACGA  >  NC_000913/2905398‑2905651

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