Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,380:1 +CG 100% intergenic (+586/+56) gltP → / ← yjcO glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C92.0% 66.3 / ‑0.7 25intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (8/15);  total (9/16)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0009134,296,3802.G92.0% 64.4 / ‑0.7 25intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (8/15);  total (9/16)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGC  >  NC_000913/4296285‑4296511
                                                                                                ||                                                                                                                                         
aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTATA‑‑CGACAAGTGATTTA‑‑‑TCAAGGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         >  2:454865/1‑143 (MQ=11)
                  gaCGCGTCTGATCGGGCCTACGCCGGGCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGt                                                                      <  2:640566/142‑1 (MQ=11)
                   aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGtt                                                                     <  1:24535/143‑1 (MQ=255)
                           tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAAtt                                                             <  1:89428/143‑1 (MQ=255)
                              cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAAttttt                                                          <  1:454865/143‑1 (MQ=255)
                               aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTc                                                         <  1:181388/143‑1 (MQ=255)
                                ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCa                                                        <  1:164143/143‑1 (MQ=255)
                                 gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTATTA‑ACCGGTGGTTGATTTTTcga                                                       >  1:637417/1‑141 (MQ=11)
                                             gaccgcAATAGCCTGATTTAGCGTGATTTTGGAGTTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                           <  2:528875/140‑1 (MQ=17)
                                              aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                          >  1:252824/1‑143 (MQ=255)
                                                cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                        >  1:545976/1‑143 (MQ=255)
                                                      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTA‑ACCGCTGGTTACCTTTTCAAACTCTTCACACCCGGCATCaa                                  >  2:4333/1‑143 (MQ=17)
                                                      aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGCTGGTTAATTTTTCAAACTCTTCCCACCCGGGATCac                                  >  1:351226/1‑142 (MQ=17)
                                                       gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                 <  1:328606/143‑1 (MQ=255)
                                                             tttAGCGTGATTTTGTAGTTCGGATAAGGCGTTTA‑‑‑CGCAGCTGA‑‑CAACTTCAAGTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTTTTT‑ACCGCCGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                           <  2:594649/143‑1 (MQ=17)
                                                             tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                           <  2:3818/143‑1 (MQ=255)
                                                              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                          <  1:413342/143‑1 (MQ=255)
                                                              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                          <  1:219026/143‑1 (MQ=255)
                                                                    tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                    <  1:396862/143‑1 (MQ=255)
                                                                               gTCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa         <  2:501719/143‑1 (MQ=255)
                                                                                tCGGATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCATACAAAGTTGG‑TATGTTGTCGAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc        >  1:471096/1‑143 (MQ=17)
                                                                                  ggATAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa      <  1:249124/143‑1 (MQ=255)
                                                                                    aTAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAAGCCTTCCATCCAGCTCAGGTTCatacaat    >  1:505804/1‑137 (MQ=25)
                                                                                    aTAAGGCGTTTACGCCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt    <  2:10072/143‑1 (MQ=255)
                                                                                      aaGGCGTTTACGCCGCATCCAA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGACCTATCTTA‑ACCGTTGGTTAATCTTTCAAACACTTCACTCCCGGTATCACACCCCTCCATACAGCTCAGGATCCACCAGtcc  >  2:639316/1‑141 (MQ=21)
                                                                                                ||                                                                                                                                         
ACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGA‑‑CA‑‑‑TCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTA‑ACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGC  >  NC_000913/4296285‑4296511

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