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.C93.3% 82.0 / ‑0.7 30intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (10/18);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_tpiA4,296,2962.G93.3% 79.0 / ‑0.6 30intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (10/18);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_pae_tpiA/4296212‑4296430
                                                                                     ||                                                                                                                                         
gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                          >  1:472395/1‑149 (MQ=9)
   gCTTAGCGCGTTTTATCAGGCCGACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCTGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                       <  2:199611/149‑1 (MQ=11)
              ctgataaGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGTATAAGGCGTTTACGC‑C‑GCATCAGATATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                            <  2:423575/143‑1 (MQ=14)
                 aTCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                         <  1:490636/149‑1 (MQ=255)
                 aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTGGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                         >  1:331171/1‑149 (MQ=255)
                    aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                      >  1:298560/1‑149 (MQ=255)
                    aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                      >  1:321638/1‑149 (MQ=255)
                        cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  <  1:160329/149‑1 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                              >  1:256764/1‑149 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    <  2:7584/141‑1 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    >  1:7584/1‑141 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                            >  1:349633/1‑149 (MQ=255)
                                    gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                      <  1:15397/149‑1 (MQ=255)
                                      gccaTAGCCTGATTTAGCGGGATTTTGTAGGGCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTGGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                    <  1:43807/146‑1 (MQ=18)
                                      gcAATAGCCTGATTTAGCGTGATTTTGCAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                    <  1:383188/149‑1 (MQ=255)
                                          tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAACGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCt                                <  2:288848/149‑1 (MQ=255)
                                           aGCCTGATTGAGGGTGATTTTTTAGATCGGATAAGGCTTTTACGC‑C‑GCATTAGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                               <  2:264837/149‑1 (MQ=18)
                                              cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                            <  1:214107/149‑1 (MQ=255)
                                                      gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                    <  2:269478/149‑1 (MQ=255)
                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                  <  1:207160/149‑1 (MQ=255)
                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                  >  2:492389/1‑149 (MQ=255)
                                                        gTGATTTTGTAGGTCGCATAAGCCGTTAA‑‑CGCAGCATCCGACATGAATGCC‑TGACGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                  <  1:318605/149‑1 (MQ=25)
                                                          gATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                <  1:492389/149‑1 (MQ=255)
                                                              ttGTAGGTCGGATAAGGCGTTTA‑‑CGCAGCACATGAGCACATTACCTTGATG‑GACGCTTTCCCCGTCGTAGCAGACCTATCTTAACCTTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa            <  2:418736/149‑1 (MQ=255)
                                                                   ggTCGGATAAGTCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt       <  2:414241/149‑1 (MQ=255)
                                                                        gATAAGGCGTTTACGC‑C‑GCATCCGACATCACTCGA‑AGCGTCTATGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGagcagc  >  1:493160/1‑149 (MQ=37)
                                                                        gATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc  <  1:159732/149‑1 (MQ=255)
                                                                        gATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc  <  2:256764/149‑1 (MQ=255)
                                                                           aaGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              >  1:98169/1‑118 (MQ=255)
                                                                           aaGGCGTTTACGC‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              <  2:98169/118‑1 (MQ=255)
                                                                                     ||                                                                                                                                         
GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCC‑TGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_pae_tpiA/4296212‑4296430

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