Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 765,135 T→C 29.4% intergenic (+86/+18) sucD → / ← mngR succinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit/transcriptional repressor for the mannosyl‑D‑glycerate catabolic operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913765,1350TC29.4% 21.3 / 10.3 17intergenic (+86/+18)sucD/mngRsuccinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit/transcriptional repressor for the mannosyl‑D‑glycerate catabolic operon
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base C (3/2);  total (9/8)
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.05e-01

TGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAA  >  NC_000913/765042‑765274
                                                                                             |                                                                                                                                           
tGAAATAAATCTCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGTAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTatgat                                                                                                                                          >  3:374053/1‑97 (MQ=255)
tGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGTCTGTCGGTTTCGACATGGTTGGCCATCGTatgat                                                                                                                                          <  4:374053/97‑1 (MQ=255)
     taaataTCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTatga                                                                                                                                           >  4:1059151/1‑91 (MQ=255)
     taaataTCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTatga                                                                                                                                           <  3:1059151/91‑1 (MQ=255)
                   aaGAAATAGCCCTCGCCGCTTCCCTCTACAGGAGTGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCATACGATGGc                                                                                                                                       >  1:145405/1‑81 (MQ=255)
                   aaGAAATAGCCCTCGCCGCTTCCCTCTACAGGAGTGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCATACGATGGc                                                                                                                                       <  2:145405/81‑1 (MQ=255)
                                      ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTtgtg                                                                                                                             <  4:88230/72‑1 (MQ=255)
                                      ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTtgtg                                                                                                                             >  3:88230/1‑72 (MQ=255)
                                                                 gTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGAc                                      >  3:131546/1‑132 (MQ=255)
                                                                 gTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGAc                                      <  4:131546/132‑1 (MQ=255)
                                                                    ggTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTtat                                                                                   <  2:586759/84‑1 (MQ=255)
                                                                    ggTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTtat                                                                                   >  1:586759/1‑84 (MQ=255)
                                                                      tttCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTgcgc                                                                                                    <  3:10902/65‑1 (MQ=255)
                                                                      tttCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTgcgc                                                                                                    >  4:10902/1‑65 (MQ=255)
                                                                                     ccATCATACGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGc                                                        <  2:113980/94‑1 (MQ=255)
                                                                                     ccATCATACGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGc                                                        >  1:113980/1‑94 (MQ=255)
                                                                                     ccATCATACGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGaa  >  1:908163/1‑148 (MQ=255)
                                                                                      cATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCAcg                   <  3:703307/130‑1 (MQ=255)
                                                                                      cATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCAcg                   >  4:703307/1‑130 (MQ=255)
                                                                                           tATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAat                          >  3:34927/1‑118 (MQ=255)
                                                                                           tATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAat                          <  4:34927/118‑1 (MQ=255)
                                                                                             |                                                                                                                                           
TGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAA  >  NC_000913/765042‑765274

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