Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,074,894	G→T	13.0%	intergenic (+14/+26)	rutR → / ← putA	rut operon transcriptional repressor for/fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,074,894	0	G	T	13.0%	147.3 / 12.8	46	intergenic (+14/+26)	rutR/putA	rut operon transcriptional repressor for/fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base T (3/3);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

CCACTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCA  >  NC_000913/1074747‑1075022                                                                                                                                                    |                                                                                                                                 ccACTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGaggagg                                                                                                                                 <  3:953799/149‑1 (MQ=255)  cACTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGt                                                                                                                                >  4:366454/1‑149 (MQ=255)   aCTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGtt                                                                                                                               <  3:449319/149‑1 (MQ=255)   aCTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGtt                                                                                                                               <  3:366454/149‑1 (MQ=255)    cTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGgaggag                                                                                                                                  <  4:971122/145‑1 (MQ=255)    cTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGgaggag                                                                                                                                  >  3:971122/1‑145 (MQ=255)      cAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTa                                                                                                                            <  3:119313/149‑1 (MQ=255)       aaCATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTaa                                                                                                                           >  4:655424/1‑149 (MQ=255)       aaCATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTaa                                                                                                                           <  1:187784/149‑1 (MQ=255)       aaCATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTaa                                                                                                                           <  3:655424/149‑1 (MQ=255)                 cGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                   <  1:1073886/147‑1 (MQ=255)                 cGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                   >  2:1073886/1‑147 (MQ=255)                 cGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAAc                                                                                                                          <  4:614295/140‑1 (MQ=255)                 cGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAAc                                                                                                                          >  3:614295/1‑140 (MQ=255)                        gCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACa                                                                                                                         <  2:1437616/134‑1 (MQ=255)                        gCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACa                                                                                                                         >  1:1437616/1‑134 (MQ=255)                          cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                   <  4:169396/138‑1 (MQ=255)                          cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                   >  3:169396/1‑138 (MQ=255)                          cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACAtcctc                                                                                                                    <  4:1319813/137‑1 (MQ=255)                          cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACAtcctc                                                                                                                    >  3:1319813/1‑137 (MQ=255)                                         ggTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCgg                                                                                                                 <  4:505106/125‑1 (MQ=255)                                         ggTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCgg                                                                                                                 >  3:505106/1‑125 (MQ=255)                                           tGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTaa                                                                                       >  1:76990/1‑149 (MQ=255)                                               aGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTg                                                                                   >  3:360705/1‑149 (MQ=255)                                                 gcgcGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGc                                                                                 <  2:76990/149‑1 (MQ=255)                                                     gACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTa                                                                             <  4:360705/149‑1 (MQ=255)                                                      aCGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCtata                                                                                               <  1:259999/130‑1 (MQ=255)                                                      aCGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCtata                                                                                               >  2:259999/1‑130 (MQ=255)                                                                aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                               >  1:908379/1‑136 (MQ=255)                                                                aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                               <  2:908379/136‑1 (MQ=255)                                                                                       gTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGttac                                                                                                                             <  2:471675/67‑3 (MQ=255)                                                                                       gTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGttac                                                                                                                             >  1:471675/1‑65 (MQ=255)                                                                                        ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                           >  3:60783/1‑116 (MQ=255)                                                                                        ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                           <  4:60783/116‑1 (MQ=255)                                                                                          gAAAACGTGCGGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCt                                        >  3:89396/1‑149 (MQ=255)                                                                                              aCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAAc                                                                                                    <  3:292007/85‑1 (MQ=255)                                                                                              aCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAAc                                                                                                    >  4:292007/1‑85 (MQ=255)                                                                                                                ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCAt                                                                                          <  3:48158/77‑1 (MQ=255)                                                                                                                ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCAt                                                                                          >  4:48158/1‑77 (MQ=255)                                                                                                                   aGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCta                                                                                                 <  3:736094/67‑1 (MQ=255)                                                                                                                   aGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCta                                                                                                 >  4:736094/1‑67 (MQ=255)                                                                                                                     gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                   >  1:240087/1‑143 (MQ=255)                                                                                                                     gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                   <  2:240087/143‑1 (MQ=255)                                                                                                                      ggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt            <  4:89396/149‑1 (MQ=255)                                                                                                                            gACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGc                                                                         <  3:582335/82‑1 (MQ=255)                                                                                                                            gACCACGTTAAAGACGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCAAACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTc      >  3:100581/1‑149 (MQ=255)                                                                                                                            gACCACGTTAAACATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGc                                                                         >  4:582335/1‑82 (MQ=255)                                                                                                                              ccACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGc    >  2:104733/1‑149 (MQ=255)                                                                                                                                aCGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa  >  1:708150/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                 CCACTCAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCA  >  NC_000913/1074747‑1075022

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