Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 4,296,411 +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_000913_3_vch_tpiA4,296,4101.C95.2% 54.9 / ‑3.3 21intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (9/11);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
*NC_000913_3_vch_tpiA4,296,4102.G95.2% 57.9 / ‑3.3 21intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (9/11);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_vch_tpiA/4296327‑4296548
                                                                                    ||                                                                                                                                            
aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTATCTTAACCGt                                                                             >  1:304727/1‑149 (MQ=255)
              ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCa                                                               >  2:252008/1‑149 (MQ=255)
                aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                             <  1:385790/149‑1 (MQ=255)
                 tgaGGCCTGCGGCGGGCAGGGCAATAGCATGATTTAGCGTGGTTTTGTAGGTGGGGTAAGGGGTTTA‑‑CGCGGGCTCCGACATCAATGACAGATGAGACGCTTGCCGCGTCTTATCAGGACTGTCTTAACCGTTGGTTGATTTTTCAAAc                                                            >  2:43958/3‑149 (MQ=9)
                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                  >  1:136646/1‑149 (MQ=255)
                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                  >  1:411901/1‑149 (MQ=255)
                             cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                >  1:325991/1‑149 (MQ=255)
                              aGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATACGGCGATTACGC‑C‑GCAACCGACATCGATGCCTGATGCGACGCTTGCTGCGTCTTATCAGGCCTATCTTAACCGTTGCTTCATTTTTCAAACTCTTCCCACCCgg                                               >  2:391804/1‑149 (MQ=17)
                              aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                               >  2:340091/1‑149 (MQ=255)
                                     gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                        <  1:340091/149‑1 (MQ=255)
                                     gcAATAGCCTGATTGAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                        <  2:195491/149‑1 (MQ=18)
                                       aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                      <  2:136646/149‑1 (MQ=255)
                                           gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  1:78140/149‑1 (MQ=255)
                                             cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GACTACGGCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                <  1:290832/149‑1 (MQ=21)
                                                 tttAGGGTGATTTTGTAGGTCGGATAAGGCGTTCACGC‑C‑CCCTCTGACATTAATGCCTGATGCGCTGCATCATGCCTCCAATCACCGGTATCTGAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                            <  1:369688/149‑1 (MQ=16)
                                                     gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCTTTCTCTTTACATCCGGTATCTAACCCTTCCCTACAGCTc                        >  2:419271/1‑149 (MQ=18)
                                                      cGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                       <  2:164002/149‑1 (MQ=255)
                                                       gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                      <  1:334948/149‑1 (MQ=255)
                                                              tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc               <  1:291918/149‑1 (MQ=255)
                                                                  ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGCTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt           <  1:439113/149‑1 (MQ=255)
                                                                           aGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt  >  2:196687/1‑149 (MQ=255)
                                                                                    ||                                                                                                                                            
ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_vch_tpiA/4296327‑4296548

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