Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,899 T→C 8.7% intergenic (+19/+21) 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,8990TC 8.7% 139.7 / 4.5 46intergenic (+19/+21)rutR/putArut operon transcriptional repressor for/fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base C (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-02

CAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAA  >  NC_000913/1074752‑1075043
                                                                                                                                                   |                                                                                                                                                
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                                                                                                                                              >  1:1437616/1‑134 (MQ=255)
                  gCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACa                                                                                                                                              <  2:1437616/134‑1 (MQ=255)
                    cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                                        <  4:169396/138‑1 (MQ=255)
                    cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGATGTTACAACCtcctcc                                                                                                                                        >  3:169396/1‑138 (MQ=255)
                    cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACAtcctc                                                                                                                                         >  3:1319813/1‑137 (MQ=255)
                    cccTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACAtcctc                                                                                                                                         <  4:1319813/137‑1 (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                                                                                                    <  2:908379/136‑1 (MQ=255)
                                                          aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                                                    >  1:908379/1‑136 (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                                                                                                                      >  4:736094/1‑67 (MQ=255)
                                                                                                             aGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCta                                                                                                                      <  3:736094/67‑1 (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)
                                                                                                                                               gTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaa  >  2:435367/1‑149 (MQ=255)
                                                                                                                                                ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                 >  1:1060133/1‑53 (MQ=255)
                                                                                                                                                ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                 <  2:1060133/53‑1 (MQ=255)
                                                                                                                                                ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                   >  1:1102167/1‑115 (MQ=255)
                                                                                                                                                ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                   <  2:1102167/115‑1 (MQ=255)
                                                                                                                                                   tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                    <  1:789841/127‑1 (MQ=255)
                                                                                                                                                   tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                    >  2:789841/1‑127 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CAACATTACGCCGATTTCGCCCCTCAGGTGGAGGCGGTGACAGGCGCGACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAA  >  NC_000913/1074752‑1075043

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: ‑