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

GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  NC_000913_3_bme_pgi/4296297‑4296519
                                                                                          ||                                                                                                                                            
gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                              >  2:852863/1‑149 (MQ=12)
 aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTATAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                             >  1:885910/1‑149 (MQ=255)
  cGCTTAACGCGTCTTATCAGGCCCACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                            <  1:90534/149‑1 (MQ=255)
   gCTTAACCCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTCG‑CGCGATTTTGTAGGTCGGATAAGGCGTT‑‑TCCG‑‑C‑C‑GCATCCGACATCAATGCCTCATGCGCAGCTTCACGCGTTTCATCATTCCTATCTTAACCTTTg                                                                           >  2:711869/1‑149 (MQ=11)
   ccTTAACGTCTCTTAACAGGCCTAAGCCAAACACCGCAATAGCCTGATTTAG‑CGTGATTTTATAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                           <  1:294246/148‑1 (MQ=11)
    cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTgg                                                                          >  2:108936/1‑149 (MQ=255)
     ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTCGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                         >  1:814467/1‑149 (MQ=255)
              cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                <  1:361625/149‑1 (MQ=255)
               ttATCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAG‑CGTCATTTTGTAGGTCGGATAAGGCGTT‑‑TGCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGAGTCTTATAAGGCCTATCTTATCCGTTGGTTAATTTTTCa                                                               >  1:766466/1‑149 (MQ=14)
                    aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                          >  2:977778/1‑149 (MQ=255)
                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  1:892917/149‑1 (MQ=255)
                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                        <  2:773398/149‑1 (MQ=14)
                        cTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTGTGTTAATTTTTCAAtctcgcca                                                      >  1:228980/1‑145 (MQ=11)
                        cTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTATTca                                                      >  2:355560/1‑149 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                  >  2:872690/1‑149 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                  >  2:144207/1‑149 (MQ=255)
                            gCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAATCTTCACAcc                                                  >  2:805798/1‑149 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                <  1:193574/149‑1 (MQ=255)
                              cAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTCGGTTACTTTTTCACACTCTTCACCCCCg                                                >  2:555401/1‑149 (MQ=14)
                               aGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                               <  1:83110/149‑1 (MQ=255)
                                gACAGCGTAAAAGCCTAATATAG‑CGTGATTTTGTAGGCCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                              <  2:432476/149‑1 (MQ=17)
                                 aCAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                             <  1:946845/149‑1 (MQ=255)
                                       taaTAGCTTGATTTAG‑CGTGATTTTGTAGGTTGGATTAGGCGTT‑‑TTCG‑‑C‑C‑GCATACGACATCAATGACTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                       <  2:344450/148‑1 (MQ=18)
                                        aaTAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                      <  1:977778/149‑1 (MQ=255)
                                        aaTAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑CGCAGCATCCGACATCAATGCCTGATGCAACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGCTACTTTTTCAAACTCTTCACACCCGGTACCAAAcc                                      >  2:38371/1‑149 (MQ=17)
                                            gCCTGATTTAGCCGTG‑TTTTTTAGGTCATATAAGGCGTT‑‑T‑‑ACGC‑C‑GCTTCCGACATCAATGCCGTATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  2:650311/149‑1 (MQ=18)
                                            gCCTGATTTAG‑CGTTATTTTGTGGGTCCGACCCGGGGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  2:810355/149‑1 (MQ=21)
                                            gCCTGATTCAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                  <  1:241010/149‑1 (MQ=255)
                                               tGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAAATCTTCACACCCGGTATCAAACCCTTCCAt                               >  2:237866/1‑149 (MQ=18)
                                                   ttAG‑CGGGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAGCGCTTCACACCCGGTATCAAACCCTTCCATACAg                           <  2:263750/149‑1 (MQ=21)
                                                         gTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCAACGATTCCGGTATCAAACCCTTCCATACAGCTCAg                      <  1:498252/149‑1 (MQ=18)
                                                             ttttGTAGGTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAACGCATGATGCGACGCTTCCCGCGTCTAATCAGGCCTATCTTAACCGTTGGTGAATTTTTCAAACTCTTCTCACCCGGTATCAAACCCTGCCATACACCTCAgctca                  >  2:373151/1‑145 (MQ=18)
                                                               ttGAAGGTCCGATAAGGCGTG‑‑T‑‑ACGC‑C‑GCCTCCGACATCAATGCCCGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa                <  2:257142/149‑1 (MQ=25)
                                                                tGTAGGTCGGATAAGGCGATACC‑‑A‑‑C‑C‑GCATCCGACATCAACGCCTCATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc               <  1:182781/149‑1 (MQ=25)
                                                                     gTCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg          >  2:387442/1‑149 (MQ=255)
                                                                      tCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACACCCGTTATCAAACCCTTCCATACAACTCAGGTTCAACCatttc         >  2:91380/1‑145 (MQ=25)
                                                                      tCGGATAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGCTGGTTAATTTTTCAAACTCTTCACACCCCGTATCAAACCCTTACATAAAGCTCCGGTTCAACAAGTGc         >  1:464053/1‑149 (MQ=21)
                                                                         gATAAGGCGTT‑‑T‑‑ACGC‑C‑GTATCCGACATCAATGCCTGATGCGAAGCTGGCCGTCACTTATCAGGCCTATCGTTACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc      <  1:445765/149‑1 (MQ=255)
                                                                           tAAGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc    <  2:550751/149‑1 (MQ=255)
                                                                             aGGCGTT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt  >  1:802869/1‑149 (MQ=255)
                                                                                          ||                                                                                                                                            
GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAG‑CGTGATTTTGTAGGTCGGATAAGGCGTT‑‑T‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCT  >  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: ‑