Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C92.0% 61.0 / 1.3 25intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (12/11);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01
*NC_0009134,296,3802.G92.0% 61.1 / 1.4 25intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (12/11);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGC  >  NC_000913/4296284‑4296511
                                                                                                   ||                                                                                                                                     
atcggcTGTTGCGCCGCTGAACGCGTCGTATCAGGCCTACGCGATTGTGTGTTATTGCCTGTTGTAGAGGGTTTTTGTAGTTAGGATAATTCGGTTTTG‑‑C‑‑CGCATCGGAGATCTTTGGCTGTTGCGACGCTTGCCGCGTCTTATCAGGc                                                                                   <  4:16604/144‑1 (MQ=2)
            gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGGTAGCGTGAGTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCCTCCGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTGTCTTAACCg                                                                       >  4:111845/1‑149 (MQ=1)
             aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                      >  6:86110/1‑149 (MQ=12)
              cGCTTAACGCGTCTTATCAGGCCTTCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATATGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                     <  1:38426/149‑1 (MQ=11)
              cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                     <  2:96677/149‑1 (MQ=12)
                   aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                 >  5:14003/1‑148 (MQ=255)
                    aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                               >  8:22760/1‑149 (MQ=255)
                             aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                      >  2:105792/1‑149 (MQ=255)
                              tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                     >  2:137254/1‑149 (MQ=255)
                                   acTACTCCACACAGCGCAATAGCCTGATTTAGCGTGATGAGGTACGTCGGATAAGGCG‑‑CTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                <  3:83468/148‑1 (MQ=17)
                                       cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                              <  7:22760/147‑1 (MQ=255)
                                          cagacagaGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                         <  6:121578/149‑1 (MQ=255)
                                            gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                       <  2:106748/149‑1 (MQ=255)
                                                           tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                 <  4:84703/124‑1 (MQ=255)
                                                           tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                 >  3:84703/1‑124 (MQ=255)
                                                            gATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                         >  7:58437/1‑99 (MQ=1)
                                                            gATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                         <  8:58437/99‑1 (MQ=1)
                                                             aTTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                       >  5:58277/1‑148 (MQ=255)
                                                                   cGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                  >  7:99894/1‑147 (MQ=255)
                                                                            gTAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc        <  5:86110/148‑1 (MQ=255)
                                                                            gTAGGTCGGATAAGGCG‑‑TTTA‑‑CGACTCATCCGACATCAATTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc        <  7:63737/148‑1 (MQ=37)
                                                                             tAGGTCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc        >  4:86587/1‑147 (MQ=255)
                                                                                 tCGGATAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc  <  1:43776/149‑1 (MQ=255)
                                                                                     aTAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                <  1:32298/131‑1 (MQ=255)
                                                                                     aTAAGGCG‑‑TTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                >  2:32298/1‑131 (MQ=255)
                                                                                                   ||                                                                                                                                     
AACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑‑TTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGC  >  NC_000913/4296284‑4296511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: