Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,293 G→A 10.3% 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,2930GA10.3% 1786.7 / 159.8 604E90E (GAG→GAAyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (342/198);  new base A (55/7);  total (397/207)
Fisher's exact test for biased strand distribution p-value = 2.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 619 total aligned reads displayed.

TCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGT  >  CP009273/2762227‑2762360
                                                                  |                                                                   
  aGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                                   >  1:5641599/1‑67 (MQ=255)
     cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                                 <  1:712007/66‑1 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              >  1:11855398/1‑68 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              >  1:2761273/1‑68 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              >  1:2932678/1‑68 (MQ=255)
      acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                              <  1:28963075/68‑1 (MQ=255)
            aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                         <  1:21455720/67‑1 (MQ=255)
            aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                         <  1:28203190/67‑1 (MQ=255)
            aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                         <  1:15215867/67‑1 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGGAAATTACCGGt                                                      >  1:11078389/1‑68 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      <  1:26615006/68‑1 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      <  1:20619257/68‑1 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      >  1:17447716/1‑68 (MQ=255)
              gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      >  1:28776015/1‑68 (MQ=255)
               tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                     >  1:25265406/1‑68 (MQ=255)
                gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:3363042/1‑68 (MQ=255)
                gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:3971947/1‑68 (MQ=255)
                  cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaatc                                                  >  1:34873354/1‑66 (MQ=255)
                        tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                             >  1:17210782/1‑67 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:13468686/68‑1 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:34525169/68‑1 (MQ=255)
                             gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:17726481/68‑1 (MQ=255)
                             gtGAGCATACAAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:11604433/68‑1 (MQ=255)
                                 gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                   >  1:4775688/1‑68 (MQ=255)
                                    tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                  >  1:2594949/1‑66 (MQ=255)
                                     aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                  >  1:30180890/1‑65 (MQ=255)
                                      cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                               >  1:18486239/1‑67 (MQ=255)
                                      cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                              <  1:937904/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                               >  1:1788949/1‑66 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:16098467/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:14772896/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:14346317/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:3185385/68‑1 (MQ=255)
                                       tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAATATCCt                               >  1:33170662/1‑66 (MQ=255)
                                        aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:20907736/64‑1 (MQ=255)
                                        aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                            >  1:29572198/1‑68 (MQ=255)
                                        aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                            >  1:31541674/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:10034764/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:20130125/1‑68 (MQ=255)
                                          gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:7819397/1‑68 (MQ=255)
                                          acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                          >  1:8947264/2‑68 (MQ=255)
                                          acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                          >  1:27887648/2‑68 (MQ=255)
                                          acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCACGAAGTTCCTGAAATTATCCTGCTTa                          >  1:5137769/2‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:16771914/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:32528368/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:30905236/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:12977463/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:270432/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:8126553/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:19571084/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:6671434/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:4188290/1‑68 (MQ=255)
                                           cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa                         >  1:15632952/1‑68 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTa                          >  1:16956021/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:33494169/1‑67 (MQ=255)
                                           cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:24078713/1‑67 (MQ=255)
                                            atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:18288468/1‑66 (MQ=255)
                                            atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:14172434/1‑66 (MQ=255)
                                            atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:21226390/1‑66 (MQ=255)
                                             tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                       <  1:34045003/68‑1 (MQ=255)
                                             tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                        >  1:1538417/1‑67 (MQ=255)
                                              atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      <  1:13700081/68‑1 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGTTTCAGGAAATTATCCTGCTTAActc                      >  1:20404366/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:1655718/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:53538/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:22767423/1‑67 (MQ=255)
                                               tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:22431676/1‑67 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:26733988/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:2978188/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:12680870/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:12185400/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:7452361/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:26284667/1‑68 (MQ=255)
                                                gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:3226367/1‑68 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:8759827/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:9800992/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:31203378/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:10718794/1‑67 (MQ=255)
                                                 ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:22042334/1‑67 (MQ=255)
                                                  tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                   <  1:10902952/67‑1 (MQ=255)
                                                  tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                   <  1:19845283/67‑1 (MQ=255)
                                                     gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt               >  1:8420731/1‑68 (MQ=255)
                                                      ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg              >  1:21688594/1‑68 (MQ=255)
                                                         aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa           <  1:10311697/68‑1 (MQ=255)
                                                         aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa           <  1:30599487/68‑1 (MQ=255)
                                                          cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa           >  1:2838349/1‑67 (MQ=255)
                                                          cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa           >  1:11301362/1‑67 (MQ=255)
                                                            tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa        >  1:24646520/1‑68 (MQ=255)
                                                            tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa        >  1:31354006/1‑68 (MQ=255)
                                                             cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg       <  1:6407376/68‑1 (MQ=255)
                                                              gggATGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:25518389/68‑1 (MQ=255)
                                                              gggAGGGGCAATTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:18671489/68‑1 (MQ=255)
                                                              gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:6003980/68‑1 (MQ=255)
                                                              gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:28594722/68‑1 (MQ=255)
                                                              gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTACCCTTCTTAACTCTCACGATGGCACCAGt      <  1:24929233/68‑1 (MQ=255)
                                                                gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc    >  1:2342197/1‑68 (MQ=255)
                                                                gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc    >  1:23115721/1‑68 (MQ=255)
                                                                gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc    >  1:26394414/1‑68 (MQ=255)
                                                                gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc    >  1:1911712/1‑68 (MQ=255)
                                                                 aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg   >  1:26847937/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: