Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,750,073 G→T 10.2% intergenic (+13/+42) yfjD → / ← grpE putative inner membrane protein/nucleotide exchange factor GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,750,0730GT10.2% 71.7 / 2.3 38intergenic (+13/+42)yfjD/grpEputative inner membrane protein/nucleotide exchange factor GrpE
Reads supporting (aligned to +/- strand):  ref base G (13/19);  new base T (2/2);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

CCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCA  >  NC_000913/2749935‑2750214
                                                                                                                                          |                                                                                                                                             
ccTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGCATCCGCTTGCCGAGAGTGTGGCGGAGTAACGAAAACGGCCACGCAt                                                                                                                                           >  1:21039/1‑143 (MQ=255)
ccTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCAt                                                                                                                                           >  1:232586/1‑143 (MQ=255)
         gggaCCGGCGGGCGTGTGGGAGAGTACGACATCGATACTCTCGACTGACAGGACAATATGATTAAGCAGGTAGACGTTTTTCCTGTGCAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCa                                                                                                                                  <  2:360074/140‑1 (MQ=255)
                  gtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTGTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTACCGAAAACGGACCGGCAATCGCATGCAGGGCCCTGa                                                                                                                         >  2:401340/1‑143 (MQ=255)
                        attctcGAGTAAGATAAGGATATCCGGGACGTACAGGACAATAGGAATAAGCAGGAAAAAGCGTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAattatt                                                                                                                   <  1:186552/138‑1 (MQ=255)
                         gtGGCGAGTACGAGATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGATGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                  <  1:53796/142‑1 (MQ=255)
                             cGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaa                                                                                                              >  1:498083/1‑143 (MQ=255)
                                 tACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCGGGGCCGGGAATTATTACGAAACCa                                                                                                          >  1:700894/1‑143 (MQ=255)
                                     atagcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGaaa                                                                                                      <  1:242443/143‑1 (MQ=255)
                                      tatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGGATGATTACGAAAGCAGAAAt                                                                                                     >  2:734181/1‑143 (MQ=255)
                                         cgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                  <  2:670398/143‑1 (MQ=255)
                                            tatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATTCGCATGCAGGGCCGTGAATTATTCCGAAAGCAGAAATTAAGCt                                                                                               >  1:700892/1‑143 (MQ=255)
                                               tctcGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTCCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCtttt                                                                                            <  1:707601/143‑1 (MQ=255)
                                                ctcGACGTCCAGCACACTATGATTAAGCAGGTAACAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTACCGAAAACGGCCCGACATTCGCATGCAGGGCCGTGAGCTACTACGAAAGCAGAAACTAAGCTGTTg                                                                                           >  2:833774/1‑143 (MQ=255)
                                                 tgtaggTACAGGACAATATGATTAGGCAGGTAAAAGTTGTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGc                                                                                          <  2:587973/138‑1 (MQ=255)
                                                         cagggaaaTATGCTTAACAAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTa                                                                                  <  2:329066/137‑1 (MQ=255)
                                                           tggcAATATGAGTAAGCAGGTACGAGTTTTTCTTGTGAGACCGCTGCGCGAGAGTGTGGCGGAGTAACTAAAACGGCCCGGTATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACa                                                                                <  2:474895/140‑1 (MQ=255)
                                                            gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGCGTAACGAGCAAGGCCCAGCATTCGCACGCAGGGCCGTGAATTATTACGAAAGCAGCAATTAAGCTTTTGCTTTCGCTACAg                                                                               >  1:390210/1‑143 (MQ=255)
                                                               aaTATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTaa                                                                            <  2:521909/143‑1 (MQ=255)
                                                                atagGATTAAGCAGGGAAAAGTTTTTCCGGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAAc                                                                           <  2:155934/143‑1 (MQ=255)
                                                                     atctaGCAGGTAAATGTTGTTCCTGGTAAACCGCTGCGCGAGAGTGTGGGTGAGTCACGAAAACGGCCCTGCAGTCGCATGCAGGGCCGTGAATTATTACTAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATcg                                                                      <  2:231473/139‑1 (MQ=255)
                                                                     aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACCGTAACCATcg                                                                      >  1:293068/1‑143 (MQ=255)
                                                                      ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATcgc                                                                     >  2:299037/1‑143 (MQ=255)
                                                                                   aGTTTTTCCTGTGAAGCCGCTGCGCGAGAGTGGGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTa                                                        <  2:50856/143‑1 (MQ=255)
                                                                                      ttttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGCATCGTACGa                                                     >  1:183566/1‑143 (MQ=255)
                                                                                                 aaCCGCTGCGACATGATGTTGCGGTGTAACTAACATGGCCCGGCATTCGCATGCAGGGCCTTGAATTATTACGAAAGCATAAATTTAGCTTTTGCTTTCGTTACAGTAACCATCGCCGCACGAATCTTACGACCATTCAGCGt                                          <  2:837779/143‑1 (MQ=255)
                                                                                                  aCCGCTGCGCGAGCGGGTGGCGGAGTAACGAAAACGGCCCGCCATTCGCATGCAGCGCCGCGACTTATTACGAAAGCAGAAATTAAGCTCTTGCTTTCGCTACAGGAACCATCCCCACACGAATCGTACGACCAGTCAGcaga                                         >  2:424851/1‑140 (MQ=255)
                                                                                                  aCCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACCGTAACCATCGCCGCACGACTCGTACGACCATTCAGCGtc                                         >  2:317423/1‑142 (MQ=255)
                                                                                                  aCCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGta                                         >  2:572527/1‑143 (MQ=255)
                                                                                                      cTGCGCGCGAGTCTGGCGTAGTAACGAAAACGGCCCGGCAGTCGCATGCAGGGCCGTGAATTATTACGAACGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAc                                     <  1:8000/143‑1 (MQ=255)
                                                                                                            gagaGTGTGCCGGAGTAACGAACACGGCCCGGCATTCGCATGCAGGGCCGTGAATTCTTTCGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCt                               <  2:786235/143‑1 (MQ=255)
                                                                                                                 tgtgTCTGAGTAACGAAAACGGCCCGGCATTCACATGCAGGGCCGTGAATTAGTACGAAAGCAGAAATTCCGCCCTCGCTTTCGCTACAGTCACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATa                          <  1:641853/143‑1 (MQ=255)
                                                                                                                      cGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGcc                     <  1:822625/143‑1 (MQ=255)
                                                                                                                        gAGTAACGAAAACGGCCCTGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCa                   <  2:700892/143‑1 (MQ=255)
                                                                                                                           tAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTa                <  2:232586/143‑1 (MQ=255)
                                                                                                                              aattcaCGGCCCTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGt             >  1:418000/6‑143 (MQ=255)
                                                                                                                                   ccGGCCGGGCATTCCCATGCAGGGCCGTGAATTATCACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCt        <  1:110484/142‑1 (MQ=255)
                                                                                                                                         cGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCa  >  1:494999/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                             
CCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCA  >  NC_000913/2749935‑2750214

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: