Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 4,296,381:1 +CG 100% intergenic (+586/+56) 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_bme_pgi4,296,3811.C93.2% 94.5 / 2.1 44intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (12/29);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_000913_3_bme_pgi4,296,3812.G93.2% 92.5 / 2.1 44intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (12/29);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_bme_pgi/4296296‑4296515
                                                                                        ||                                                                                                                                        
cGACGCCTGCCGCGTCTTATGAGGCCGAAGCCAGGCAGCGCAATGGCCAGATTGCGGGTGATTTTTTAGGTCGGGCCAGGGGTTT‑‑ACGC‑C‑GCATACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                           <  2:369265/149‑1 (MQ=2)
  aCGCTTAACGCGTCTTATCAGGCCTACCCCAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTT‑‑A‑‑CGCAGAATCCGACATCACTCCCCGATGCGACGCTTCCCGCGTCCTCTCCGGCCTCTCTTAAcctc                                                                         >  2:180259/1‑147 (MQ=1)
    gCTTGACGCGTCTTATCAGGACTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTT‑‑ACGC‑C‑GTATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                       <  1:758579/149‑1 (MQ=11)
    gCTTAACGCGTCTTATTATGCTTACGACAGTCCGGCGAACAGCCTGATCTAGCCCGATCTTGTCGGTCTGACACGGCGCTTCGC‑‑C‑C‑GTATCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                       <  2:290954/149‑1 (MQ=2)
      ttAACGCGTCTTATCAGGCCTACGCCAGACAGTGCAATAGCTTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                     <  1:219212/149‑1 (MQ=11)
       tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                    <  1:846817/149‑1 (MQ=255)
         aCGCGTTTTATCGGGCCTCCGCCAGACACCGTAATAGTCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTT‑‑CCGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                  <  2:261516/149‑1 (MQ=11)
            cgttctttCAGGCCTACGGCAGACAGCGCAATAGACTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                               <  2:460739/142‑1 (MQ=14)
            cgTCTTATCTGTCTTACGCCAGACAGCGCCATCGCGTTATTTAGCCTGATTTTGTAGGTCGGATAAGGCGTTTCCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                               <  2:460780/149‑1 (MQ=14)
                     aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAActg                                                      >  2:133326/1‑148 (MQ=255)
                      ggCCTACGCCAGACAGCGCAATACCCTGATTTAGTGTGATTTTTTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                     <  1:284554/149‑1 (MQ=17)
                       gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    <  1:265867/149‑1 (MQ=255)
                       gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    >  1:476119/1‑149 (MQ=255)
                         cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  <  1:841919/149‑1 (MQ=255)
                         cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  >  1:598050/1‑149 (MQ=255)
                          tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCACTGCCTGATGCGTCGATTGCCGCGTCTTATCAGGCTTATCTTAACCGTTGGTTAATTTTTCAAACTGTTcac                                                 >  1:113077/1‑149 (MQ=17)
                            cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                               <  1:136746/149‑1 (MQ=255)
                            cGCCAGACAGCCCAATAGCCTGATTTAGCGGCATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAAGGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                               <  1:297533/149‑1 (MQ=17)
                               cataccgcgcAATAGGCTGTTTGAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑C‑GCACCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                            <  2:22258/143‑1 (MQ=17)
                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  >  1:376805/1‑142 (MQ=255)
                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                           <  1:628668/149‑1 (MQ=255)
                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCGTca                                                  <  2:376805/142‑1 (MQ=14)
                                  gaaGCGCAATATCTTGATTTTGCGTGATTTATTAGGTCGGAGAAGGCTTTT‑‑ACGC‑C‑GCATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATTAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                         <  2:395713/147‑1 (MQ=18)
                                  caaGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                         <  1:464663/147‑1 (MQ=18)
                                        cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                   <  1:30765/149‑1 (MQ=255)
                                         aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                  <  1:131858/149‑1 (MQ=255)
                                          attgCCTGTTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                 <  1:91036/146‑1 (MQ=18)
                                           gtgCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCt                                <  2:442225/147‑1 (MQ=255)
                                             gCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCGCACCCGGTATCAACCCCttg                              >  2:180435/1‑148 (MQ=18)
                                             gCCTGATTGAGCGTGATTTTGTAGGTCGGATAAGTCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCGGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              <  1:798753/149‑1 (MQ=21)
                                              ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                             >  2:582648/1‑149 (MQ=255)
                                               cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACTCTTGCTGCGTCTTATCAGGCCTTTCTTATCCGTTGATTAATTTTTCAAACTCTTCATACCCGGTCTAAAAACCTTCCa                            >  1:628029/1‑149 (MQ=17)
                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                          <  1:630505/149‑1 (MQ=255)
                                                          tGATTTTGTAGGTCGGAT‑AGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                >  2:641263/1‑149 (MQ=25)
                                                             ttttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTATCTTAACCGTTGGTTAATTTATCAAACACTTCACACCCGGTATCAAACCCTTACATACAGATCAGGTTc              >  1:285737/1‑149 (MQ=21)
                                                              tttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTATCTTAACCGTTGGTTAATTTTTAAAACTCTTCACACCCGATATCAAAACCTTACATACCTCTCAGGTTCa             >  1:124725/1‑149 (MQ=21)
                                                               ttGTAGGTCGGATAAGGCATTG‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTAATTTGTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa            <  2:324431/149‑1 (MQ=25)
                                                                tctAGGTCGGATAAGTCGTTT‑‑ACGC‑C‑GTATCCGACATCAATGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc           <  2:634003/147‑1 (MQ=25)
                                                                  tAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATTCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa         <  1:370016/149‑1 (MQ=255)
                                                                  gaagTCGGATAAGGCGTTT‑‑ACGC‑C‑GTATCCGCAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa         <  2:221408/146‑1 (MQ=25)
                                                                     gTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCGCAGGTTTCAAACCAATCCATACAGCACAAGTTCAAcaactg      >  1:387835/1‑144 (MQ=18)
                                                                     gTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg      <  1:286604/149‑1 (MQ=255)
                                                                        gaaaaaGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACACTTGCCCCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTATTCAAACGCCGCACACCCAGTACCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg   <  2:410820/145‑1 (MQ=21)
                                                                         gATAAGGCGTTT‑‑ACGC‑C‑GCATCCGAGATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGGCTATCTTAACCGTTGGCTAATTTTTCAAACTCTTCACGCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc  <  1:812657/149‑1 (MQ=25)
                                                                                        ||                                                                                                                                        
CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_bme_pgi/4296296‑4296515

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