Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,296,382 +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_bme_pgi4,296,3811.C87.5% 96.4 / 11.3 40intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (17/18);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_bme_pgi4,296,3812.G87.5% 98.2 / 11.4 40intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (17/18);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_bme_pgi/4296297‑4296519
                                                                                          ||                                                                                                                                            
gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                              >  2:852863/1‑149 (MQ=12)
 aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTATAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                             >  1:885910/1‑149 (MQ=255)
  cGCTTAACGCGTCTTATCAGGCCCACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                            <  1:90534/149‑1 (MQ=255)
   gCTTAACCCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTCGCGCGATTTTGTAGGTCGGATAAGGCGTT‑‑TCCG‑‑CCGCATCCGACATCAATGCCTCATGCGCAGCTTCACGCGTTTCATCATTCCTATCTTAACCTTTg                                                                           >  2:711869/1‑149 (MQ=11)
   ccTTAACGTCTCTTAACAGGCCTAAGCCAAACACCGCAATAGCCTGATTTAGCGTGATTTTATAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                           <  1:294246/148‑1 (MQ=11)
    cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTgg                                                                          >  2:108936/1‑149 (MQ=255)
     ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                         >  1:814467/1‑149 (MQ=255)
              cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                <  1:361625/149‑1 (MQ=255)
               ttATCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAGCGTCATTTTGTAGGTCGGATAAGGCGTT‑‑TGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGAGTCTTATAAGGCCTATCTTATCCGTTGGTTAATTTTTCa                                                               >  1:766466/1‑149 (MQ=14)
                    aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                          >  2:977778/1‑149 (MQ=255)
                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  1:892917/149‑1 (MQ=255)
                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  2:773398/149‑1 (MQ=14)
                        cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTGTGTTAATTTTTCAAtctcgcca                                                      >  1:228980/1‑145 (MQ=11)
                        cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTATTca                                                      >  2:355560/1‑149 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                  >  2:872690/1‑149 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                  >  2:144207/1‑149 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAATCTTCACAcc                                                  >  2:805798/1‑149 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                <  1:193574/149‑1 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTCGGTTACTTTTTCACACTCTTCACCCCCg                                                >  2:555401/1‑149 (MQ=14)
                               aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                               <  1:83110/149‑1 (MQ=255)
                                gACAGCGTAAAAGCCTAATATAGCGTGATTTTGTAGGCCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                              <  2:432476/149‑1 (MQ=17)
                                 aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                             <  1:946845/149‑1 (MQ=255)
                                       taaTAGCTTGATTTAGCGTGATTTTGTAGGTTGGATTAGGCGTT‑‑TTCG‑‑CCGCATACGACATCAATGACTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                       <  2:344450/148‑1 (MQ=18)
                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                      <  1:977778/149‑1 (MQ=255)
                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑CGCAGCATCCGACATCAATGCCTGATGCAACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGCTACTTTTTCAAACTCTTCACACCCGGTACCAAAcc                                      >  2:38371/1‑149 (MQ=17)
                                            gCCTGATTTAGCCGTGTTTTTTAGGTCATATAAGGCGTT‑‑T‑‑ACGCCGCTTCCGACATCAATGCCGTATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  2:650311/149‑1 (MQ=18)
                                            gCCTGATTTAGCGTTATTTTGTGGGTCCGACCCGGGGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  2:810355/149‑1 (MQ=21)
                                            gCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  1:241010/149‑1 (MQ=255)
                                               tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAAATCTTCACACCCGGTATCAAACCCTTCCAt                               >  2:237866/1‑149 (MQ=18)
                                                   ttAGCGGGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAGCGCTTCACACCCGGTATCAAACCCTTCCATACAg                           <  2:263750/149‑1 (MQ=21)
                                                         gTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCAACGATTCCGGTATCAAACCCTTCCATACAGCTCAg                      <  1:498252/149‑1 (MQ=18)
                                                             ttttGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAACGCATGATGCGACGCTTCCCGCGTCTAATCAGGCCTATCTTAACCGTTGGTGAATTTTTCAAACTCTTCTCACCCGGTATCAAACCCTGCCATACACCTCAgctca                  >  2:373151/1‑145 (MQ=18)
                                                               ttGAAGGTCCGATAAGGCGTG‑‑T‑‑ACGCCGCCTCCGACATCAATGCCCGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa                <  2:257142/149‑1 (MQ=25)
                                                                tGTAGGTCGGATAAGGCGATACC‑‑A‑‑CCGCATCCGACATCAACGCCTCATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc               <  1:182781/149‑1 (MQ=25)
                                                                     gTCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg          >  2:387442/1‑149 (MQ=255)
                                                                      tCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACACCCGTTATCAAACCCTTCCATACAACTCAGGTTCAACCatttc         >  2:91380/1‑145 (MQ=25)
                                                                      tCGGATAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGCTGGTTAATTTTTCAAACTCTTCACACCCCGTATCAAACCCTTACATAAAGCTCCGGTTCAACAAGTGc         >  1:464053/1‑149 (MQ=21)
                                                                         gATAAGGCGTT‑‑T‑‑ACGCCGTATCCGACATCAATGCCTGATGCGAAGCTGGCCGTCACTTATCAGGCCTATCGTTACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc      <  1:445765/149‑1 (MQ=255)
                                                                           tAAGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc    <  2:550751/149‑1 (MQ=255)
                                                                             aGGCGTT‑‑T‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt  >  1:802869/1‑149 (MQ=255)
                                                                                          ||                                                                                                                                            
GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_bme_pgi/4296297‑4296519

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: