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‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_bme_pgi/4296296‑4296515
                                                                                        ||                                                                                                                                        
cGACGCCTGCCGCGTCTTATGAGGCCGAAGCCAGGCAGCGCAATGGCCAGATTGCGGGTGATTTTTTAGGTCGGGCCAGGGGTTT‑‑ACGCCGCATACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                           <  2:369265/149‑1 (MQ=2)
  aCGCTTAACGCGTCTTATCAGGCCTACCCCAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTT‑‑A‑‑CGCAGAATCCGACATCACTCCCCGATGCGACGCTTCCCGCGTCCTCTCCGGCCTCTCTTAAcctc                                                                         >  2:180259/1‑147 (MQ=1)
    gCTTGACGCGTCTTATCAGGACTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTT‑‑ACGCCGTATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                       <  1:758579/149‑1 (MQ=11)
    gCTTAACGCGTCTTATTATGCTTACGACAGTCCGGCGAACAGCCTGATCTAGCCCGATCTTGTCGGTCTGACACGGCGCTTCGC‑‑CCGTATCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                       <  2:290954/149‑1 (MQ=2)
      ttAACGCGTCTTATCAGGCCTACGCCAGACAGTGCAATAGCTTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                     <  1:219212/149‑1 (MQ=11)
       tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                    <  1:846817/149‑1 (MQ=255)
         aCGCGTTTTATCGGGCCTCCGCCAGACACCGTAATAGTCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTT‑‑CCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                  <  2:261516/149‑1 (MQ=11)
            cgttctttCAGGCCTACGGCAGACAGCGCAATAGACTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                               <  2:460739/142‑1 (MQ=14)
            cgTCTTATCTGTCTTACGCCAGACAGCGCCATCGCGTTATTTAGCCTGATTTTGTAGGTCGGATAAGGCGTTTCCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                               <  2:460780/149‑1 (MQ=14)
                     aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAActg                                                      >  2:133326/1‑148 (MQ=255)
                      ggCCTACGCCAGACAGCGCAATACCCTGATTTAGTGTGATTTTTTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                     <  1:284554/149‑1 (MQ=17)
                       gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    <  1:265867/149‑1 (MQ=255)
                       gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    >  1:476119/1‑149 (MQ=255)
                         cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  <  1:841919/149‑1 (MQ=255)
                         cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  >  1:598050/1‑149 (MQ=255)
                          tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCACTGCCTGATGCGTCGATTGCCGCGTCTTATCAGGCTTATCTTAACCGTTGGTTAATTTTTCAAACTGTTcac                                                 >  1:113077/1‑149 (MQ=17)
                            cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                               <  1:136746/149‑1 (MQ=255)
                            cGCCAGACAGCCCAATAGCCTGATTTAGCGGCATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAAGGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                               <  1:297533/149‑1 (MQ=17)
                               cataccgcgcAATAGGCTGTTTGAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGCCGCACCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                            <  2:22258/143‑1 (MQ=17)
                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  >  1:376805/1‑142 (MQ=255)
                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                           <  1:628668/149‑1 (MQ=255)
                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCGTca                                                  <  2:376805/142‑1 (MQ=14)
                                  gaaGCGCAATATCTTGATTTTGCGTGATTTATTAGGTCGGAGAAGGCTTTT‑‑ACGCCGCATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATTAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                         <  2:395713/147‑1 (MQ=18)
                                  caaGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                         <  1:464663/147‑1 (MQ=18)
                                        cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                   <  1:30765/149‑1 (MQ=255)
                                         aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                  <  1:131858/149‑1 (MQ=255)
                                          attgCCTGTTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                 <  1:91036/146‑1 (MQ=18)
                                           gtgCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCt                                <  2:442225/147‑1 (MQ=255)
                                             gCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCGCACCCGGTATCAACCCCttg                              >  2:180435/1‑148 (MQ=18)
                                             gCCTGATTGAGCGTGATTTTGTAGGTCGGATAAGTCGTTT‑‑ACGCCGCATCCGACATCAATGCCGGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              <  1:798753/149‑1 (MQ=21)
                                              ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                             >  2:582648/1‑149 (MQ=255)
                                               cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACTCTTGCTGCGTCTTATCAGGCCTTTCTTATCCGTTGATTAATTTTTCAAACTCTTCATACCCGGTCTAAAAACCTTCCa                            >  1:628029/1‑149 (MQ=17)
                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                          <  1:630505/149‑1 (MQ=255)
                                                          tGATTTTGTAGGTCGGATAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                >  2:641263/1‑149 (MQ=25)
                                                             ttttGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTATCTTAACCGTTGGTTAATTTATCAAACACTTCACACCCGGTATCAAACCCTTACATACAGATCAGGTTc              >  1:285737/1‑149 (MQ=21)
                                                              tttGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTATCTTAACCGTTGGTTAATTTTTAAAACTCTTCACACCCGATATCAAAACCTTACATACCTCTCAGGTTCa             >  1:124725/1‑149 (MQ=21)
                                                               ttGTAGGTCGGATAAGGCATTG‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTAATTTGTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa            <  2:324431/149‑1 (MQ=25)
                                                                tctAGGTCGGATAAGTCGTTT‑‑ACGCCGTATCCGACATCAATGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc           <  2:634003/147‑1 (MQ=25)
                                                                  tAGGTCGGATAAGGCGTTT‑‑ACGCCGCATTCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa         <  1:370016/149‑1 (MQ=255)
                                                                  gaagTCGGATAAGGCGTTT‑‑ACGCCGTATCCGCAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa         <  2:221408/146‑1 (MQ=25)
                                                                     gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCGCAGGTTTCAAACCAATCCATACAGCACAAGTTCAAcaactg      >  1:387835/1‑144 (MQ=18)
                                                                     gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg      <  1:286604/149‑1 (MQ=255)
                                                                        gaaaaaGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACACTTGCCCCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTATTCAAACGCCGCACACCCAGTACCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg   <  2:410820/145‑1 (MQ=21)
                                                                         gATAAGGCGTTT‑‑ACGCCGCATCCGAGATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGGCTATCTTAACCGTTGGCTAATTTTTCAAACTCTTCACGCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc  <  1:812657/149‑1 (MQ=25)
                                                                                        ||                                                                                                                                        
CGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  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: