Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,177,806 G→T V50F (GTT→TTT)  nusG → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,177,8060GT97.4% 118.3 / NA 39V50F (GTT→TTT) nusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (23/15);  total (24/15)
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

TGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGA  >  NC_000913_3_pae_tpiA/4177660‑4177951
                                                                                                                                                  |                                                                                                                                                 
tGTCTGCAGCTCCTAAACAGCGCTTGTACGGCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGttt                                                                                                                                                 <  2:233004/149‑1 (MQ=255)
tGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGttt                                                                                                                                                 <  2:202381/149‑1 (MQ=255)
 gTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTg                                                                                                                                                <  2:359030/149‑1 (MQ=255)
    tGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTTAATTCATGTTACCAACCGAAGAAGTGTTTGaaa                                                                                                                                             <  2:123959/149‑1 (MQ=255)
        gCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGGGTTTGAAATCCg                                                                                                                                         >  1:530791/1‑149 (MQ=255)
           ccTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGAATATCAAATTACACAACATGGAAGATTTGTTTGGTGAGGTCATGGTACCAACGTAGGAAGTGTTTGAAATccatag                                                                                                                                      >  1:270976/1‑145 (MQ=255)
              aaaaaGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCgg                                                                                                                                   >  2:421333/1‑149 (MQ=255)
                     gcTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGt                                                                                                                            >  1:108981/1‑149 (MQ=255)
                      cTGGTACGTCGGTCACGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGGGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTc                                                                                                                           <  1:8736/149‑1 (MQ=255)
                             gtcgtTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATAACAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAg                                                                                                                    <  1:421333/149‑1 (MQ=255)
                                    aGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAAt                                                                                                                                            >  1:475965/1‑118 (MQ=255)
                                    aGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAAt                                                                                                                                            <  2:475965/118‑1 (MQ=255)
                                         ttttACGATTTTGGAGGCCGGGAAGCGACGTCACTGCGGGAGCAGGGCAAGTTACACAACATGGAAGAATTGGTTGGGGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAAtt                                                                                                        >  2:340212/1‑149 (MQ=255)
                                            tCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAAttctt                                                                                                     >  2:292224/1‑149 (MQ=255)
                                               ggTTGTGAAGGCCGCGGAGGAACGTCGCGGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTccc                                                                                                  <  2:366139/149‑1 (MQ=255)
                                                gTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAAGATGGCAGATTTGTTTGGTGAAGTCGTGGTACAAAGCGAAGAAGCGTTTGAAAACAGTGGCGTTCAGCGTCGCAAGAACGAACGTACATTCATCCCt                                                                                                 >  2:186604/1‑149 (MQ=255)
                                                            gggTAGCAACGTCGCTGCGTGACCATCTTAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGGGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTc                                                                                     <  2:112322/147‑1 (MQ=255)
                                                            gcgTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTc                                                                                     >  1:316407/1‑149 (MQ=255)
                                                                   aaCGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTCGCTACGTCCTCGTTCCGa                                                                              >  2:415438/1‑149 (MQ=255)
                                                                      gTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGGTgg                                                                           <  2:108981/149‑1 (MQ=255)
                                                                       tCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATggg                                                                          >  2:142713/1‑148 (MQ=255)
                                                                           tGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGtgatg                                                                      >  1:384245/1‑149 (MQ=255)
                                                                                       tCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGc                                                          <  1:572481/149‑1 (MQ=255)
                                                                                        cAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCt                                                         <  1:331891/149‑1 (MQ=255)
                                                                                                  aacaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGt                                               <  2:335470/149‑1 (MQ=255)
                                                                                                  aacaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGt                                               >  1:95173/1‑149 (MQ=255)
                                                                                                       ggAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCa                                          >  1:300538/1‑149 (MQ=255)
                                                                                                       ggAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTCGTgcgcc                                          >  2:397418/1‑148 (MQ=255)
                                                                                                                 tttgGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACcgcg                                >  2:430235/1‑149 (MQ=255)
                                                                                                                 tttgGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACcgcg                                <  1:468234/149‑1 (MQ=255)
                                                                                                                   tgGTGAAGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCgtg                              <  1:427998/149‑1 (MQ=255)
                                                                                                                         aGTCATGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATgg                        >  1:55855/1‑149 (MQ=255)
                                                                                                                              tGGTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTc                   <  1:492871/149‑1 (MQ=255)
                                                                                                                                gTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                 >  1:39366/1‑149 (MQ=255)
                                                                                                                                gTACCAACCGAAGAAGTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                 >  1:327865/1‑149 (MQ=255)
                                                                                                                                     aaCCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGGGATGGGCTTCATcggcc            >  1:472251/1‑148 (MQ=255)
                                                                                                                                         gaagacGTGTTTAAAATCCGTGGCGGTCAGGATCGCAAAAGCGTACGTAAATTCTTTCCTGGCTCCGTCCTCGTTCCGACGGTGAAGAACGACGTGGGCTGGCCTCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTAc        <  2:417201/149‑1 (MQ=255)
                                                                                                                                          aagaagTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACt       >  1:450914/1‑149 (MQ=255)
                                                                                                                                             aagTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAATGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCAGTACTTcc    >  1:509207/1‑149 (MQ=255)
                                                                                                                                             aagTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTcc    >  2:76822/1‑149 (MQ=255)
                                                                                                                                             aagTGTTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGAGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTCCTTcc    >  2:362432/1‑149 (MQ=255)
                                                                                                                                               gtgtTTGAAATCCGCGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCGTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa  <  2:450198/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATCAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGA  >  NC_000913_3_pae_tpiA/4177660‑4177951

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: ‑