Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,907 A→C 14.0% intergenic (+14/+21) mog → / ← satP molybdopterin adenylyltransferase/acetate/succinate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,9070AC14.0% 65.2 / 8.4 41intergenic (+14/+21)mog/satPmolybdopterin adenylyltransferase/acetate/succinate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base C (2/4);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

TTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCG  >  NC_000913/9778‑10036
                                                                                                                                  |                                                                                                                                 
tttttGCCAGCGTAGCGTACTGCATTCTGTGCGTG‑GAAGTGCCATACGTTGAAACGGCACCGTAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                                                      <  1:383909/143‑4 (MQ=255)
tttttGCCAGCGTACCGTACTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                                                      <  1:507229/143‑4 (MQ=255)
tttttGCCAGCGTACCGTACTGCATTCAGGTGCTG‑GAAGTGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                                                      <  2:239007/143‑4 (MQ=255)
 ttttGCTAGCTTTCCGTACTGTATT‑ACTTTCTGAAAAGTTCTATACGTTGAAACGGCAACGGAAGTGGGTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttttt                                                                                                                     <  1:517003/143‑5 (MQ=255)
 ttttGCTAGCGTATTGTACTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCCAGACGTGCCGGTTGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                                                     <  1:525339/143‑1 (MQ=255)
 tgttGGCAGCGTGTCGTGCGGCGTTCAGGTGCGG‑GAAGAGGTATACGTTGAAACGGGACCGGAAGTGGTTGCAGCATTCAGAGCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                                                     <  2:756897/141‑1 (MQ=255)
  tttGCCAGCGTACCGTACTGCATTCAGTTGCAT‑GAAGGCCCATACGTTGAAACTGCACCGGAAGTGATTGCAGCATTCAGACCGAAGAGTGAAAGACGCGACGTTAGCGAATAAACAAAAGCCCCCGAACGGGGGGAACTCaa                                                                                                                    >  2:372612/1‑143 (MQ=255)
   ttGCCAGCGTACCGTACTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaa                                                                                                                   <  1:652897/143‑1 (MQ=255)
     gCCAGCGTACCGTACTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaaa                                                                                                                 >  1:373759/1‑142 (MQ=255)
                gtacTGCATTCAGTTGCTG‑GAAGTGCCATACGTTGTAACGGTACCGGAAGTGGTTGCAGCATTCAGACCAAAGAGTGCAAGACGCGACGTGAGCGAAGAAAAAAAGCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                                                      <  1:114003/143‑1 (MQ=255)
                gtacTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                                                      <  2:56924/143‑1 (MQ=255)
                 tccttcATTGAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCTAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                                                     <  2:62193/138‑1 (MQ=255)
                 tagtGCATTCGGTTTGTG‑GATGTGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCATACCGAAGCGTGCAAGACGGGACATTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACCATTAGTGGGAt                                                                                                     <  2:436676/140‑1 (MQ=255)
                     gCATTCAGTTGCTG‑GAAGAGCCATACGTTGAAACGGCAGCGGAAGTGGTTGCAGCATTCAGACCGCAGCGAGCAAGACGCGACGTTAGCGACAAAAAAAATCCCCACGAGCGGAGGGATCTCAAGACAATTAGTggggatcga                                                                                                 >  2:597206/1‑137 (MQ=255)
                      cATTCAGTTGCTG‑GAAGGGACGACGGGGGAAACGGGACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAACAATCCCCCCGAGCAGGGGGATCTCAAAACAAATAGCGAGACTCaac                                                                                                >  2:464994/1‑141 (MQ=255)
                         tCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                                                             >  1:316612/1‑143 (MQ=255)
                            gTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATCAAAACATCCCCCCGCGCGGGGGGATCTCAAAACAATAAGTGGGATTCACCAATCga                                                                                          >  1:592214/1‑142 (MQ=255)
                                    gAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACTCGTGGTTACCGAATAAAAAAATCAACCAGAGAGGGGGGATATAAAAACAATTAGTGGGATTAAACAATCAGCAGAACg                                                                                   >  1:436891/1‑143 (MQ=255)
                                           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGCCGCGCCGGTAGCGAATCAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACACTCAGTGGCCTTCACCACTCGGCCGCACGACGCGAc                                                                            >  2:537522/1‑143 (MQ=255)
                                            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGAGGAGCTCAAAACAATCAGTGGGATTCACCAAGCGGCAGAACGGTGGGAcc                                                                           >  2:540696/1‑143 (MQ=255)
                                                gTTGAAACGGCACCGGAAGTGGTTGCAGCAACCTTTTGGTAGTTTGCATGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                                                       <  2:155643/143‑1 (MQ=255)
                                                     cgcGGCACCGGAAGTGGGTGCAGCAGTCAGGCCGCAGAGTGAAAGACGCGACGGTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                                                  <  2:121016/141‑1 (MQ=255)
                                                     aaCGGCACCGGATGTTTTTTCAGCATTCAGACCGAAGAGTGCAAGGCGGGACGGGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTc                                                                  <  1:65718/143‑1 (MQ=255)
                                                             cGGAAGGGGGTGCAGCAGTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAAGCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                                                          <  1:22977/143‑1 (MQ=255)
                                                               gAAGTGGTTGCAGCATGCAGGCCGAAGTGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATCAGTGGGATTCACCAATCGACAGTACAGTGGGCCCAAACTGCTCGTTCAGTACt                                                        <  2:464263/143‑1 (MQ=255)
                                                                     gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGCCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGCTCTCAAAACAATTAGTGGGATTCACCTGTTTGCGGGACGCAGCCACCAAACTGCTCGTTCAGTACTTCAccc                                                  >  1:137903/1‑143 (MQ=255)
                                                                      ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCGCCGAGCGGGGGCACCTCCAAACAATCAGTGGGACTCACCAATCGGAAGCACGGAGCGACCAAACCGCTCGTCCAGTACGTCACCCa                                                 >  2:709536/1‑143 (MQ=255)
                                                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGGTTCACCAATCGGCAGGACGGGGAGGCCAAACTGCTCGTTCAGTACCTCTCCCCTCg                                              >  1:316628/1‑143 (MQ=255)
                                                                         aagcaTTCAGACAGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAACACCCGCCCGAGCGGGGGGATCGCAAAACGATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                                              <  1:229575/142‑1 (MQ=255)
                                                                                cAGACCGCGGACTGCAAGACGCGACGTGAGCGAATAACGAAGTGCCCCCGAGCGGGGGGATCGCAAAACACGTAGTGGGAGTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata                                       <  1:537522/143‑1 (MQ=255)
                                                                                  caCCGGAGAGTGGAGGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACATCTGCTGGGATTCACCAATCGGCAGAACGGTGCGACCAATCTGCTCGTTCATTACTTCACCCATCGCCagataga                                     <  2:623967/142‑1 (MQ=255)
                                                                                   aCCGAAGAGTGCAAGACGCGACGATAGCGAATAAAAAAATCCCCCCGAGAGGGAGGATCTCACAACAATTAGTGAGATTCACCAATCGGCAGCACGGTGCGACCAAACTGCTCGTGCAGCACTACACCCAAAGACagatagat                                    >  2:174902/1‑143 (MQ=255)
                                                                                       aaGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACCGCTCCCCCCGCACTTCACCGATCGCCagacagatcccg                                >  1:313995/1‑139 (MQ=255)
                                                                                           gTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTgg                            <  2:316612/143‑1 (MQ=255)
                                                                                               aaGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCAcc                        >  1:692254/1‑143 (MQ=255)
                                                                                                        cGTTACTGAATAAAAGAATGCCCCCGAGCGGCGGGATCTCAAAACAATGAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAg               <  2:651956/143‑1 (MQ=255)
                                                                                                              cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGTCACCGCAGATCAGCCAAAt         >  1:281239/1‑143 (MQ=255)
                                                                                                               gAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGAACCGCAGATCAGCCAAATc        >  1:700634/1‑143 (MQ=255)
                                                                                                                  ttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGATCGCTCAGTACTTCACCCATCGCCAGATAGATTGCGCCGGCACCGCAGATCCGCCAAATCCAg     >  2:303521/8‑143 (MQ=255)
                                                                                                                     aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCACAGCTCAGCCCAATCCAGCCg  >  1:119157/1‑143 (MQ=255)
                                                                                                                     aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGCTCAGTACTTCACCCATCGCCAGATAGATTGCGCAGGCACAGCAGATAAGCCAAATCCAGCCg  >  1:92742/1‑143 (MQ=255)
                                                                                                                                  |                                                                                                                                 
TTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTG‑GAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCG  >  NC_000913/9778‑10036

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