Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,293,960 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate‑aspartate symport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,293,9591.C94.4% 18.6 / ‑2.8 18intergenic (+586/+56)gltP/yjcOglutamate‑aspartate symport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (5/12);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
*NC_0009134,293,9592.G93.8% 16.4 / ‑6.2 17intergenic (+586/+56)gltP/yjcOglutamate‑aspartate symport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (5/10);  minor base C (0/1);  total (5/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGG  >  NC_000913/4293887‑4294025
                                                                           ||                                                                  
gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑ccc                                                                 >  2:2507496/1‑75 (MQ=14)
                                                       ggTCAGAAAAGGCGTTT‑‑ACCCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCt            <  1:2151494/76‑1 (MQ=1)
                                                          cGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa         <  1:1404558/76‑1 (MQ=1)
                                                           ggATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc        >  1:1687161/1‑76 (MQ=11)
                                                           ggATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc        <  2:2316062/76‑1 (MQ=1)
                                                           ggATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCATCTTATCAGGCCTATCTTAAc        >  1:5742938/1‑76 (MQ=11)
                                                            gATAAGGCGTTT‑‑ACGCCGCATCCGACATTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc       >  1:3368722/1‑76 (MQ=11)
                                                            gATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc       <  2:4632592/76‑1 (MQ=11)
                                                             aTAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg      <  2:4778750/76‑1 (MQ=11)
                                                               agggCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt    <  2:4871635/74‑1 (MQ=11)
                                                               aaGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt     <  1:491361/75‑1 (MQ=11)
                                                               aaGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt     >  2:491361/1‑75 (MQ=11)
                                                                aGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg   >  2:2836345/1‑76 (MQ=14)
                                                                aGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg   <  2:5216430/76‑1 (MQ=11)
                                                                aGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg   <  2:5263377/76‑1 (MQ=11)
                                                                 ggCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg  <  1:4744569/76‑1 (MQ=14)
                                                                 ggCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg  <  1:282287/76‑1 (MQ=14)
                                                                 ggCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg  <  2:4917137/76‑1 (MQ=14)
                                                                 ggCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg  <  2:893000/76‑1 (MQ=14)
                                                                           ||                                                                  
GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGG  >  NC_000913/4293887‑4294025

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