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.C90.9% 46.2 / ‑0.6 22intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (9/10);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
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.G90.9% 45.8 / ‑0.7 22intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (9/11);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAG  >  NC_000913/4296301‑4296513
                                                                                  ||                                                                                                                                       
ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                          >  1:793138/1‑143 (MQ=2)
 tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                         >  1:893820/1‑143 (MQ=9)
   aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATGAGGCCTATCTTaaacctg                                                                       >  1:73036/1‑138 (MQ=11)
               tggCGTACGCCAGAGAGGGCACTAGCCTGATTTGGCGTGATTTTGTAGGCCGGATAGGGCGAT‑‑TA‑‑CGAAGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCCTCTGATCTGGCTTATCTTAACCGTTGGTTAATTTTTc                                                           <  2:697100/142‑1 (MQ=21)
               tggCCGACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                           <  1:785675/142‑1 (MQ=14)
                ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                          <  1:793662/143‑1 (MQ=255)
                  ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGTG‑‑G‑‑CGGCTCCGACACCAATGCCTGGGTCCACGCGTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                        <  1:914817/143‑1 (MQ=17)
                              aGCGGAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                            <  1:805391/143‑1 (MQ=17)
                                     tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                     >  1:647848/1‑143 (MQ=255)
                                       gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                   <  1:85549/143‑1 (MQ=255)
                                       gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                   <  1:837500/143‑1 (MQ=255)
                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                  >  2:85048/1‑143 (MQ=255)
                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                <  2:639966/143‑1 (MQ=255)
                                             tttCGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                             <  1:832181/143‑1 (MQ=255)
                                             ttcaGCCTGATTTCGTATGTCCGTTAATGCGTT‑‑TATGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCGTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                             <  1:484838/140‑1 (MQ=21)
                                              ttAGCGTGATTGTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                            <  1:80116/143‑1 (MQ=255)
                                                 gCGGGATTTGGTAGGTCGGATGAGGCGTC‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGGTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                         <  2:317909/143‑1 (MQ=21)
                                                  cGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCCTGCCGCGTCTGATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                        >  2:407732/1‑143 (MQ=21)
                                                           gTAGGTTGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt               <  1:194005/143‑1 (MQ=255)
                                                              ggTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCGCACTCTTCACACGCGGTATCAAACCCTTCCATACAGCTCAGGtcta            >  1:259405/1‑140 (MQ=21)
                                                               gTCGGATAAGGCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa           >  1:264855/1‑143 (MQ=255)
                                                                        ggCGTT‑‑TACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCTTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTCCAg  >  1:65837/1‑143 (MQ=37)
                                                                                  ||                                                                                                                                       
TTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAG  >  NC_000913/4296301‑4296513

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: