Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,623 C→A 6.9% V640L (GTA→TTA)  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,6230CA 6.9% 143.1 / 4.0 58V640L (GTA→TTA) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base C (33/21);  new base A (2/2);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-02

CCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGC  >  NC_000913/2761482‑2761764
                                                                                                                                             |                                                                                                                                             
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)
                                                                                                                                            aCAGCGTGGTGTTTAACGACTTCCTCCGTCAATCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAATAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGc  <  2:253948/143‑1 (MQ=255)
                                                                                                                                            aCAGCGTCGTGTTTAACTACATCCTCCAGCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGc  <  2:433579/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGC  >  NC_000913/2761482‑2761764

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: ‑