Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,388 G→T 8.5% M1M (ATG→ATT) † yfjY → CP4‑57 prophage; putative DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,3880GT 8.5% 147.8 / 3.4 47M1M (ATG→ATT) †yfjYCP4‑57 prophage; putative DNA repair protein
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base T (2/2);  total (21/26)
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.06e-01

TTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGAC  >  NC_000913/2776256‑2776526
                                                                                                                                    |                                                                                                                                          
ttGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAAc                                                                                                                                     <  2:138815/140‑1 (MQ=255)
ttGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAAc                                                                                                                                     <  2:135070/140‑1 (MQ=255)
         gCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCa                                                                                                                            >  2:396322/1‑140 (MQ=255)
             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATccc                                                                                                                        >  1:14301/1‑140 (MQ=255)
             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATccc                                                                                                                        >  2:414453/1‑140 (MQ=255)
             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATccc                                                                                                                        >  2:242676/1‑140 (MQ=255)
               aCGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCcac                                                                                                                      >  1:428999/1‑140 (MQ=255)
                  cctgcctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                   <  2:723170/140‑1 (MQ=255)
                  cctgcctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                   <  2:260453/140‑1 (MQ=255)
                  cctgcctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                   <  2:175977/140‑1 (MQ=255)
                   ctgcctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGac                                                                                                                  <  1:750317/140‑1 (MQ=255)
                    tgcctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGaca                                                                                                                 >  2:56473/1‑140 (MQ=255)
                    tgcctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGaca                                                                                                                 >  2:56455/1‑140 (MQ=255)
                        tCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGTGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCgg                                                                                                             <  2:855271/140‑1 (MQ=255)
                        tCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCgg                                                                                                             <  2:854599/140‑1 (MQ=255)
                        tCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCgg                                                                                                             <  2:827032/140‑1 (MQ=255)
                          aCCGAATGCGATGCCATGTCGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTa                                                                                                           <  1:386912/140‑1 (MQ=255)
                           ccGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATTATGGAACAGTCACTCATCCCACAGACACCGGTAc                                                                                                          >  1:1172073/1‑140 (MQ=255)
                           ccGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATTATGGAACAGTCACTCATCCCACAGACACCGGTAc                                                                                                          >  1:1171626/1‑140 (MQ=255)
                                     tGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGAc                                                                                                >  1:1063517/1‑140 (MQ=255)
                                                 gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCAc                                                                                    >  1:813409/1‑140 (MQ=255)
                                                     tATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTa                                                                                >  1:514252/1‑140 (MQ=255)
                                                        tACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTaaaa                                                                             >  2:555471/1‑140 (MQ=255)
                                                          ccGCCTGACGGATTACGCTCTGCAGCACCCTGAATCCAGCCCCATAATGCGCATCATTGATTAAGGAGACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                           <  2:723300/140‑1 (MQ=255)
                                                              cTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCcaca                                                                                                                     <  1:1025119/94‑1 (MQ=255)
                                                              cTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCcaca                                                                                                                     >  2:1025119/1‑94 (MQ=255)
                                                               tGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCt                                                                      >  1:650849/1‑140 (MQ=255)
                                                                  gggATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTaa                                                                   <  2:411708/140‑1 (MQ=255)
                                                                         cgcccgccagccCCCCGCAGCCCGCCCCATGAGGCGCATTAGTCATAAAGGATCCCCATGAGGGACCAGTCCCTCAGCCCTCCGACCCCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCATCACgctgct                                                            <  2:800052/133‑1 (MQ=255)
                                                                         cGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                            <  2:14301/140‑1 (MQ=255)
                                                                                          aaTCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGTCACCGGTACTTCCACTGAACGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGcgcg                                           >  2:738934/1‑140 (MQ=255)
                                                                                                gcgcCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAg                                     >  1:739694/1‑140 (MQ=255)
                                                                                                  gcCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGc                                   <  1:400419/140‑1 (MQ=255)
                                                                                                   ccATTATGCGCATTATTGATTAAGGATACTCATTATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCg                                  <  2:1172073/140‑1 (MQ=255)
                                                                                                   ccATTATGCGCATTATTGATTAAGGATACTCATTATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCg                                  <  2:1171626/140‑1 (MQ=255)
                                                                                                   ccATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCg                                  <  2:507143/140‑1 (MQ=255)
                                                                                                   ccATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCg                                  <  2:543882/140‑1 (MQ=255)
                                                                                                   ccATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCg                                  <  1:69620/140‑1 (MQ=255)
                                                                                                          gcgcATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAtt                           <  2:784721/140‑1 (MQ=255)
                                                                                                           cgcATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTc                          >  1:255590/1‑140 (MQ=255)
                                                                                                           cgcATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTc                          >  1:665020/1‑140 (MQ=255)
                                                                                                              attattGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAcc                       >  2:375850/1‑140 (MQ=255)
                                                                                                              attattGATTAAGGATACTCATGATGGAACAGTCAATCATCACACAGAGACCGGTACGTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAcc                       <  2:456680/140‑1 (MQ=255)
                                                                                                                tattGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCTGTACTTCCACTTACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTc                     <  2:355947/140‑1 (MQ=255)
                                                                                                                    gATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACt                 <  1:21600/140‑1 (MQ=255)
                                                                                                                       tAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAg              >  2:1040338/1‑140 (MQ=255)
                                                                                                                        aaGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAgg             <  1:957702/140‑1 (MQ=255)
                                                                                                                                   tgatgGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGAc  <  1:603930/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                          
TTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGAC  >  NC_000913/2776256‑2776526

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: