Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,260,901 C→T 15.1% intergenic (‑98/+27) dauA ← / ← prs C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9010CT15.1% 139.0 / 17.6 53intergenic (‑98/+27)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (4/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

AGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTC  >  NC_000913/1260760‑1261045
                                                                                                                                             |                                                                                                                                                
agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCg                                                                                                                                            >  4:1045389/1‑148 (MQ=255)
                 cacaTGTGAGGAAAATATTTTGTTCACAATTATGTTTCATACGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAAGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                                                                           <  3:1045389/148‑1 (MQ=255)
                                           caATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCtt                                                                                                                                      >  3:847632/1‑111 (MQ=255)
                                           caATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCtt                                                                                                                                      <  4:847632/111‑1 (MQ=255)
                                             ataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATc                                                                                               >  4:273081/1‑148 (MQ=255)
                                                  gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAg                                                                                                                                              >  4:391321/1‑96 (MQ=255)
                                                  gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAg                                                                                                                                              <  3:391321/96‑1 (MQ=255)
                                                        tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                                                                    >  4:902745/1‑148 (MQ=255)
                                                            tGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGa                                                                                                                                   <  1:1074277/97‑1 (MQ=255)
                                                            tGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGa                                                                                                                                   >  2:1074277/1‑97 (MQ=255)
                                                                  catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                                                                          >  4:935501/1‑84 (MQ=255)
                                                                  catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                                                                          <  3:935501/84‑1 (MQ=255)
                                                                          gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAGACCCGCTGCGGCGGGTttttttg                                                                                                                                  <  3:785211/84‑4 (MQ=37)
                                                                          gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                                                             <  1:796139/105‑1 (MQ=255)
                                                                          gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                                                             >  2:796139/1‑105 (MQ=255)
                                                                          gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACACACAAAAAGACCCGCTGCGGCGGGTttttttg                                                                                                                                  >  4:785211/1‑81 (MQ=25)
                                                                           tCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCt                                                                 <  3:273081/148‑1 (MQ=255)
                                                                            cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                           <  3:430188/121‑1 (MQ=255)
                                                                            cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                           >  4:430188/1‑121 (MQ=255)
                                                                              aCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtgt                                                                                                                 >  4:592578/1‑97 (MQ=255)
                                                                              aCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtgt                                                                                                                 <  3:592578/97‑1 (MQ=255)
                                                                                      gCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGc                                                                                                                                 >  4:288490/1‑73 (MQ=255)
                                                                                      gCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGc                                                                                                                                 <  3:288490/73‑1 (MQ=255)
                                                                                           tGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTg                                                                                >  2:695452/1‑117 (MQ=255)
                                                                                           tGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTg                                                                                <  1:695452/117‑1 (MQ=255)
                                                                                                gaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAgg                                            <  3:902745/148‑1 (MQ=255)
                                                                                                    cATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                        <  4:917529/100‑1 (MQ=255)
                                                                                                    cATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                        >  3:917529/1‑100 (MQ=255)
                                                                                                       aGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                                                                      <  1:556258/67‑1 (MQ=255)
                                                                                                       aGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                                                                      >  2:556258/1‑67 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAg                                                                                                                    <  3:1036692/66‑1 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAg                                                                                                                    >  4:1036692/1‑66 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                                                                <  3:468542/70‑1 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                                                                >  4:468542/1‑70 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGaa                                                                      <  4:862439/112‑1 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGaa                                                                      >  3:862439/1‑112 (MQ=255)
                                                                                                             tACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc                                                        >  2:34973/1‑123 (MQ=255)
                                                                                                             tACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc                                                        <  1:34973/123‑1 (MQ=255)
                                                                                                                             aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                        <  1:331801/75‑1 (MQ=255)
                                                                                                                             aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                        >  2:331801/1‑75 (MQ=255)
                                                                                                                               aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGt                                          >  4:467923/1‑119 (MQ=255)
                                                                                                                               aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGt                                          <  3:467923/119‑1 (MQ=255)
                                                                                                                                cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc                                                            >  2:660968/1‑100 (MQ=255)
                                                                                                                                cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc                                                            <  1:660968/100‑1 (MQ=255)
                                                                                                                                  caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCag                                       >  3:496996/2‑119 (MQ=255)
                                                                                                                                  caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCag                                       <  4:496996/118‑1 (MQ=255)
                                                                                                                                  aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACg                             >  3:42155/1‑129 (MQ=255)
                                                                                                                                  aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACg                             <  4:42155/129‑1 (MQ=255)
                                                                                                                                    aagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCag                                       <  2:930620/114‑1 (MQ=255)
                                                                                                                                    aagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCag                                       >  1:930620/4‑117 (MQ=255)
                                                                                                                                     agaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                        >  2:958417/3‑67 (MQ=37)
                                                                                                                                     agaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                        <  1:958417/65‑1 (MQ=37)
                                                                                                                                     aaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCAt       >  4:409357/1‑148 (MQ=255)
                                                                                                                                          cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACg                                 <  2:302455/117‑1 (MQ=255)
                                                                                                                                          cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACg                                 >  1:302455/1‑117 (MQ=255)
                                                                                                                                          cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCgg                        >  4:334391/1‑126 (MQ=255)
                                                                                                                                          cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCgg                        <  3:334391/126‑1 (MQ=255)
                                                                                                                                          cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTc  >  1:424941/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                
AGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTC  >  NC_000913/1260760‑1261045

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: