Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,097,658 T→G intergenic (‑210/‑78) rhaB ← / → rhaS rhamnulokinase/DNA‑binding transcriptional activator RhaS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,097,6580TG97.7% 137.7 / NA 43intergenic (‑210/‑78)rhaB/rhaSrhamnulokinase/DNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (13/29);  minor base A (1/0);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGT  >  NC_000913/4097533‑4097785
                                                                                                                             |                                                                                                                               
cATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACCCTGTATCTAGAAGAATGGACGTTTTATACGTGGGTTTCCGt                                                                                                                          >  2:145404/1‑133 (MQ=255)
 aTTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGACAAATCGAGGTTGTTGACGTGGGTTTCCGTc                                                                                                                         >  1:200899/1‑133 (MQ=255)
 aTTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTc                                                                                                                         >  1:444669/1‑133 (MQ=255)
 aTTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGAAGTTTTTGACGTGGATTTCCGTc                                                                                                                         >  1:192425/1‑133 (MQ=255)
         ccAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAAttt                                                                                                                 <  2:444669/133‑1 (MQ=255)
             ggAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGATTTGTACGTGGGTTTCTGTCGAAAATTTAAgg                                                                                                             >  2:190860/1‑133 (MQ=255)
               aaaGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTa                                                                                                           <  1:86708/133‑1 (MQ=255)
               aaaGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTa                                                                                                           <  1:323553/133‑1 (MQ=255)
                aaGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTaa                                                                                                          <  1:44600/133‑1 (MQ=255)
                      aaCGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCt                                                                                                    <  1:83763/133‑1 (MQ=255)
                         gTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTATTACGTGGGTTTCCGTGGAAAATTTGAGGTGAGAACCTGAc                                                                                                 >  2:136739/1‑133 (MQ=255)
                          tgatgaTGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGTATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGAcc                                                                                                <  2:88274/133‑1 (MQ=255)
                           gatgatGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCt                                                                                               >  1:338285/1‑133 (MQ=255)
                            atgatgTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTc                                                                                              <  2:400179/133‑1 (MQ=255)
                              gatgTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGt                                                                                            <  1:269305/133‑1 (MQ=255)
                                 gTTCACAACTTGCTGCATTGTGCTGATGTGCTGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGAt                                                                                         <  1:145404/133‑1 (MQ=255)
                                   tCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTa                                                                                       <  1:483724/133‑1 (MQ=255)
                                           tGCTGAATTGTGGTGTTGTGATGCTCACTGCATTTCCTGAATATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCg                                                                               <  2:200899/133‑1 (MQ=255)
                                           tGCTGAATTGAGGTGATGTGATGCTCACTGCACTTCCTGAATATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCg                                                                               <  2:192425/133‑1 (MQ=255)
                                            gCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGc                                                                              >  2:307106/1‑133 (MQ=255)
                                            gCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGGAAGAACCTGACATCGTGATTAATATTTCGc                                                                              >  1:166454/1‑133 (MQ=255)
                                                    gtggtgATGTGATGCTCACCGCAATTTCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTGACGTGGGTTTCCGTCGAAAATTTAAGGTAAGATCCTGACCTCGTGATGACTATTTCGCCGTGTTga                                                                      <  2:255465/133‑1 (MQ=255)
                                                     tggtgATGTGATGCTAACCGCATTTCCTCAAAATTCACGCTTTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTgac                                                                     <  2:166454/133‑1 (MQ=255)
                                                                  cTCACCGCTTTTCCTCAAGATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGc                                                        <  1:192753/133‑1 (MQ=255)
                                                                  cTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGc                                                        >  1:321289/1‑133 (MQ=255)
                                                                  cTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACATGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGc                                                        <  2:109154/133‑1 (MQ=255)
                                                                     aCCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAg                                                     >  1:487501/1‑133 (MQ=255)
                                                                        gCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTAt                                                  <  2:343110/133‑1 (MQ=255)
                                                                          aTTTCCTTAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                <  2:453075/133‑1 (MQ=255)
                                                                          aTTTCCTGTAAATTCCCGCTGTCTCTTGAATAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                <  1:136739/133‑1 (MQ=255)
                                                                          aTTTCCTGAAAATTCACGCTGTATTTTGCAAACTCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                <  2:191851/133‑1 (MQ=255)
                                                                          aTTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                <  2:338285/133‑1 (MQ=255)
                                                                          aTTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                <  1:278712/133‑1 (MQ=255)
                                                                          aTTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGa                                                <  2:278518/133‑1 (MQ=255)
                                                                            ttCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGAcc                                              <  2:332030/133‑1 (MQ=255)
                                                                            ttCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGAcc                                              <  1:22724/133‑1 (MQ=255)
                                                                            ttCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGAcc                                              <  1:5462/133‑1 (MQ=255)
                                                                            ttCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGAcc                                              <  1:269329/133‑1 (MQ=255)
                                                                            ttCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGAcc                                              >  1:114917/1‑133 (MQ=255)
                                                                                  aaaaTTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTa                                        <  2:280501/133‑1 (MQ=255)
                                                                                                         aaTCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTc                 >  1:89059/1‑133 (MQ=255)
                                                                                                         aaTCGACGTTTTTTACGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTc                 >  2:349890/1‑133 (MQ=255)
                                                                                                                      tACGTGGGTTTCCGTTGAAAATCTAAGGTACGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTcc    <  1:252492/133‑1 (MQ=255)
                                                                                                                        cGTGGGTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGt  <  1:349890/133‑1 (MQ=255)
                                                                                                                             |                                                                                                                               
CATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGT  >  NC_000913/4097533‑4097785

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