Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 765,132 G→A 23.8% intergenic (+83/+21) 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,1320GA23.8% 34.1 / 10.3 21intergenic (+83/+21)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 G (8/8);  new base A (3/2);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-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)
                   aaGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTa                                                                                                                                              >  3:855440/1‑74 (MQ=255)
                   aaGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTa                                                                                                                                              <  4:855440/74‑1 (MQ=255)
                                      ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTtgtg                                                                                                                             >  3:88230/1‑72 (MQ=255)
                                      ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTtgtg                                                                                                                             <  4:88230/72‑1 (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)
                                                                                          |                                                                                                                                              
TGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAA  >  NC_000913/765042‑765274

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