Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,155 C→A T11N (ACC→AAC)  fabR → transcriptional repressor of fabA and fabB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1550CA97.8% 152.5 / ‑3.0 45T11N (ACC→AAC) fabRtranscriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (28/16);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

AATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCT  >  NC_000913/4161007‑4161302
                                                                                                                                                    |                                                                                                                                                   
aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAacc                                                                                                                                                   >  1:3065699/1‑151 (MQ=255)
aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAacc                                                                                                                                                   <  2:2173203/151‑1 (MQ=255)
   tGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAAccgcc                                                                                                                                                >  2:4000963/1‑151 (MQ=255)
   tGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAAccgcc                                                                                                                                                >  2:1900566/1‑151 (MQ=255)
               gTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCtggtgg                                                                                                                                    >  2:3642276/1‑151 (MQ=255)
                ttttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCtggtgg                                                                                                                                    >  2:1089562/1‑150 (MQ=255)
                 tttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAAACCGCCGTTCGCTGGTGGa                                                                                                                                   <  2:1377052/150‑1 (MQ=255)
                     ttGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCg                                                                                                                              >  1:3328337/1‑151 (MQ=255)
                               gTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                                     >  1:1348132/1‑150 (MQ=255)
                                   attcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                                 >  2:2766956/1‑150 (MQ=255)
                                    ttcaAAATACTGGAGCAATCAAGTATTTTCATTCTCTGGTATGGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAACAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa                                                                                                                <  1:167395/150‑1 (MQ=255)
                                          aTACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                         >  2:1132845/1‑151 (MQ=255)
                                            aCTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTGGGAAAGGATTCAGACATGGTGATGGGCGTAAGAGCGAAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                                                                                                       <  1:3846437/151‑1 (MQ=255)
                                                        cAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc                                                                                           >  1:61707/1‑151 (MQ=255)
                                                          gTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa                                                                                          <  2:2378251/150‑1 (MQ=255)
                                                            aTGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                       >  2:4051360/1‑151 (MQ=255)
                                                              gTTCATTCTCTGGTATAGTGCCGGCAGAACTTTTGGCAAGGATTCAGACATCGTGGTGGGCGTAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                                     >  1:2392415/1‑151 (MQ=255)
                                                                tCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                       >  1:787479/1‑147 (MQ=255)
                                                                tCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                       <  2:787439/147‑1 (MQ=255)
                                                                tCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                                     >  1:3285183/1‑149 (MQ=255)
                                                                tCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                                     >  1:2283667/1‑149 (MQ=255)
                                                                      tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCg                                                                             >  1:1587118/1‑151 (MQ=255)
                                                                        tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATTGGCGCAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGt                                                                            >  1:478999/1‑150 (MQ=255)
                                                                           tataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGa                                                                          >  1:3370999/1‑149 (MQ=255)
                                                                            ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAg                                                                        >  2:2693658/1‑150 (MQ=255)
                                                                             taGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg                                                                      >  1:1066072/1‑151 (MQ=255)
                                                                               gTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgc                                                                    <  2:2794598/151‑1 (MQ=255)
                                                                                tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGGGCGAAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg                                                                      <  2:2889030/148‑1 (MQ=255)
                                                                                    agcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTg                                                               >  1:3603768/1‑151 (MQ=255)
                                                                                         tACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCg                                                          >  2:1560954/1‑151 (MQ=255)
                                                                                           cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCggg                                                        >  1:3595254/1‑151 (MQ=255)
                                                                                            ttttGGCAAGATTTCAGACATCGTGATGGGCGTAAGAGCCCAACAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGca                                                          <  1:2162893/148‑2 (MQ=255)
                                                                                                ggCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                                    <  2:34177/150‑1 (MQ=255)
                                                                                                ggCAAGGATTCAGACATCGTGAGGGGCGAAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAt                                                     <  1:28369/149‑1 (MQ=255)
                                                                                                                 cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTAt                                  >  1:3309583/1‑151 (MQ=255)
                                                                                                                  gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATc                                 >  1:2714691/1‑151 (MQ=255)
                                                                                                                      tGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGc                              >  2:3087332/1‑150 (MQ=255)
                                                                                                                           gTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                        >  1:328048/1‑151 (MQ=255)
                                                                                                                           gTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTc                         <  2:3065901/150‑1 (MQ=255)
                                                                                                                             aaGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                       >  1:117847/1‑150 (MQ=255)
                                                                                                                                    caacaaAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGt                                                                >  1:2178674/1‑102 (MQ=255)
                                                                                                                                    caacaaAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                 >  2:2019814/1‑149 (MQ=255)
                                                                                                                                     aacaaAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                        <  1:1907914/141‑1 (MQ=255)
                                                                                                                                     aacaaAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                   <  2:3285390/146‑1 (MQ=255)
                                                                                                                                     aaaaaaaaGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGt                                                                <  2:2178905/98‑1 (MQ=255)
                                                                                                                                     aaaaaaaaGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTa                <  1:1751442/146‑1 (MQ=255)
                                                                                                                                        aataaaaaaaaaACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                   <  1:773319/137‑1 (MQ=255)
                                                                                                                                               aaaaaaCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt    >  1:2285868/1‑151 (MQ=255)
                                                                                                                                                 aaaaCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt  <  2:428398/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCT  >  NC_000913/4161007‑4161302

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