Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,951,448 G→T G18G (GGG→GGT)  yecE → DUF72 domain‑containing protein YecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,951,4480GT87.0% 46.8 / ‑1.9 23G18G (GGG→GGT) yecEDUF72 domain‑containing protein YecE
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (9/11);  total (10/13)
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.11e-01

TGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCC‑GCATGTCACCGTTGGCTCCGCGCATTG  >  NC_000913/1951205‑1951680
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
tgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTTGAAGAGATCCT‑TCACGCGGTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCAc                                                                                                                                                                                                                                      <  1:101796/249‑1 (MQ=255)
               cGACCAATCTCGGTGGTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTTGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑GCACGAGTTATGATTTACATCGGTCTGCCTCAATGGTCGCAT‑CCTAAATTGGTTCGGGTGGGTATCACCATCCTTGAAGAGTAt                                                                                                                                                                                                                      <  1:49294/250‑1 (MQ=255)
               cGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg                                                                                                                                                                                                                     >  1:28264/1‑251 (MQ=255)
                gACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc                                                                                                                                                                                                                    >  1:26449/1‑251 (MQ=255)
                 aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑GCACGCGTTATGATTTACATCTGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGcc                                                                                                                                                                                                                   <  1:69877/251‑1 (MQ=255)
                  ccAAGATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCAc                                                                                                                                                                                                                                      <  1:72487/231‑1 (MQ=255)
                  ccAAGATCGGTGTTGAAGCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGA‑GACGACTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccct                                                                                                                                                                                                                 <  1:15784/251‑2 (MQ=255)
                                         gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAATAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCAcc                                                                                                                                                                                                                                     <  1:125776/209‑1 (MQ=255)
                                              cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGAGTgacgcgggcatt                                                                                                                                                                                      >  1:48018/1‑240 (MQ=255)
                                                   ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGCAGAGTATGCCCGCCACTTTAACTGCGTgacgcgggcattttaa                                                                                                                                                                                  >  1:151208/1‑235 (MQ=255)
                                                           cTTGGATTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTgacgcgggcattttaaaaatcact                                                                                                                                                                          >  1:4047/1‑227 (MQ=255)
                                                                  ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAAGCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTg                                                                                                                                                                                                                              <  1:142672/191‑1 (MQ=255)
                                                                       tggtggATTGCCAAAGACGCCGGTAGGCCCGCCAGCGCCGAGCAGCACAAAAACAGTAGCAACCATACCCCCCCGCGCATCGCCGGTGTGCGTAGCGGGGAAGATACCCT‑CACCGCGTTATGTTTTACACCGTTCTCCCGCAAGGGTC‑CATTCCTAAAGGGGTGGGTTTGGGGACCACCAGCCTTGAAGAGTATCCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAAcc                                                                                                                                                              <  1:4390/245‑1 (MQ=255)
                                                                                          ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAag                                                                                                                                                                                                                           <  1:139612/170‑1 (MQ=255)
                                                                                           cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGGTGTGTATCACCAGCCTTGAAGAGTATGct                                                                                                                                                                                                                   <  1:51332/177‑2 (MQ=255)
                                                                                                                 naacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagcg                                                                                                                                                                                                                         >  1:19354/4‑147 (MQ=255)
                                                                                                                                    atcATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAA‑GCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCa                                                                                                                                                                                                                                    >  1:40797/3‑119 (MQ=255)
                                                                                                                                    atcATATCTACCCGCGCATCGCCCGTGTGCNTAGCGTGGANNAGNTCCT‑TNANNNGNNNTGNTTTACATCGGTCTGCCGCANTGNNCGCAT‑CNTAAATGGGTGCGGTTGGGTATCACCa                                                                                                                                                                                                                                    <  2:40797/117‑1 (MQ=255)
                                                                                                                                                  cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGAAT‑CATAAATGGGTGCGGTTGGGTATCACCAg                                                                                                                                                                                                                                   >  1:128772/1‑106 (MQ=255)
                                                                                                                                                  cgcATCGCCCGNGTGCGTAGCGNGGANNAGATNCN‑TNACGCGTTATGATTTACATCGGTCTGCCGCNATGGTCGCAT‑CCTAAATGGGTGCNNTTNNNTANNACCAg                                                                                                                                                                                                                                   <  2:128772/106‑1 (MQ=255)
                                                                                                                                                              gtgCGTAGCGGGGAAGAGGTCCT‑TCACGCGTTATGATTTGCATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATc                                                                                                                                                                                                                                        <  1:160080/89‑1 (MQ=255)
                                                                                                                                                                                                                              aT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTgacgcggg                                                                                                                                                                                          >  2:9605/1‑65 (MQ=255)
                                                                                                                                                                                                                                    aaaTGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTCACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCCGCATGTCACCGTTGGCTCCGCGCATTg  >  1:74773/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
TGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCC‑GCATGTCACCGTTGGCTCCGCGCATTG  >  NC_000913/1951205‑1951680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑