Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 400,157 T→C Y256C (TAC→TGC)  ddlA ← D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913400,1570TC97.5% 150.5 / ‑6.2 40Y256C (TAC→TGC) ddlAD‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/16);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGG  >  NC_000913/400010‑400304
                                                                                                                                                   |                                                                                                                                                   
aaaaaCACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAgg                                                                                                                                                  >  1:168011/1‑151 (MQ=255)
  aaaCACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCa                                                                                                                                                <  1:642053/151‑1 (MQ=255)
   aaCACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCAt                                                                                                                                               >  1:31158/1‑151 (MQ=255)
        gTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGaaa                                                                                                                                          <  1:165331/151‑1 (MQ=255)
         tCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAAt                                                                                                                                         >  1:452938/1‑151 (MQ=255)
         tCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGGGTCGCAGGCATAGAAAt                                                                                                                                         >  1:131116/1‑151 (MQ=255)
            acacGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGc                                                                                                                                      >  2:120122/1‑151 (MQ=255)
                  gcCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTg                                                                                                                                 >  2:300593/1‑150 (MQ=255)
                        ccTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGa                                                                                                                          <  1:396482/151‑1 (MQ=255)
                             gcATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGAtctct                                                                                                                     <  2:592445/151‑1 (MQ=255)
                                        gttGTTAAGCCTGACCGGCAATGGACCGGATCTTATCGGTGATTTCTGGCGCAATGGCTGCTGGAACTACCACTTTCGCGCAATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCTACAGGTGc                                                                                                           >  1:376357/2‑150 (MQ=255)
                                         ttGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTg                                                                                                         >  2:600320/1‑151 (MQ=255)
                                             tAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCtt                                                                                                     >  2:266439/1‑151 (MQ=255)
                                              aaGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCtt                                                                                                     <  1:258362/150‑1 (MQ=255)
                                                gCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCg                                                                                                  <  2:60758/151‑1 (MQ=255)
                                                 ccTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCgg                                                                                                 >  1:548561/1‑151 (MQ=255)
                                                  cTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCgg                                                                                                 >  1:177463/1‑150 (MQ=255)
                                                   tGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGAt                                                                                               >  1:412175/1‑151 (MQ=255)
                                                           aaTCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTg                                                                                       <  1:624892/151‑1 (MQ=255)
                                                               gCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGccc                                                                                    >  2:36922/1‑150 (MQ=255)
                                                               gCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCa                                                                                   <  2:165146/151‑1 (MQ=255)
                                                                    gATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGCAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAAc                                                                               >  1:66102/1‑150 (MQ=255)
                                                                        ttATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTgcgc                                                                          >  2:217926/1‑151 (MQ=255)
                                                                              ttGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCg                                                                     >  1:186421/1‑150 (MQ=255)
                                                                                     cTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACg                                                             >  2:320714/1‑151 (MQ=255)
                                                                                                   gCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATccc                                                 <  1:148599/149‑1 (MQ=255)
                                                                                                      ggAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGc                                             >  1:62995/1‑150 (MQ=255)
                                                                                                       gAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTc                                           >  1:628810/1‑151 (MQ=255)
                                                                                                       gAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTc                                           >  2:321131/1‑151 (MQ=255)
                                                                                                                cTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCAc                                  >  1:99493/1‑151 (MQ=255)
                                                                                                                    cgcgCCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTa                              <  2:325892/151‑1 (MQ=255)
                                                                                                                       gcCATCTTCGTCAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGa                           <  2:49870/151‑1 (MQ=255)
                                                                                                                                  cAATGTACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAtcgaacgcga                 <  1:78842/150‑1 (MQ=255)
                                                                                                                                       tACTTGGTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGcc            >  2:442397/1‑150 (MQ=255)
                                                                                                                                            ggTGTCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATc       <  2:526559/150‑1 (MQ=255)
                                                                                                                                               gtCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGAc    <  1:229528/150‑1 (MQ=255)
                                                                                                                                               gtCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGAc    <  1:502626/150‑1 (MQ=255)
                                                                                                                                               gtCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACg   <  2:520074/151‑1 (MQ=255)
                                                                                                                                                tCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACgg  >  1:395642/1‑151 (MQ=255)
                                                                                                                                                tCGCAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACg   <  1:617614/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGG  >  NC_000913/400010‑400304

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