Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,996 G→T 16.2% intergenic (+98/‑372) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT16.2% 102.1 / 5.2 46intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (5/3);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGC  >  NC_000913/2560848‑2561143
                                                                                                                                                    |                                                                                                                                                   
aGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                                     <  8:414114/149‑1 (MQ=255)
aGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                                     <  8:339645/149‑1 (MQ=255)
  ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGtt                                                                                                                                                   <  1:332039/149‑1 (MQ=255)
  ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTGGCCGtt                                                                                                                                                   <  1:203319/149‑1 (MQ=255)
        cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                             <  8:188867/149‑1 (MQ=255)
           cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTAAAAGGGTGCCGTTAAGGTGGGt                                                                                                                                          >  2:55543/1‑149 (MQ=255)
                  cAGTTATTATTAAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTTTGGTTTGTAAAAAGTTACAGTTAACGTGTGAAAAAAGt                                                                                                                                   >  3:305009/1‑149 (MQ=255)
                    gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTg                                                                                                                                 >  3:30180/1‑149 (MQ=255)
                      tattatTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTGGGGTGGGTCAAATGTTGCCGTTAAGGTGGGTGAAAAGTTg                                                                                                                                 >  1:31230/1‑147 (MQ=255)
                          attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGACGtg                                                                                                                            >  3:322186/1‑147 (MQ=255)
                          attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                                                >  8:82080/1‑128 (MQ=255)
                          attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                                                <  7:82080/128‑1 (MQ=255)
                                gCCGTACAACATGAAATAACGGAATCAGGCCGGCAATGTTCCCGCCTTTTCTTCGCCTTAAATCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGCGGAGGCTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                     <  2:75277/149‑1 (MQ=255)
                                 ccGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                              >  6:13565/1‑123 (MQ=255)
                                 ccGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                              <  5:13565/123‑1 (MQ=255)
                                    tACAGCATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTAAAATGTAGCCGTAAAGggggaac                                                                                                                  >  3:133527/1‑145 (MQ=255)
                                       aacaTGAAATAAGGGACTCAGGCGGGCATTGTTACCGCATTTACTTTGCCTTAATTCAATGCATTCAGGTGGGTAAATATTTGCCCTTAGGGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                              <  1:55543/149‑1 (MQ=25)
                                        acaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCAATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGAAACAGTTGCCGTTAACGTGGCTCCAAAGTAGCCCGCTAGGTGGGTCAAAAg                                                                                                             >  5:15440/1‑149 (MQ=255)
                                          aTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                           >  7:67734/1‑149 (MQ=255)
                                                  aCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAATGTTGGGTCAAAAGTTGCCGTTAAGGTGGCTCAAATTTTGCAGTATAGGTGTCGCAAAAGTTgactcga                                                                                                    >  8:294862/1‑142 (MQ=25)
                                                           ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTAAAAAGTTGCCTTTTAGGTGGGTCAAAATTTGAAGTTAAGGTGTGTGAAATATTTAAGAGTAAGTAttat                                                                                           >  5:47065/1‑146 (MQ=25)
                                                                         ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTAAAATGTTGCCGATAAGGTGGGGGAAAATTTGCGGAGCATGAAGTCGACAATTATAAATTc                                                                             >  5:148559/1‑148 (MQ=21)
                                                                          cGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                              <  6:285934/82‑1 (MQ=255)
                                                                          cGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                              >  5:285934/1‑82 (MQ=255)
                                                                          cGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAATGTTGCCTTTAAGGAGGGTCAAAATTTTAAGTATAGGGGTCGCTAAATTTGCAGAGAAAGAATAGTTCTATTAAGAAtgc                                                                             >  4:7603/1‑145 (MQ=16)
                                                                                     tGCCTTAAATCCCGGCATTCCGTTGGGTAAATAGCAGCCGTTCCGGTTGGGTAAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACggg                                                                  <  5:325388/147‑1 (MQ=21)
                                                                                         ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                                         <  2:43661/72‑1 (MQ=255)
                                                                                             attaaTGTATTCAGGTTGGTTTAAATTTGCCTTTACGGGTGGGTAAAAACTTGCCTTTAAGGTGGGTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                         <  8:410707/144‑1 (MQ=25)
                                                                                                      ttCAGGTGGGTAAATATCTGGCTTTAGTGGGGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggt                                               <  3:222726/149‑1 (MQ=37)
                                                                                                       tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCTAATTTTGCCTAGTAAGAATTCTACTATTAGGAATTCTAGGCGTCGGGATTTTTTTAGCCTCGcggtc                                              >  7:232958/1‑149 (MQ=25)
                                                                                                         aGGTGGGTCAAAAGTCGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGGGTAAAAATTTGCAGTTAAGGTGGGGCAAAAGTTGACGAGAAAGAATTCTCAAATTAATAATTAAAGGCGTCGGGAGTTTTTTATCATAAcggtcga                                            >  6:237663/1‑148 (MQ=255)
                                                                                                           gTGGGGCAAAAGCCGCCGTTAGTGGTGCGTCAAAACCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                          <  1:412698/149‑1 (MQ=38)
                                                                                                                              ttgtggtgGGTCAAAAGTTGCCTTTCAGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga                         <  3:7603/145‑1 (MQ=255)
                                                                                                                                gtggtgGCTAAAATCTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt                      <  2:339298/148‑1 (MQ=39)
                                                                                                                                    tggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTACCGAGTAAGAATTCTAAAATTAAGATTACTAGGCTACCGGATATAATTAGCCTCCCGGGCGGACATAGAGAGAAGCCGAGAGTTAAACa                 >  2:355863/1‑149 (MQ=21)
                                                                                                                                     gggTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa                   >  7:208646/1‑146 (MQ=255)
                                                                                                                                       gTCAAAAGTTGCCGTTAAGGTTGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGGTGCCGAGTAAGAAATCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa              <  1:355863/149‑1 (MQ=255)
                                                                                                                                       gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa              <  8:208646/149‑1 (MQ=255)
                                                                                                                                       gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                >  8:395256/1‑147 (MQ=255)
                                                                                                                                       gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                >  5:303758/1‑147 (MQ=255)
                                                                                                                                          aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa           >  2:357518/1‑149 (MQ=255)
                                                                                                                                            aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCTGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAat          >  1:407630/1‑148 (MQ=255)
                                                                                                                                            aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAat          >  6:241229/1‑148 (MQ=255)
                                                                                                                                              gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat       >  7:57757/1‑149 (MQ=255)
                                                                                                                                                  ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                                          >  7:65447/1‑62 (MQ=255)
                                                                                                                                                  ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                                          <  8:65447/62‑1 (MQ=255)
                                                                                                                                                  ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc    <  7:395256/148‑1 (MQ=255)
                                                                                                                                                   cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc  <  1:357518/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGC  >  NC_000913/2560848‑2561143

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