Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,925 T→G 5.2% intergenic (+80/‑137) yfjQ → / → yfjR CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,9250TG 5.2% 1486.3 / 8.5 445intergenic (+80/‑137)yfjQ/yfjRCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (51/365);  new base G (26/3);  total (77/368)
Fisher's exact test for biased strand distribution p-value = 6.88e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 452 total aligned reads displayed.

GCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCAT  >  CP009273/2762858‑2762991
                                                                   |                                                                  
gCCACGTTGTTAATATCTGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGt                                                                    <  1:26351811/68‑1 (MQ=255)
    cGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTtata                                                                <  1:30737990/68‑1 (MQ=255)
    cGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAAGGCCCCGGGtagg                                                                >  1:8988180/1‑66 (MQ=21)
     gttgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATGCGCCGGGTGTAg                                                               >  1:30381477/1‑68 (MQ=21)
     gttgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTtata                                                                >  1:3751548/1‑67 (MQ=255)
     gttgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGGTagag                                                               >  1:6231026/1‑68 (MQ=255)
       tgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGc                                                              <  1:25558780/67‑1 (MQ=255)
       tgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGc                                                              <  1:11524434/67‑1 (MQ=255)
       tgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCt                                                             <  1:2576860/68‑1 (MQ=255)
       tgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCt                                                             <  1:3165333/68‑1 (MQ=255)
       tgttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCt                                                             <  1:8214702/68‑1 (MQ=255)
              atCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                                                       >  1:11351114/1‑67 (MQ=255)
               tCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGGTATAGCTTTTAAcc                                                      >  1:2658429/1‑67 (MQ=255)
                cGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAcccc                                                    <  1:7584565/68‑1 (MQ=255)
                 ggACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAccccc                                                   <  1:31181722/68‑1 (MQ=255)
                 ggACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAccccc                                                   <  1:17198991/68‑1 (MQ=255)
                 ggACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAccccc                                                   <  1:25313468/68‑1 (MQ=255)
                   acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCg                                                  <  1:28955208/67‑1 (MQ=255)
                   acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCg                                                  <  1:6488971/67‑1 (MQ=255)
                       caccTGTCCGCATCGCTATGTGCTCGTGTGCCTCTATCCCCCGGTTATAGCTTTTAACCCCCGTTACa                                             <  1:28726318/68‑1 (MQ=255)
                       caccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACa                                             <  1:10466611/68‑1 (MQ=255)
                       caccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACa                                             <  1:22555766/68‑1 (MQ=255)
                         ccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACa                                             >  1:9811552/1‑66 (MQ=255)
                                  aTCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACa                                             >  1:13767164/1‑57 (MQ=255)
                                         tgtgCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAAt                           <  1:21195154/68‑1 (MQ=255)
                                            gCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATaa                         <  1:8078405/67‑1 (MQ=255)
                                            gCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATaa                         <  1:36992906/67‑1 (MQ=255)
                                            gCTCGTGTGCCTCAATCCCCCGGGTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATaa                         <  1:28565269/67‑1 (MQ=255)
                                             cTCGTGTGCCTCAATCCCCCGGTTAGAGCTTTTAGCCGCCGGTACTGCTGGGTTTTGGAGAAATaaca                       >  1:34433013/1‑66 (MQ=255)
                                             cTCGTGTGCCTCAATCCCCCGGGTAGCGCGTTTGACGGCCGTTGCCGCGGGCTGTTGCagaaagaaga                       >  1:13262401/1‑66 (MQ=255)
                                              tCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATaaaa                       <  1:19676177/67‑1 (MQ=255)
                                               cGTGTGCCTCAATCCCCCGGTTATAGCTTTTTGCGGCGGTTACATGTGGGGTTTGCAGATATaaaaaa                     >  1:18845293/1‑68 (MQ=255)
                                               cGTGTGCCTCAATCCCCCGGGTATAGCTTTTAACCCGCGTTGCATGTGGGGTTTGCTGATATaaaaaa                     >  1:32976709/1‑68 (MQ=255)
                                                 tgtgCCTCAATCCCCCGGTTATGGCGTTTGTGCCCCGTTCGGTCTGGCGGTTGGAGAAGTaataaagt                   >  1:3502652/1‑66 (MQ=255)
                                                 tgtgCCTCAATCCCCCGGTTATAGCTTTTACCCGCCGGTACATCTGGCGGTTGCTGAATTTAataata                   >  1:19198898/1‑68 (MQ=255)
                                                 tgtgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGGTAGATCTGGGTTTTGCAGAAATAAAAAATa                   >  1:30186702/1‑68 (MQ=38)
                                                 tgtgCCTCAATCCCCCGGGTATAGCTTTTAAGCCCCGTTACGTCTGGGGTTTGCTGAATTAAAAAATa                   >  1:27720158/1‑68 (MQ=255)
                                                         aaTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGc           <  1:24045629/68‑1 (MQ=255)
                                                          aTCCCCCGGGTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCg          >  1:24588755/1‑68 (MQ=255)
                                                             ccccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTg       <  1:32627086/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:7364127/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:29680407/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:14294785/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:1418013/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:5908276/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:16410024/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:12478378/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:17333726/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:17793045/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:24914768/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:108436/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:7883487/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:32178086/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:22242927/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:21926236/68‑1 (MQ=255)
                                                              cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCGTTTGCAGAAATAAAAAATAGTTTCTGCGTTGt      <  1:28156794/68‑1 (MQ=255)
                                                               cctgTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:21742497/65‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTTGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:4807062/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGGTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:10082519/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:19519711/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:20107633/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:36607587/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:4091984/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:439392/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:3516559/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:4608941/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:18516395/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:1314774/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:34594514/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:20333464/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:8378349/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:8800930/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:34007169/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:27397157/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:22006596/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:17971594/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:23051355/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:23288228/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:23439587/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:24270002/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:24801769/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:25156905/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:16744744/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:16000705/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:27097888/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:33546851/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:28331512/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:21547911/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:15724811/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:15160204/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:29113268/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:29401149/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:14853311/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:31388775/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:20831055/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:31969389/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:33413820/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTCGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:36030296/68‑1 (MQ=255)
                                                               ccGGTTATAGCTTTTAACCCCCGTTACAGCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTc     <  1:5430207/68‑1 (MQ=255)
                                                                 ggTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCa   >  1:10598989/1‑68 (MQ=255)
                                                                   |                                                                  
GCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCAT  >  CP009273/2762858‑2762991

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