Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,293 G→A 8.2% E90E (GAG→GAA yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2930GA 8.2% 2009.2 / 130.3 648E90E (GAG→GAAyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (432/163);  new base A (38/15);  total (470/178)
Fisher's exact test for biased strand distribution p-value = 8.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 673 total aligned reads displayed.

TCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGT  >  CP009273/2762227‑2762360
                                                                  |                                                                   
 cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                                   <  1:32171344/68‑1 (MQ=255)
 cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                                    <  1:19767802/67‑1 (MQ=255)
     cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                                 <  1:9099362/66‑1 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              >  1:34935117/1‑68 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              >  1:15762944/1‑68 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              >  1:15207833/1‑68 (MQ=255)
             cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                       <  1:5258266/68‑1 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTTAGCATATCCTGCGCCTACGTCGGGAGGGCCAAAATACCGGt                                                      >  1:20118287/1‑68 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      >  1:26310230/1‑68 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      >  1:5762144/1‑68 (MQ=255)
              gTGACCCCCGGCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      <  1:2779193/68‑1 (MQ=255)
                gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:35400497/1‑68 (MQ=255)
                gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:17826409/1‑68 (MQ=255)
                 aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:18160099/1‑67 (MQ=255)
                 aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:15987214/1‑67 (MQ=255)
                     ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                               <  1:35191997/68‑1 (MQ=255)
                        tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                             >  1:13170270/1‑67 (MQ=255)
                         cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                           >  1:15492339/1‑68 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:30149167/68‑1 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:32656537/68‑1 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:34714029/68‑1 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:8120639/68‑1 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:36729434/68‑1 (MQ=255)
                                 gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                   >  1:18940874/1‑68 (MQ=255)
                                    tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                  >  1:5978048/1‑66 (MQ=255)
                                     aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                  >  1:11611314/1‑65 (MQ=255)
                                      cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                               >  1:10123082/1‑67 (MQ=255)
                                      cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                               >  1:23684233/1‑67 (MQ=255)
                                      cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                              <  1:17212806/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:12766302/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:25412480/68‑1 (MQ=255)
                                        aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:10872406/64‑1 (MQ=255)
                                        aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:35842517/64‑1 (MQ=255)
                                         aGCATATGCTTCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                           >  1:12386781/1‑68 (MQ=255)
                                         aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                           >  1:38647340/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:7061655/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:28972316/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:31670131/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:16904944/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:34112183/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa                          >  1:18432790/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:21753957/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:23952795/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:11459902/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:4063398/1‑68 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:38131500/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:27237331/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:873338/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:28409686/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:37260337/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:19604682/1‑67 (MQ=255)
                                            atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                        <  1:4748301/68‑1 (MQ=255)
                                             tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                       <  1:21014522/68‑1 (MQ=255)
                                              atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      <  1:26505972/68‑1 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct                     <  1:8536071/68‑1 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:16280174/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:11986069/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:21579813/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:32420721/1‑67 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:9627290/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:4389475/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:38355018/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:19308863/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:21254811/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:13606898/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:29536983/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:24685584/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:33064941/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:837377/1‑68 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:3607579/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:17467332/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:29298502/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:25024754/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:30494063/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:14635447/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:30996047/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:24357782/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:7693984/1‑67 (MQ=255)
                                                   gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                 <  1:32038288/68‑1 (MQ=255)
                                                   gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                 <  1:17964910/68‑1 (MQ=255)
                                                   gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                 <  1:2255698/68‑1 (MQ=255)
                                                     gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                <  1:14032474/67‑1 (MQ=255)
                                                        tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc            <  1:11171146/68‑1 (MQ=255)
                                                        tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc            <  1:26815094/68‑1 (MQ=255)
                                                        tACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc            <  1:6493028/68‑1 (MQ=255)
                                                          cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc           >  1:10014583/1‑66 (MQ=255)
                                                            tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa        >  1:23120303/1‑68 (MQ=255)
                                                             cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg       <  1:37781716/68‑1 (MQ=255)
                                                             cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg       <  1:13000159/68‑1 (MQ=255)
                                                              tggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:20608368/67‑1 (MQ=255)
                                                              gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:38817270/68‑1 (MQ=255)
                                                              gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:7428150/68‑1 (MQ=255)
                                                              gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:34264651/68‑1 (MQ=255)
                                                               ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      >  1:39451025/1‑67 (MQ=255)
                                                               ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      >  1:27959606/1‑67 (MQ=255)
                                                                 aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg   >  1:1501323/1‑68 (MQ=255)
                                                                 aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg   >  1:36408732/1‑68 (MQ=255)
                                                                 aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg   >  1:33775378/1‑68 (MQ=255)
                                                                  ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt  >  1:22279399/1‑68 (MQ=255)
                                                                  ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt  >  1:25791867/1‑68 (MQ=255)
                                                                  |                                                                   
TCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGT  >  CP009273/2762227‑2762360

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: