Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,320 G→T 6.0% P99P (CCG→CCT)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3200GT 6.0% 2240.3 / 89.4 643P99P (CCG→CCT) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (309/294);  new base T (26/13);  total (336/307)
Fisher's exact test for biased strand distribution p-value = 6.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 645 total aligned reads displayed.

GTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGG  >  CP009273/2762253‑2762387
                                                                   |                                                                   
gtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCg                                                                     >  1:3592073/1‑68 (MQ=255)
  cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                                   >  1:26104866/1‑68 (MQ=255)
   gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                                  <  1:12010672/68‑1 (MQ=255)
       gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                              >  1:29228381/1‑68 (MQ=255)
          tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                             >  1:31835362/1‑66 (MQ=255)
           aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                          <  1:33399659/68‑1 (MQ=255)
            cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                          >  1:10054902/1‑67 (MQ=255)
            cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                          >  1:17935102/1‑67 (MQ=255)
             tAAGCTTATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                        <  1:22183774/68‑1 (MQ=255)
             tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                        <  1:5937424/68‑1 (MQ=255)
             tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                        <  1:18388551/68‑1 (MQ=255)
             tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                        <  1:21641196/68‑1 (MQ=255)
              aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                        <  1:12825980/64‑1 (MQ=255)
              aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                        <  1:33750752/64‑1 (MQ=255)
              aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                        <  1:40736/64‑1 (MQ=255)
              aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                                       >  1:34094856/1‑68 (MQ=255)
               aGCATATGCTGCGCCTACGTCGGGAAGTGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                                      >  1:14572712/1‑68 (MQ=255)
                gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:16375650/1‑68 (MQ=255)
                gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:15246836/1‑68 (MQ=255)
                gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:12625600/1‑68 (MQ=255)
                gCATATGCTGCGCCTACGTCGGGAGGGGCAAAATACCGCTAAACAGCTTCCTGAAATTATCCTGCTTa                                                     >  1:19129264/1‑68 (MQ=25)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:15817011/1‑68 (MQ=255)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:425119/1‑68 (MQ=255)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:13967956/1‑68 (MQ=255)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:26497361/1‑68 (MQ=255)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:20760862/1‑68 (MQ=255)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:9702951/1‑68 (MQ=255)
                 cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                    >  1:2030406/1‑68 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:30909816/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:18866903/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:25305116/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:26809866/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:32529350/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:28596817/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:11653901/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:17575252/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:5073050/1‑67 (MQ=255)
                 cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATTATCCTGCTTa                                                     >  1:4745800/1‑67 (MQ=255)
                  atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:33546426/1‑66 (MQ=255)
                  atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:24655138/1‑66 (MQ=255)
                  atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                     >  1:1959787/1‑66 (MQ=255)
                   tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                                  <  1:27994078/68‑1 (MQ=255)
                    atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                                 <  1:32811250/68‑1 (MQ=255)
                     tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                                 >  1:24387524/1‑67 (MQ=255)
                     tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                                 >  1:29015933/1‑67 (MQ=255)
                     tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                                 >  1:31333213/1‑67 (MQ=255)
                     tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                                 >  1:33064135/1‑67 (MQ=255)
                      gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:30062130/1‑68 (MQ=255)
                      gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:34377669/1‑68 (MQ=255)
                      gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:29731348/1‑68 (MQ=255)
                      gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:23001747/1‑68 (MQ=255)
                      gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:17185550/1‑68 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:6684042/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:9901523/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:7428873/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:9378806/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:16534433/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:11121851/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:22761565/1‑67 (MQ=255)
                       ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                               >  1:12308943/1‑67 (MQ=255)
                        tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                             >  1:21516314/1‑68 (MQ=255)
                        tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                             >  1:30775073/1‑68 (MQ=255)
                           gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                          >  1:19370414/1‑68 (MQ=255)
                              tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                       <  1:27436333/68‑1 (MQ=255)
                              tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                       <  1:2407107/68‑1 (MQ=255)
                              tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                       <  1:8403013/68‑1 (MQ=255)
                              tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                       <  1:27263505/68‑1 (MQ=255)
                               aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                      <  1:26370762/68‑1 (MQ=255)
                               aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                      <  1:10480971/68‑1 (MQ=255)
                               aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                      <  1:13320837/68‑1 (MQ=255)
                                   cGGGAGGTGCAAATTACATGTAAACAGGTTCCGGAAATTATCATGCTTAACTCTCACGATGGCACCAg                                  <  1:13719214/68‑1 (MQ=21)
                                    gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                                 <  1:32285148/68‑1 (MQ=255)
                                    gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                                 <  1:5688436/68‑1 (MQ=255)
                                    gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                                 <  1:5247479/68‑1 (MQ=255)
                                     ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTAAATCTCACGATGGCACCAGtt                                >  1:25077349/1‑68 (MQ=255)
                                      gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                               >  1:30392864/1‑68 (MQ=255)
                                       aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                               <  1:1566629/67‑1 (MQ=255)
                                       aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                              >  1:21151658/1‑68 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                             <  1:27671225/66‑1 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                             <  1:7171040/66‑1 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                             <  1:14344333/66‑1 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                             <  1:3735765/66‑1 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                             <  1:14821978/66‑1 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                            <  1:443155/67‑1 (MQ=255)
                                          ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                            <  1:16831762/67‑1 (MQ=255)
                                           gCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTAcc                          >  1:21883385/1‑68 (MQ=255)
                                                  aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCATAtgct                   <  1:16112533/68‑1 (MQ=255)
                                                  aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct                   <  1:19958454/68‑1 (MQ=255)
                                                  aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct                   <  1:23832446/68‑1 (MQ=255)
                                                   ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg                  <  1:23380454/68‑1 (MQ=255)
                                                     ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc                >  1:10708959/1‑68 (MQ=255)
                                                        aaaCAGGTTCGGGAATTTATCCTGCTTAACTCTGACGATGTCACCAGTACGTACCAGATGCTGCCggt             <  1:7758574/68‑2 (MQ=255)
                                                        aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTGACGATGGCACCAGTTCGTACCAGATGCTGCCggg             <  1:11360962/68‑1 (MQ=255)
                                                        aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg             <  1:6491351/68‑1 (MQ=255)
                                                        aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg             <  1:2698242/68‑1 (MQ=255)
                                                            aGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATg         >  1:6889362/1‑68 (MQ=255)
                                                             ggTTCCGGAAATTATCCTGCTTAACTCTCCCGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt        >  1:8937367/1‑68 (MQ=255)
                                                             ggTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt        >  1:28393473/1‑68 (MQ=255)
                                                               ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt        <  1:13057104/66‑1 (MQ=255)
                                                               ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt        <  1:25844622/66‑1 (MQ=255)
                                                                   |                                                                   
GTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGG  >  CP009273/2762253‑2762387

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