Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 581,185 C→T 10.9% L118L (CTG→TTG)  nohD → DLP12 prophage; DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913581,1850CT10.9% 154.6 / 10.2 55L118L (CTG→TTG) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base C (30/19);  new base T (3/3);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCC  >  NC_000913/581049‑581307
                                                                                                                                        |                                                                                                                          
cAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGt                                                                                                                         >  1:1191841/1‑140 (MQ=32)
         aGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCg                                                                                                                <  1:820257/140‑1 (MQ=32)
           aTCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTGCGTGCTGTCGCGGATCGCa                                                                                                              >  1:126757/1‑140 (MQ=17)
              tcCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                                                           >  1:1121517/1‑140 (MQ=32)
                   ccagGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                                                                      <  1:814122/140‑1 (MQ=32)
                   ccagGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                                                                      >  2:513951/1‑140 (MQ=32)
                   ccagGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                                                                      >  2:1616770/1‑140 (MQ=32)
                   ccagGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACGTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                                                                      >  2:513958/1‑140 (MQ=17)
                     agGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg                                                                                                    >  1:1661590/1‑140 (MQ=32)
                       ggACTATTGTGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGcc                                                                                                  >  1:880061/1‑140 (MQ=17)
                       ggACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGcc                                                                                                  >  1:1131548/1‑140 (MQ=32)
                       ggACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGcc                                                                                                  >  2:1240703/1‑140 (MQ=32)
                       ggACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGcc                                                                                                  >  2:1241155/1‑140 (MQ=32)
                        gACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCa                                                                                                 <  1:944101/140‑1 (MQ=32)
                           tATTGAGTACGAACGCCATCAACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTa                                                                                              >  1:943022/1‑140 (MQ=17)
                            aTTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAt                                                                                             >  1:845282/1‑140 (MQ=32)
                             ttGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAtt                                                                                            >  1:1001292/1‑140 (MQ=32)
                             ttGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAtt                                                                                            >  1:775954/1‑140 (MQ=32)
                               gAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATtct                                                                                          <  2:289459/140‑1 (MQ=32)
                                  tACGAACGCCATCTACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       >  1:628542/1‑140 (MQ=17)
                                  tACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       >  1:856141/1‑140 (MQ=32)
                                  tACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       >  1:539485/1‑140 (MQ=32)
                                  tACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       <  1:547037/140‑1 (MQ=32)
                                  tACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       >  2:201367/1‑140 (MQ=32)
                                  tACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       <  1:103830/140‑1 (MQ=32)
                                  tACGAACGCCATCGACTTACGCGAGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGa                                                                                       <  1:691143/140‑1 (MQ=17)
                                     gAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACgg                                                                                    >  2:103771/1‑140 (MQ=32)
                                     gAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACgg                                                                                    >  2:1086395/1‑140 (MQ=32)
                                       aCGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGa                                                                                  >  1:908797/1‑140 (MQ=32)
                                       aCGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGa                                                                                  >  2:705205/1‑140 (MQ=32)
                                                aCTTACGCGTGCGCAGGCCGATGCACCGGAGCTGAAAAATGCCAGAGCCTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGTATCCCCCTGt                                                                         <  1:1334896/140‑1 (MQ=12)
                                                  ttACGCGTGCGCAGGCCGTTGCACAGGAGCTGAAAATTGCCAGAGACTCCGCTGATGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:698179/140‑1 (MQ=12)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAATAATGCCAGAGACTCCGCTGTAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGTTGAAAGTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:109853/140‑1 (MQ=255)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:1014172/140‑1 (MQ=32)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:932569/140‑1 (MQ=32)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:698336/140‑1 (MQ=32)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:1362209/140‑1 (MQ=32)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:540493/140‑1 (MQ=32)
                                                  ttACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTCGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCg                                                                       <  2:21280/140‑1 (MQ=17)
                                                   tACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCgg                                                                      <  1:1540726/140‑1 (MQ=32)
                                                   tACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCgg                                                                      >  2:1642390/1‑140 (MQ=32)
                                                    aCGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGt                                                                     >  1:1462941/1‑140 (MQ=32)
                                                         tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAgc                                                                >  1:64758/1‑140 (MQ=32)
                                                         tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAgc                                                                >  1:1285861/1‑140 (MQ=32)
                                                          gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAgcg                                                               >  1:693748/1‑140 (MQ=32)
                                                           cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAgcgg                                                              >  1:1264540/1‑140 (MQ=32)
                                                             cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAgcggcg                                                            <  2:744179/140‑1 (MQ=32)
                                                                    aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCg                                                     >  2:98701/1‑140 (MQ=31)
                                                                            gaacTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGTTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCg                                                                       <  2:1676760/111‑1 (MQ=1)
                                                                            gaacTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGTTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCg                                                                       <  2:1676755/111‑1 (MQ=1)
                                                                            gaacTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGTTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCg                                                                       >  1:1676755/4‑114 (MQ=1)
                                                                            gaacTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGTTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCg                                                                       >  1:1676760/4‑114 (MQ=1)
                                                                            gaacTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGTTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCg                                                                       >  1:1075252/4‑114 (MQ=1)
                                                                            gaacTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGTTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCg                                                                       <  2:1075252/111‑1 (MQ=1)
                                                                                                                       cagtCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTACCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGcc  <  2:29891/137‑1 (MQ=1)
                                                                                                                                        |                                                                                                                          
CAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCC  >  NC_000913/581049‑581307

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