Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,295,641 +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_pgi4,295,6401.C93.3% 41.9 / ‑2.4 15intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (8/6);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
*NC_000913_3_pae_pgi4,295,6402.G93.3% 41.0 / ‑2.3 15intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (8/6);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913_3_pae_pgi/4295548‑4295779
                                                                                             ||                                                                                                                                              
ccTGATGCGACGCTTACCGCGTCTTATCAGGTCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCGGATAAGCCGCTCA‑‑CCGGACAGCTTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAt                                                                                        <  1:388397/149‑1 (MQ=11)
            cTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGTCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                            <  1:148831/149‑1 (MQ=11)
            cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                            >  2:116006/1‑149 (MQ=255)
               aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                         >  1:224515/1‑149 (MQ=255)
                      cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTc                                                                  >  1:74958/1‑149 (MQ=255)
                              gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         >  1:404641/1‑149 (MQ=255)
                                 tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                      <  1:116006/149‑1 (MQ=255)
                                                     ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  <  1:75510/149‑1 (MQ=255)
                                                      cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAATCCTTCCa                                 >  2:228304/1‑149 (MQ=255)
                                                          tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTATATACCGTCTTACACCCGGTATCAAACCCTTCCATACa                             <  1:306622/149‑1 (MQ=17)
                                                           ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                            <  2:333182/149‑1 (MQ=255)
                                                                gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       >  2:339751/1‑149 (MQ=255)
                                                                        gTAGGTCGGATAAGGCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGTTTCAAcc               >  1:206648/1‑149 (MQ=255)
                                                                                     ggCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCATTCCACCCCCtt  >  2:41770/1‑149 (MQ=25)
                                                                                     ggCGTTTACGCCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt  <  1:114910/149‑1 (MQ=255)
                                                                                             ||                                                                                                                                              
CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCG‑‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913_3_pae_pgi/4295548‑4295779

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