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‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_pgi/4295557‑4295777
                                                                                      ||                                                                                                                                            
aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGt                                                                            >  1:392719/1‑149 (MQ=2)
        cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                    <  2:615179/149‑1 (MQ=255)
          cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                  >  1:284819/1‑149 (MQ=255)
           gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGGTAAtttt                                                                 >  2:962599/1‑149 (MQ=14)
               tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAGTTTTTCaa                                                             >  1:547454/1‑149 (MQ=255)
                 tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCCACGCTTGCCGCCTCTTATCCGGCCTATCTTCACCGTTGTTTAATTTTTCacac                                                           >  2:243855/1‑149 (MQ=11)
                    ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                        >  2:1058226/1‑149 (MQ=255)
                       cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGGCGCGTCGCATCAGGCCTCTCTTACCCGCTGGTTAATTTTCCAACCTCTGca                                                     >  1:419326/1‑149 (MQ=11)
                            ccAGACAGCGGAATAACCTCATTTAGCGTGATTTTATAGGTCGGATAAGGCGGTT‑‑ACGC‑‑CGCATCCGACATCAATTCCTCATGCGACGCTTTCCGCGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                <  2:21008/149‑1 (MQ=17)
                              aGACATCGCGATAGCCTGATTTATCGTGATTTTGTAGATCGGATAAGTCGTTT‑‑A‑‑CTGCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                              <  2:1033710/149‑1 (MQ=18)
                              aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                              <  2:558767/149‑1 (MQ=255)
                                     gcAATAGCCTTATGTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                       <  2:256401/149‑1 (MQ=18)
                                     gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                       <  2:392719/149‑1 (MQ=255)
                                         tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                             >  1:541598/1‑139 (MQ=255)
                                         tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                             <  2:541598/139‑1 (MQ=255)
                                             cTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               <  1:363140/149‑1 (MQ=255)
                                                aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                            <  2:253097/149‑1 (MQ=255)
                                                aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAACACCCGGTATCAAACCCTTCCATAc                            >  1:253097/1‑149 (MQ=255)
                                                 tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                             >  1:189104/1‑147 (MQ=255)
                                                 tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                             <  2:189104/147‑1 (MQ=255)
                                                      cGTGATTTTGTAGGTCGGGTAAGGCCTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                      <  1:728723/149‑1 (MQ=25)
                                                      cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                      >  1:895284/1‑149 (MQ=255)
                                                         gATTTTGTAGGTGGGATAAGGCGTTT‑‑A‑‑CGCGCCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                   <  2:640836/149‑1 (MQ=25)
                                                         gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACCACCGGTCTCAAACCCTTCCTTCCAGCTCCGGt                   >  1:76381/1‑149 (MQ=18)
                                                             ttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTATCTTAAACGTTGGTTAATTTTTCTAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa               >  2:62365/1‑149 (MQ=25)
                                                                  ggTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCATACCCTTGCATACAGTTTAGGTTCAACCAGt          >  1:416800/1‑149 (MQ=25)
                                                                     cGGATAAGGCGTTTCCG‑‑C‑‑CGCATGCGACATCCATGCCTGTTTCGGCGCTGGCAGTGGCCTATCACGCCTATCTTATGCGTGGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa       <  1:747387/149‑1 (MQ=255)
                                                                     cGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCTTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa       <  1:746430/149‑1 (MQ=255)
                                                                      ggATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTCCTGTTTCTGAGCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg      <  1:206181/149‑1 (MQ=37)
                                                                          aaGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc  >  1:592100/1‑149 (MQ=255)
                                                                                      ||                                                                                                                                            
ACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  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: