Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,296,297 +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_pae_tpiA4,296,2961.C92.3% 23.2 / ‑3.2 14intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (3/9);  minor base G (1/0);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01
*NC_000913_3_pae_tpiA4,296,2962.G92.9% 27.1 / ‑2.0 14intergenic (+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 (4/9);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_tpiA/4296212‑4296433
                                                                                     ||                                                                                                                                           
gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCCTATCGGGCCTTTCTTAACCg                                                                             >  1:412315/1‑149 (MQ=1)
          gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                   <  1:137113/149‑1 (MQ=255)
           cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                  <  1:510823/149‑1 (MQ=255)
                        cTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                     <  1:131367/149‑1 (MQ=255)
                         tACGCCAGACGGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGGCGCATCCGACATCAATGAATTAGGCGGTGTTCGACTCGTTATATCAGGCCCAGCTTAACCGTTGGTTAATTTTTCAAACTCttacc                                                    >  2:417438/1‑146 (MQ=9)
                                         aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGC‑‑CGCATCCGACGTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaacacc                                    >  2:87753/1‑145 (MQ=18)
                                            gCCTGATTTTGCGTGATTTTGTAGGTCGGATAAGGCCTTTACGC‑‑CGCATTCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                 <  2:63867/149‑1 (MQ=21)
                                             ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                   <  1:389171/114‑1 (MQ=255)
                                             ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                   >  2:389171/1‑114 (MQ=255)
                                                      gggTGAATTTGTAAGGCGTACAAGGCGTATACGC‑‑CCCACCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                       <  1:13870/147‑1 (MQ=21)
                                                                     tCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                     <  1:415678/136‑1 (MQ=255)
                                                                     tCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                     >  2:415678/1‑136 (MQ=255)
                                                                       gggtAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGGTGCGACGCTTGCCGCGTCTTATCAGGCTTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg      <  1:437493/146‑1 (MQ=37)
                                                                           aaGGCGTTTACGC‑‑CGCATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc  <  1:10449/149‑1 (MQ=255)
                                                                                     ||                                                                                                                                           
GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_tpiA/4296212‑4296433

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