Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 933,022 A→G D143G (GAC→GGC)  lrp → leucine‑responsive global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913933,0220AG96.4% 101.4 / ‑5.3 28D143G (GAC→GGC) lrpleucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/16);  minor base T (0/1);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCGCTGTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACC  >  NC_000913/932875‑933164
                                                                                                                                                   |                                                                                                                                              
aTACCGCTGTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGgcac                                                                                                                                             <  1:7613/151‑1 (MQ=255)
      cTGTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGa                                                                                                                                        <  1:493736/150‑1 (MQ=255)
        gTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACAt                                                                                                                                     >  2:291952/1‑151 (MQ=255)
                        aaaTTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGgaag                                                                                                                      <  2:437080/150‑1 (MQ=255)
                        aaaTTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGgaa                                                                                                                       <  1:393541/149‑1 (MQ=255)
                          aTTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGaagaag                                                                                                                   >  1:274594/1‑151 (MQ=255)
                                    gtCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCagag                                                                                                         >  2:188244/1‑151 (MQ=255)
                                                               aCCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTa                                                                              >  1:352657/1‑151 (MQ=255)
                                                               aCCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTa                                                                              >  1:45724/1‑151 (MQ=255)
                                                                    ttGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAacac                                                                          <  2:394626/150‑1 (MQ=255)
                                                                        aaaCGCGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGGCTGCCTGGCGTCTATGTCACACGGACATACGTTGTTATGGAAGAAGTCAAGCATAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAAc                                                                     <  1:533493/151‑1 (MQ=255)
                                                                        aaaCACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAAc                                                                     <  1:502592/151‑1 (MQ=255)
                                                                        aaaCACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAAc                                                                     <  1:180213/151‑1 (MQ=255)
                                                                                 tGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCaaa                                                            >  2:169830/1‑151 (MQ=255)
                                                                                   ccGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAAt                                                          <  2:151593/151‑1 (MQ=255)
                                                                                               gCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTg                                              >  2:11111/1‑151 (MQ=255)
                                                                                                              cTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTa                               <  1:513953/151‑1 (MQ=255)
                                                                                                                   ggAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATcc                           <  1:262997/150‑1 (MQ=255)
                                                                                                                       aCCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACa                      >  2:81804/1‑151 (MQ=255)
                                                                                                                          ctgctgCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCa                   >  1:60961/1‑151 (MQ=255)
                                                                                                                           tgctgcGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCaa                  >  2:103722/1‑151 (MQ=255)
                                                                                                                                  tCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAc            <  2:364028/150‑1 (MQ=255)
                                                                                                                                    tGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgg         <  1:506352/151‑1 (MQ=255)
                                                                                                                                      ccTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg       <  2:274599/151‑1 (MQ=255)
                                                                                                                                      ccTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg       <  2:436264/151‑1 (MQ=255)
                                                                                                                                      ccTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg       <  1:304619/151‑1 (MQ=255)
                                                                                                                                      ccTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg       <  1:23637/151‑1 (MQ=255)
                                                                                                                                          gCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTaa    >  1:386992/1‑150 (MQ=255)
                                                                                                                                           cGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAcc  <  2:20671/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ATACCGCTGTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACC  >  NC_000913/932875‑933164

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