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.C90.0% 76.1 / 3.0 30intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (15/12);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_pgi4,295,6402.G90.0% 73.2 / 2.7 30intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (15/12);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_pgi/4295557‑4295777
                                                                                      ||                                                                                                                                            
aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGt                                                                            >  1:392719/1‑149 (MQ=2)
        cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                    <  2:615179/149‑1 (MQ=255)
          cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                  >  1:284819/1‑149 (MQ=255)
           gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGGTAAtttt                                                                 >  2:962599/1‑149 (MQ=14)
               tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAGTTTTTCaa                                                             >  1:547454/1‑149 (MQ=255)
                 tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCCACGCTTGCCGCCTCTTATCCGGCCTATCTTCACCGTTGTTTAATTTTTCacac                                                           >  2:243855/1‑149 (MQ=11)
                    ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                        >  2:1058226/1‑149 (MQ=255)
                       cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGGCGCGTCGCATCAGGCCTCTCTTACCCGCTGGTTAATTTTCCAACCTCTGca                                                     >  1:419326/1‑149 (MQ=11)
                            ccAGACAGCGGAATAACCTCATTTAGCGTGATTTTATAGGTCGGATAAGGCGGTT‑‑ACGC‑‑CG‑CATCCGACATCAATTCCTCATGCGACGCTTTCCGCGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                <  2:21008/149‑1 (MQ=17)
                              aGACATCGCGATAGCCTGATTTATCGTGATTTTGTAGATCGGATAAGTCGTTT‑‑A‑‑CTGCG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                              <  2:1033710/149‑1 (MQ=18)
                              aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                              <  2:558767/149‑1 (MQ=255)
                                     gcAATAGCCTTATGTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                       <  2:256401/149‑1 (MQ=18)
                                     gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                       <  2:392719/149‑1 (MQ=255)
                                         tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                             >  1:541598/1‑139 (MQ=255)
                                         tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                             <  2:541598/139‑1 (MQ=255)
                                             cTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               <  1:363140/149‑1 (MQ=255)
                                                aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                            <  2:253097/149‑1 (MQ=255)
                                                aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAACACCCGGTATCAAACCCTTCCATAc                            >  1:253097/1‑149 (MQ=255)
                                                 tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                             >  1:189104/1‑147 (MQ=255)
                                                 tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                             <  2:189104/147‑1 (MQ=255)
                                                      cGTGATTTTGTAGGTCGGGTAAGGCCTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                      <  1:728723/149‑1 (MQ=25)
                                                      cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                      >  1:895284/1‑149 (MQ=255)
                                                         gATTTTGTAGGTGGGATAAGGCGTTT‑‑A‑‑CG‑CGCCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                   <  2:640836/149‑1 (MQ=25)
                                                         gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACCACCGGTCTCAAACCCTTCCTTCCAGCTCCGGt                   >  1:76381/1‑149 (MQ=18)
                                                             ttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTATCTTAAACGTTGGTTAATTTTTCTAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa               >  2:62365/1‑149 (MQ=25)
                                                                  ggTCGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCATACCCTTGCATACAGTTTAGGTTCAACCAGt          >  1:416800/1‑149 (MQ=25)
                                                                     cGGATAAGGCGTTTCCG‑‑C‑‑CG‑CATGCGACATCCATGCCTGTTTCGGCGCTGGCAGTGGCCTATCACGCCTATCTTATGCGTGGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa       <  1:747387/149‑1 (MQ=255)
                                                                     cGGATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCTTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa       <  1:746430/149‑1 (MQ=255)
                                                                      ggATAAGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTCCTGTTTCTGAGCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg      <  1:206181/149‑1 (MQ=37)
                                                                          aaGGCGTTT‑‑ACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc  >  1:592100/1‑149 (MQ=255)
                                                                                      ||                                                                                                                                            
ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_pgi/4295557‑4295777

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: ‑