Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 7.2% K641Q (AAG→CAG)  abpB ← CP4‑57 prophage; putative helicase YfjK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6200TG 7.2% 138.2 / 3.1 56K641Q (AAG→CAG) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base T (34/18);  new base G (2/2);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

TTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTC  >  NC_000913/2761480‑2761760
                                                                                                                                            |                                                                                                                                            
ttCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGACAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACCTTCCTGCGCATACCACGCTTTa                                                                                                                                            >  2:753702/1‑143 (MQ=255)
  cggtcgtgTCGACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAGGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGTCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACa                                                                                                                                          <  2:255123/136‑1 (MQ=255)
      cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGATTGCCACTCTCTTCTAAGGCATTGCCTACAAAGACTATTTGTAAACGCTCGTACTTTAGCTGAACTTTTGTTTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                      <  1:941977/143‑1 (MQ=255)
       aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                     <  2:259283/143‑1 (MQ=255)
          tcgaCCAGTGGGTGTAACGTCTGGATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                  <  2:362073/140‑1 (MQ=255)
           tcaCCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtat                                                                                                                                 <  1:679008/141‑3 (MQ=255)
            tACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGtt                                                                                                                                <  2:363579/143‑1 (MQ=255)
             aCCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTCCAGCGTGGTGttt                                                                                                                               >  2:641050/1‑143 (MQ=255)
                aGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtatacgca                                                                                                                            <  1:509032/143‑8 (MQ=255)
                 gTCGGTGTAACGTCTGAATGGGAATGCCGATCTCTTCTAAGGCATTGACACCTGGCGGTAGATGGAAACTTTCGAAAGCCAATTTGACATGTGTGTAATCGACTTTGCGGCGTCTAGCAGGATTTACAGAGTGGTGTGTGACt                                                                                                                           >  2:809455/1‑143 (MQ=255)
                  tCGGTGTAACGGCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                          <  2:206232/143‑1 (MQ=255)
                           cGTCTGAATTGGAATTCCAATCTTTTCTAATGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTGCACATGTGTGTATTCGACTTTGCTGCGTATACCACGTTTTACAGCGTGGTGTTTAACTACATCCTCCa                                                                                                                 <  2:224465/143‑1 (MQ=255)
                            gTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACGTCATCCAt                                                                                                                >  2:906356/1‑143 (MQ=255)
                            gTCTGAATGGAAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACCATAGAACGCCAATGTGACATGTGTGTAATTGACTTGGCTGCGTATACCACGCTTTACAGCGTGGTGTGTAGCTAGAtcatccct                                                                                                                >  1:592225/1‑141 (MQ=255)
                               tGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCaa                                                                                                             >  1:790613/1‑143 (MQ=255)
                                gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCTTCGAg                                                                                                            >  1:204387/1‑143 (MQ=255)
                                gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAAATTTCGAACGCCAATTTGACCTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAg                                                                                                            >  1:416728/1‑143 (MQ=255)
                                  aTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGCGGTGTTTAACTACATCCTCCATCAAGCt                                                                                                          >  2:625074/1‑143 (MQ=255)
                                   tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTg                                                                                                         <  2:281662/143‑1 (MQ=255)
                                            cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAGTCGACTTTGCTGCGGATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCtga                                                                                                >  2:561874/1‑141 (MQ=255)
                                            cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                >  1:453273/1‑143 (MQ=255)
                                            cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATACTCCATCACGCTGCGCACTTta                                                                                                >  2:617622/1‑143 (MQ=255)
                                             aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTGCAGCGTGGTGTTTAACTACATCCTCCATCAACCTGAGCAGCTtag                                                                                               >  2:307837/1‑143 (MQ=255)
                                                ctctTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtgggta                                                                                            >  1:243959/1‑143 (MQ=255)
                                                 tcttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGCCTTTGCTGCGTATACCACGCGTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCgtaggtag                                                                                           >  2:896440/1‑143 (MQ=255)
                                                   ttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCGTAGGTAggg                                                                                         >  1:119164/1‑143 (MQ=255)
                                                       aaGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACCGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtaca                                                                                     >  2:71374/1‑141 (MQ=255)
                                                          gCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGcc                                                                                  >  1:358441/1‑143 (MQ=255)
                                                               gACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAgg                                                                             >  1:192409/1‑143 (MQ=255)
                                                                    cgtgtggtAGATGGAAACCTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGGGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTGAGGTAGGGTATAGCCGAAGGTATTg                                                                        <  2:21412/140‑1 (MQ=255)
                                                                    cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCACCCAGCTGCGCACCGTGGGTAGGGTATAGCCGCAGGTATTg                                                                        >  2:911892/1‑143 (MQ=255)
                                                                     tggtggGAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTgg                                                                       <  2:856847/143‑1 (MQ=255)
                                                                       gtggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGtt                                                                     >  1:662754/1‑143 (MQ=255)
                                                                         ggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTCTACTACGCTTTACAGCGTGGTGTTTAACTACATACTCCATCCACCTGAGCACCTTAGGTACGGTATAGCCGCAGGTATTgctgat                                                                   >  2:708140/1‑138 (MQ=255)
                                                                           tAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTAGAGCCGAAGGTATTGGTTATGa                                                                 >  1:783067/1‑143 (MQ=255)
                                                                            agagGGAAACTTTCGAACTCCAATTTGACATGCGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAg                                                                <  1:88839/143‑1 (MQ=255)
                                                                               tgggtaCTGTCGAACGCCAATTTGCCATGTGGGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATGCGCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCtt                                                             <  2:349335/138‑1 (MQ=255)
                                                                                ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCCAGCTGAGCCACTTAGGTACGGTCTAGCCGCCGGCACTGGCTCTGAGCCTc                                                            >  2:739135/1‑143 (MQ=255)
                                                                                   aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGGTATGAgactcacc                                                         >  2:19304/1‑136 (MQ=255)
                                                                                    aCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAAt                                                        <  1:854334/143‑1 (MQ=255)
                                                                                      tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGACGGTATTGGTTATGAGCTTCAAATCa                                                      >  2:513183/1‑143 (MQ=255)
                                                                                              gCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAg                                              <  1:61232/143‑1 (MQ=255)
                                                                                              gCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGCTACAGGTTATGAGCTCCAAATCATTGCTGAg                                              >  2:240910/1‑143 (MQ=255)
                                                                                                cAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                            >  1:315676/1‑143 (MQ=255)
                                                                                                  agttGTCATGTGTGGAATCGGCTGTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCa                                          <  2:263829/141‑1 (MQ=255)
                                                                                                    ttgttaTGTGTGTAATCGACTTTGCTGCGTATTCCACGCTTTTCAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAAt                                        <  2:92783/138‑1 (MQ=255)
                                                                                                           gtgtgtAGTGGACTTGGCGGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGGTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAaga                                 <  2:768946/143‑1 (MQ=255)
                                                                                                           gtgtgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGGAATGGAAaga                                 >  2:484372/1‑143 (MQ=255)
                                                                                                                tAATCGACCTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACg                            >  1:613212/1‑143 (MQ=255)
                                                                                                                   tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCGTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAACTCATTGTTGAGGGCAATGGAAAGAGTACGttt                         >  2:704876/1‑143 (MQ=255)
                                                                                                                      aCTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtctc                      >  1:665359/1‑143 (MQ=255)
                                                                                                                            ctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAACTCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTc                >  1:630381/1‑143 (MQ=255)
                                                                                                                              gcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTCTTGGTTCTGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTGCAg              >  2:432526/1‑143 (MQ=255)
                                                                                                                                  acaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCATCAGTAGGGAAAAGCCGAAGGTACTGGTCACCAGCTTCAAATAATTGTTGAGCGCAATGGAAAGAGTACGTTTATCCCCAACAGAGAt          >  2:156083/3‑143 (MQ=255)
                                                                                                                                     ccACGCTTTTCAGCGTGGTGTTTAACTTCACCCGCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGGGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa       <  2:28569/143‑1 (MQ=255)
                                                                                                                                     ccACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCATAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCCTGAAATCAGTGTTGAGGGCAATCGAAAGAGTACGGTTCTCCCCTTCCGacacccc       >  2:44135/1‑137 (MQ=255)
                                                                                                                                          gtttACAGCTTGGTGTTTAACTACATCCTCCATCAACCTGAGCATCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTc  <  2:279197/142‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTC  >  NC_000913/2761480‑2761760

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