Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,446,972 Δ1 bp coding (312/906 nt) feaR ← transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,9690G.95.6% 192.0 / 3.1 45coding (315/906 nt)feaRtranscriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (21/22);  total (22/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 = 9.28e-01

ATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGCCTGACCACTAAGCTGAAAAAC  >  NC_000913/1446877‑1447064
                                                                                            |                                                                                               
aTAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGa                                                                                         <  1:10951/100‑1 (MQ=255)
 tAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAg                                                                                        <  1:1086035/100‑1 (MQ=255)
 tAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAg                                                                                        <  1:299756/100‑1 (MQ=255)
  aGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAgg                                                                                       >  2:735082/1‑100 (MQ=255)
    gTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCa                                                                                     >  1:1237289/1‑100 (MQ=255)
    gTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCa                                                                                     >  1:656413/1‑100 (MQ=255)
       tttGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCg                                                                                  >  2:873222/1‑100 (MQ=255)
       tttGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCg                                                                                  <  2:486276/100‑1 (MQ=255)
       tttGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCg                                                                                  <  1:774464/100‑1 (MQ=255)
         tGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAg                                                                                >  2:945724/1‑100 (MQ=255)
         tGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAg                                                                                <  1:627801/100‑1 (MQ=255)
              gggAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACg                                                                           >  2:118152/1‑100 (MQ=255)
               ggAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGt                                                                          >  2:1184521/1‑100 (MQ=255)
                 naaaTATTGTTCCAGCAGTGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTaa                                                                        >  2:807704/2‑100 (MQ=255)
                     tatTGGTCCAGCATATTGCGTCGCACAAGTAATGAGATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACACGGCGTGAGGCATCGAGTAACGTa                                                                         <  2:1246048/95‑1 (MQ=255)
                          ttCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGCCTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAg                                                               <  1:216883/100‑1 (MQ=255)
                           tCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAgc                                                              <  1:1123279/100‑1 (MQ=255)
                            ccAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAgcg                                                             <  1:40965/100‑1 (MQ=255)
                            ccAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAgcg                                                             >  2:850146/1‑100 (MQ=255)
                             cagcagAGTCCGTGGCAAAAGTAATGAAAACTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAgcg                                                             <  1:996090/99‑1 (MQ=255)
                              agcagAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGcc                                                           >  2:392102/1‑100 (MQ=255)
                                        gTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCAcc                                                 <  2:437377/100‑1 (MQ=255)
                                        gTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCAcc                                                 <  2:295175/100‑1 (MQ=255)
                                         tGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCt                                                >  2:1144810/1‑100 (MQ=255)
                                         tGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCt                                                >  2:670631/1‑100 (MQ=255)
                                         tGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCt                                                >  1:728068/1‑100 (MQ=255)
                                          ggCAAAAGTAATGAAATNTGTNTAGAAGACTCCTGCCAGTAAAGCGAACANGGCGTGAGGCATCGAGTAACGNAATATCNCCAGCGCCAATCTGCACCTg                                               <  2:675739/100‑1 (MQ=255)
                                          ggCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTg                                               >  2:738875/1‑100 (MQ=255)
                                          ggCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTg                                               >  2:910611/1‑100 (MQ=255)
                                            cAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGa                                              >  1:462913/1‑100 (MQ=255)
                                                   aaTGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATc                                      >  2:318225/1‑100 (MQ=255)
                                                      gAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTg                                   >  1:491408/1‑100 (MQ=255)
                                                      gAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTg                                   >  2:251867/1‑100 (MQ=255)
                                                       aaaTCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGc                                  <  1:945732/100‑1 (MQ=255)
                                                              tttAGAAGACTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCatta                           <  1:236350/100‑1 (MQ=255)
                                                                 agaagaCTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTg                        <  2:446142/100‑1 (MQ=255)
                                                                  gaagaCTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGc                       >  2:605700/1‑100 (MQ=255)
                                                                   aagaCTCCTGCCAGTAAAGCGAACAGGTCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGcc                      <  2:1258392/100‑1 (MQ=255)
                                                                   aagaCTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGcc                      <  2:706849/100‑1 (MQ=255)
                                                                   aagaCTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGcc                      <  1:1114181/100‑1 (MQ=255)
                                                                   aagaCTCCTGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATAGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGcc                      <  2:723498/100‑1 (MQ=255)
                                                                         ccTGCCAGTAAAGCGAACAGGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGCCTGAcc                 <  1:117868/100‑1 (MQ=255)
                                                                           tGCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCATCTGACGCTCATCCTGCTCCATTATTGCCTGACCACt              >  1:76765/1‑100 (MQ=255)
                                                                            gCCAGTAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGCCTGACCACTa             <  1:404089/100‑1 (MQ=255)
                                                                                 tAAAGCGAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGCCTGACCACTAAGCTg        <  2:284718/100‑1 (MQ=255)
                                                                                       gAACAGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGCCTGACCACTAAGCTGAAAAAc  >  1:682229/1‑100 (MQ=255)
                                                                                            |                                                                                               
ATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAAATCTGTTTAGAAGACTCCTGCCAGTAAAGCGAACAGGGGCGTGAGGCATCGAGTAACGTAATATCGCCAGCGCCAATCTGCACCTGACGCTCATCCTGCTCCATTATTGCCTGACCACTAAGCTGAAAAAC  >  NC_000913/1446877‑1447064

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: