Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 12.6% 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,6200TG12.6% 113.2 / 4.9 48K641Q (AAG→CAG) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base G (3/4);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTC  >  NC_000913/2761480‑2761760
                                                                                                                                            |                                                                                                                                            
ttCCAGTCAGTCTACCATTCGTTGTAACGTCTGAATGGGAATGCCAATCTCTGCTAAGGCATGGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTa                                                                                                                                            <  1:395118/143‑1 (MQ=255)
    aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGc                                                                                                                                        >  1:51360/1‑143 (MQ=255)
      cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                      <  1:608636/143‑1 (MQ=255)
                aGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAAc                                                                                                                            <  1:214699/143‑1 (MQ=255)
                  tCGGTGGAACGGCTGAATGGTCATGCCAATCTCGTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACTTGTGTGTAATCTACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                          <  2:339516/143‑1 (MQ=255)
                    ggggtAACGTCTGAATGGGAAACTCAAGCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACa                                                                                                                        <  2:375181/140‑1 (MQ=255)
                    ggTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACa                                                                                                                        <  1:26757/143‑1 (MQ=255)
                      tgtAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGTCATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                      <  1:54531/143‑1 (MQ=255)
                           cGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCa                                                                                                                 <  1:734534/143‑1 (MQ=255)
                             tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtatac                                                                                                                               >  2:617926/1‑123 (MQ=255)
                             tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtatac                                                                                                                               <  1:617926/127‑5 (MQ=255)
                                gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGGGTTTAACTACATCCTCCATCAAg                                                                                                            >  2:425774/1‑143 (MQ=255)
                                   tGGGAATGCCAATCTCTTCTAAGGCATTGACTCCTGGTGGTAGATGGATACTTGCGAACGCCAATGTGACACGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTg                                                                                                         <  2:573926/143‑1 (MQ=255)
                                    gggAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGATGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                        >  1:446487/1‑143 (MQ=255)
                                      gAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGaac                                                                                                      >  1:522389/1‑141 (MQ=255)
                                               tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGGGTGTAATCGACTTTGCTGCGTATACCACGCTCTACAGCGTGGTGTATAACTCAATCGTTGCTCAAGCTGCGCACCTtaggt                                                                                             >  2:309095/1‑143 (MQ=255)
                                                 tcttcgGAGGCGTTGACAGCTGGTGGTAGCTGGAAACGTTCGGACGCCAATTTGACATGTGTGTAATCGACTTTGCTCCCTATACCACGCTTTCCAGCGTGGTGTTTAACTCCATCCTCCATCAAGCTGAGCACCTtaggtag                                                                                           <  2:340491/143‑1 (MQ=255)
                                                       aaGGCATTGACACCTGGTGGTCGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACCTCCTCCATCAAGCTGAGCACCTTAGGTAGGGtata                                                                                     >  2:768532/1‑143 (MQ=255)
                                                        aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTGGCTAGGGTATAg                                                                                    >  1:243242/1‑143 (MQ=255)
                                                            aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGACAACTTCACATGTGTGGAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                >  1:540730/1‑143 (MQ=255)
                                                                           tAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCCTCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                 <  2:277571/143‑1 (MQ=255)
                                                                                  aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCATGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaa                                                          >  1:409605/1‑143 (MQ=255)
                                                                                    aCTTGCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAAt                                                        <  1:357681/143‑1 (MQ=255)
                                                                                      tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                      >  2:746739/1‑143 (MQ=255)
                                                                                        tctaaCCCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTAGTTATGAGCTTCAAATCAtt                                                    <  2:522311/140‑1 (MQ=255)
                                                                                        tCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                    <  1:445980/143‑1 (MQ=255)
                                                                                              gCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAg                                              >  2:415912/1‑143 (MQ=255)
                                                                                               tcgagttGAGATGTGTGTAATCGACTTTGCTGTGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCCCCATCAAGCTGAGCTCCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                             <  2:249098/138‑1 (MQ=255)
                                                                                               tccaTTTGACATGTGTGTAATCGACTTTGCTGCGTATATCTCGCTTTATAGCGTGGTGTTTAACTATATCCTACATCTATCTGTGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAATTAATTGTTGAgg                                             <  2:100011/140‑1 (MQ=255)
                                                                                                cAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                            >  1:525553/1‑143 (MQ=255)
                                                                                                cAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                            >  2:191472/1‑143 (MQ=255)
                                                                                                  aTTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGGTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGGAGGTAGTGGTTGTGAGCTTCAAATCCTTGTTGAGGGCa                                          >  1:101466/1‑143 (MQ=255)
                                                                                                         aTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaaa                                   >  2:320308/1‑143 (MQ=255)
                                                                                                           gtgtgtAATCGACTTTGCTGCGTATACCACGCTGTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGAAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAaga                                 <  2:336878/143‑1 (MQ=255)
                                                                                                              tgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCACGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                              >  1:767118/1‑143 (MQ=255)
                                                                                                                   tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGttt                         >  1:232223/1‑143 (MQ=255)
                                                                                                                    cGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAAACATTGTTGAGGGCAATGGAAAGAGAACGtctc                        >  2:674910/1‑143 (MQ=255)
                                                                                                                    cGACTGAGCTGCGTATACCACGCTGTACAGCGTGGTGTTTAACTACATCCGCCATCAAGCTGGGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtc                        <  1:347625/143‑1 (MQ=255)
                                                                                                                        tttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCGTAGGTAGGGGATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCGATGGAAAGAGTACGTTTCTccc                    >  1:583773/1‑143 (MQ=255)
                                                                                                                             tgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTGCAAATCATTGTTGAGGGCAATGGAATGAGTACGTTGCTCCCCTTCa               >  2:546405/1‑143 (MQ=255)
                                                                                                                              ttaaacaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAg              >  2:149330/7‑143 (MQ=255)
                                                                                                                              gcGTATTCCTCTCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAg              <  2:121188/143‑1 (MQ=255)
                                                                                                                              gcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCATTCAg              >  2:40800/1‑143 (MQ=255)
                                                                                                                                gTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGAAAAGAGTACGTTTCTCCCCTTCatac            >  2:1687/1‑140 (MQ=255)
                                                                                                                                ctagaCCACGCTGTACAGCGTGGTGTTGAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaa            <  2:496167/139‑1 (MQ=255)
                                                                                                                                  ataCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAAt          <  1:616338/143‑1 (MQ=255)
                                                                                                                                   caCCACGCTGTAAAGCGTGGTGTTTAACTACTTCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATc         >  2:484157/2‑143 (MQ=255)
                                                                                                                                      cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaaagacac      >  1:484948/1‑138 (MQ=255)
                                                                                                                                          ttttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTc  <  2:546813/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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑