Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,903 A→C 9.8% intergenic (+23/+17) 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_0009131,074,9030AC 9.8% 132.7 / 7.8 41intergenic (+23/+17)rutR/putArut operon transcriptional repressor for/fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (20/17);  new base C (2/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

CGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGT  >  NC_000913/1074763‑1075051
                                                                                                                                            |                                                                                                                                                    
cGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                                                <  1:1073886/147‑1 (MQ=255)
cGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                                                >  2:1073886/1‑147 (MQ=255)
       gCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACa                                                                                                                                                      <  2:1437616/134‑1 (MQ=255)
       gCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACa                                                                                                                                                      >  1:1437616/1‑134 (MQ=255)
         cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                                                >  3:169396/1‑138 (MQ=255)
         cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                                                <  4:169396/138‑1 (MQ=255)
         cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACAtcctc                                                                                                                                                 <  4:1319813/137‑1 (MQ=255)
         cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACAtcctc                                                                                                                                                 >  3:1319813/1‑137 (MQ=255)
                        ggTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCgg                                                                                                                                              >  3:505106/1‑125 (MQ=255)
                        ggTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCgg                                                                                                                                              <  4:505106/125‑1 (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                                                                                                                            >  2:259999/1‑130 (MQ=255)
                                     aCGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCtata                                                                                                                            <  1:259999/130‑1 (MQ=255)
                                               aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                                                            >  1:908379/1‑136 (MQ=255)
                                               aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                                                            <  2:908379/136‑1 (MQ=255)
                                                                       ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                                                        <  4:60783/116‑1 (MQ=255)
                                                                       ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                                                        >  3:60783/1‑116 (MQ=255)
                                                                         gAAAACGTGCGGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCt                                                                     >  3:89396/1‑149 (MQ=255)
                                                                             aCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAAc                                                                                                                                 <  3:292007/85‑1 (MQ=255)
                                                                             aCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAAc                                                                                                                                 >  4:292007/1‑85 (MQ=255)
                                                                                               ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCAt                                                                                                                       >  4:48158/1‑77 (MQ=255)
                                                                                               ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCAt                                                                                                                       <  3:48158/77‑1 (MQ=255)
                                                                                                  aGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCta                                                                                                                              <  3:736094/67‑1 (MQ=255)
                                                                                                  aGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCta                                                                                                                              >  4:736094/1‑67 (MQ=255)
                                                                                                    gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                <  2:240087/143‑1 (MQ=255)
                                                                                                    gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                >  1:240087/1‑143 (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)
                                                                                                                                    gTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaa          >  2:435367/1‑149 (MQ=255)
                                                                                                                                     ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                         <  2:1060133/53‑1 (MQ=255)
                                                                                                                                     ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                         >  1:1060133/1‑53 (MQ=255)
                                                                                                                                     ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                           <  2:1102167/115‑1 (MQ=255)
                                                                                                                                     ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                           >  1:1102167/1‑115 (MQ=255)
                                                                                                                                        tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                            >  2:789841/1‑127 (MQ=255)
                                                                                                                                        tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                            <  1:789841/127‑1 (MQ=255)
                                                                                                                                          aCATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAAtt    <  1:435367/149‑1 (MQ=255)
                                                                                                                                            aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt  <  2:708150/149‑1 (MQ=255)
                                                                                                                                            aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt  >  2:351685/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGT  >  NC_000913/1074763‑1075051

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: