Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,207,942 (G)7→9 coding (442/630 nt) ycfK → e14 prophage; protein StfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,9351.G97.5% 138.2 / ‑3.2 40coding (435/630 nt)ycfKe14 prophage; protein StfP
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (22/17);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
*NC_0009131,207,9352.G97.5% 139.3 / ‑3.2 40coding (435/630 nt)ycfKe14 prophage; protein StfP
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (22/17);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGT  >  NC_000913/1207803‑1208065
                                                                                                                                       ||                                                                                                                                  
AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggggggg                                                                                                                             <  2:282608‑M1/140‑1 (MQ=255)
 GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCt‑‑gggggggg                                                                                                                            >  2:311739‑M1/1‑139 (MQ=255)
    gCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggggggat                                                                                                                         <  1:282608/140‑5 (MQ=255)
          cTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAAt                                                                                                                     >  1:331040/1‑140 (MQ=255)
            cAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATAc                                                                                                                   <  2:331040/140‑1 (MQ=255)
             aaTGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACa                                                                                                                  >  1:16534/1‑140 (MQ=255)
             aaTGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACa                                                                                                                  >  2:842550/1‑140 (MQ=255)
             aaTGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACa                                                                                                                  >  2:734674/1‑140 (MQ=255)
                         gCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCg                                                                                                      >  1:591208/1‑140 (MQ=255)
                                tAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                                                                                               <  2:62286/140‑1 (MQ=255)
                                                  tAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTg                                                                             >  2:443027/1‑140 (MQ=255)
                                                  tAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGGGGGGGCGGACTGCCGCAGAAGCATCTGCCTCT‑‑TGGGGGGATGGTAAGACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGcc                                                                           >  2:253055/1‑140 (MQ=255)
                                                       ataCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCgtt                                                                        <  2:827607/140‑1 (MQ=255)
                                                           cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    <  1:217899/140‑1 (MQ=255)
                                                             gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    >  1:575556/1‑138 (MQ=255)
                                                             gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    <  2:575556/138‑1 (MQ=255)
                                                                aTGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAAt                                                               <  1:364937/140‑1 (MQ=255)
                                                                aTGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGAGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAAt                                                               >  2:364937/1‑140 (MQ=255)
                                                                        tggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                                                       >  1:854855/1‑140 (MQ=255)
                                                                           tATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGt                                                    >  2:26588/1‑140 (MQ=255)
                                                                             ttGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa                                                  <  2:854855/140‑1 (MQ=255)
                                                                             ttGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa                                                  <  2:669717/140‑1 (MQ=255)
                                                                               gAGCAATCAGGAAATATTCAGGGGGGGGTGAATGCCCCCAAAGGTTCTTGCCCTGG‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                                <  2:663962/140‑1 (MQ=255)
                                                                                aGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                                               >  2:797514/1‑140 (MQ=255)
                                                                                  cAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                             <  1:734674/140‑1 (MQ=255)
                                                                                    aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                           <  1:797514/140‑1 (MQ=255)
                                                                                                 cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTg                              >  2:765952/1‑140 (MQ=255)
                                                                                                  aGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTgg                             <  1:443027/140‑1 (MQ=255)
                                                                                                   gTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                             <  1:159618/123‑1 (MQ=255)
                                                                                                   gTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                             >  2:159618/1‑123 (MQ=255)
                                                                                                        ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgg                       >  2:543020/1‑140 (MQ=255)
                                                                                                          tGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggg                     <  1:842550/140‑1 (MQ=255)
                                                                                                          gggAATGCCGCAAAAGTATCTGCCTTTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggg                     <  2:645682/139‑1 (MQ=255)
                                                                                                            gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCGCCTATCAATGCAGTGGCACCgggggg                   >  1:872086/1‑140 (MQ=255)
                                                                                                              aTGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGat                 >  1:174935/1‑140 (MQ=255)
                                                                                                                 ccGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGatatg              >  1:88181/1‑139 (MQ=255)
                                                                                                                 ccGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGatatg              >  2:828600/1‑139 (MQ=255)
                                                                                                                  cGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                    <  1:321086/133‑1 (MQ=255)
                                                                                                                  cGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                    >  2:321086/1‑133 (MQ=255)
                                                                                                                    cAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATGCaa           >  2:332192/1‑140 (MQ=255)
                                                                                                                    cAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATGCaa           >  1:173520/1‑140 (MQ=255)
                                                                                                                            tCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                    <  2:434190/123‑1 (MQ=255)
                                                                                                                            tCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                    >  1:434190/1‑123 (MQ=255)
                                                                                                                             cTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATGCAACCACTTCGt  <  1:765952/140‑1 (MQ=255)
                                                                                                                             cTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATGCAACCACTTCGt  <  1:332192/140‑1 (MQ=255)
                                                                                                                                       ||                                                                                                                                  
AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGT  >  NC_000913/1207803‑1208065

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