Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,951,412 A→G L6L (CTA→CTG 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,4120AG93.9% 94.9 / ‑3.1 35L6L (CTA→CTGyecEDUF72 domain‑containing protein YecE
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (11/20);  total (11/24)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACC  >  NC_000913/1951167‑1951525
                                                                                                                                                                                                                                                      |                                                                                                                  
cgatgtgcggttACGCCGCCGAGGTATCGATACAAGTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATTATGTACATCGGTCTGCCGc                                                                                                                <  1:133693/240‑1 (MQ=255)
       cggttACGCCGCCGAGTTATCGATACAATTGTGTTATGTGGGATCTCGACCAATAGCGGTGGTGAATCCTCCGCCCGCAATGCCTGGGAACTGGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAGAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCTTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCg                                                                                                        <  1:65340/248‑1 (MQ=255)
         gttACGCCGCCGAGGTATCGATACAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGc                                                                                                       >  1:82473/2‑250 (MQ=255)
                    gAGGTATCGATACAAGTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACAGCGTTCTGCCGCAATGGTCGCATCCTAAATggn                                                                                          <  1:18911/251‑2 (MQ=255)
                      ggTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTg                                                                                        <  1:138799/251‑1 (MQ=255)
                              cagactgtgctcttccgatctctcGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                                                <  1:79384/231‑1 (MQ=255)
                               tcAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCGACCCGCTCATCGCCCGTGTGCGTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                                                <  1:132636/249‑1 (MQ=255)
                               aCAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTTGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                                                <  1:73708/250‑1 (MQ=255)
                               aCAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                                                <  1:7045/250‑1 (MQ=255)
                                      tgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTTGAAGAGATCCTTCACGCGGTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAc                                                                          <  1:101796/249‑1 (MQ=255)
                                                     cGACCAATCTCGGTGGTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTTGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTGCACGAGTTATGATTTACATCGGTCTGCCTCAATGGTCGCATCCTAAATTGGTTCGGGTGGGTATCACCATCCTTGAAGAGTAt                                                          <  1:49294/250‑1 (MQ=255)
                                                     cGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg                                                         >  1:28264/1‑251 (MQ=255)
                                                      gACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc                                                        >  1:26449/1‑251 (MQ=255)
                                                       aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTGCACGCGTTATGATTTACATCTGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGcc                                                       <  1:69877/251‑1 (MQ=255)
                                                        ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCACCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                                                >  1:53542/1‑225 (MQ=255)
                                                        ccAAGATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAc                                                                          <  1:72487/231‑1 (MQ=255)
                                                        ccAAGATCGGTGTTGAAGCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAGACGACTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccct                                                     <  1:15784/251‑2 (MQ=255)
                                                                     tGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCTAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTTCCGCAATGGTCg                                                                                                        <  1:25610/189‑1 (MQ=255)
                                                                               gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAATAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                                                                         <  1:125776/209‑1 (MQ=255)
                                                                                   gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                                                                        <  1:55780/190‑1 (MQ=255)
                                                                                    cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGAGTgacgcgggcatt                          >  1:48018/1‑240 (MQ=255)
                                                                                         ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGCAGAGTATGCCCGCCACTTTAACTGCGTgacgcgggcattttaa                      >  1:151208/1‑235 (MQ=255)
                                                                                                 cTTGGATTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTgacgcgggcattttaaaaatcact              >  1:4047/1‑227 (MQ=255)
                                                                                                        ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAAGCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTg                                                                  <  1:142672/191‑1 (MQ=255)
                                                                                                             tggtggATTGCCAAAGACGCCGGTAGGCCCGCCAGCGCCGAGCAGCACAAAAACAGTAGCAACCATACCCCCCCGCGCATCGCCGGTGTGCGTAGCGGGGAAGATACCCTCACCGCGTTATGTTTTACACCGTTCTCCCGCAAGGGTCCATTCCTAAAGGGGTGGGTTTGGGGACCACCAGCCTTGAAGAGTATCCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAAcc  <  1:4390/245‑1 (MQ=255)
                                                                                                                                ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAag                                                               <  1:139612/170‑1 (MQ=255)
                                                                                                                                 cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGGTGTGTATCACCAGCCTTGAAGAGTATGct                                                       <  1:51332/177‑2 (MQ=255)
                                                                                                                                           gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTg                                                                                  <  1:9408/140‑1 (MQ=255)
                                                                                                                                             cAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                                                                          >  1:167136/1‑130 (MQ=255)
                                                                                                                                                       naacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagcg                                                             >  1:19354/4‑147 (MQ=255)
                                                                                                                                                                          atcATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAAGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                                                                        >  1:40797/3‑119 (MQ=255)
                                                                                                                                                                          atcATATCTACCCGCGCATCGCCCGTGTGCNTAGCGTGGANNAGNTCCTTNANNNGNNNTGNTTTACATCGGTCTGCCGCANTGNNCGCATCNTAAATGGGTGCGGTTGGGTATCACCa                                                                        <  2:40797/117‑1 (MQ=255)
                                                                                                                                                                                        cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGAATCATAAATGGGTGCGGTTGGGTATCACCAg                                                                       >  1:128772/1‑106 (MQ=255)
                                                                                                                                                                                        cgcATCGCCCGNGTGCGTAGCGNGGANNAGATNCNTNACGCGTTATGATTTACATCGGTCTGCCGCNATGGTCGCATCCTAAATGGGTGCNNTTNNNTANNACCAg                                                                       <  2:128772/106‑1 (MQ=255)
                                                                                                                                                                                                    gtgCGTAGCGGGGAAGAGGTCCTTCACGCGTTATGATTTGCATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                                                                            <  1:160080/89‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                  
GGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACC  >  NC_000913/1951167‑1951525

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: