Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,293 G→A 6.7% 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 6.7% 1858.3 / 85.9 579E90E (GAG→GAAyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (344/196);  new base A (26/13);  total (370/209)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 600 total aligned reads displayed.

GTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGT  >  CP009273/2762226‑2762360
                                                                   |                                                                   
  cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                                   <  1:31941420/68‑1 (MQ=255)
     acacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCa                                                                >  1:31747262/1‑68 (MQ=255)
           gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                                           >  1:22069267/1‑67 (MQ=255)
             aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                         <  1:13117900/67‑1 (MQ=255)
             aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                                        >  1:1602480/1‑68 (MQ=255)
              cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                       >  1:26441977/1‑68 (MQ=255)
               gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      >  1:3708005/1‑68 (MQ=255)
               gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                      <  1:29795873/68‑1 (MQ=255)
                tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                     >  1:28249535/1‑68 (MQ=255)
                 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTTCGTCGGGAGGGGCAAATTACCGGTaa                                                    >  1:26082806/1‑68 (MQ=255)
                    ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTTCGTCGGGAGGGGCTAATTACCGGGAATCa                                                 >  1:12099810/1‑68 (MQ=25)
                      ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                               <  1:29123151/68‑1 (MQ=255)
                       cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                              <  1:29586359/68‑1 (MQ=255)
                        gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGt                                              >  1:7039732/1‑67 (MQ=255)
                          cgtcgtGCGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                           >  1:11084832/1‑68 (MQ=255)
                          cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                           <  1:11639183/68‑1 (MQ=255)
                              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:15817757/68‑1 (MQ=255)
                              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                       <  1:8880685/68‑1 (MQ=255)
                               tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                      >  1:19151195/1‑68 (MQ=255)
                                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                   >  1:29228381/1‑68 (MQ=255)
                                    aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                  <  1:5820774/67‑1 (MQ=255)
                                     tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                  >  1:16528053/1‑66 (MQ=255)
                                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                               >  1:10054902/1‑67 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                               >  1:18276414/1‑66 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:27467311/68‑1 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:20001750/68‑1 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:14182409/68‑1 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:22578419/68‑1 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:21744938/68‑1 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                             <  1:14790800/68‑1 (MQ=255)
                                        tAAGCATATGCTGCGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTTCGGAAATTATCCTGc                             <  1:33665442/68‑1 (MQ=21)
                                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:17195611/64‑1 (MQ=255)
                                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:12825980/64‑1 (MQ=255)
                                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:13231038/64‑1 (MQ=255)
                                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                             <  1:40736/64‑1 (MQ=255)
                                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                            >  1:4805744/4‑68 (MQ=255)
                                          aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGTAATTATACTGCtt                           >  1:34140997/1‑68 (MQ=255)
                                          aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                           >  1:125985/1‑68 (MQ=255)
                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:30754710/1‑68 (MQ=255)
                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:21078742/1‑68 (MQ=255)
                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:1131059/1‑68 (MQ=255)
                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:2030406/1‑68 (MQ=255)
                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         <  1:7816833/68‑1 (MQ=255)
                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                         >  1:1976419/1‑68 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:26920337/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:7420498/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:13766170/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:28596817/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:1638561/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:17757625/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:18494544/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:5073050/1‑67 (MQ=255)
                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                          >  1:25765953/1‑67 (MQ=255)
                                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa                         <  1:32700738/67‑1 (MQ=255)
                                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                       <  1:27994078/68‑1 (MQ=255)
                                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                        >  1:25179442/1‑67 (MQ=255)
                                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                        >  1:24901544/1‑67 (MQ=255)
                                               atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      <  1:15629437/68‑1 (MQ=255)
                                               atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      <  1:3863056/68‑1 (MQ=255)
                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:33331234/1‑67 (MQ=255)
                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:4253471/1‑67 (MQ=255)
                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:6918706/1‑67 (MQ=255)
                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:15229174/1‑67 (MQ=255)
                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:31333213/1‑67 (MQ=255)
                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTACCGGAAATTATCCTGCTTAActc                      >  1:24312170/1‑67 (MQ=255)
                                                tgctgcACCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                      >  1:978422/1‑67 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATc                                 >  1:30875966/1‑55 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:30062130/1‑68 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:23001747/1‑68 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:33016740/1‑68 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:5679865/1‑68 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:34377669/1‑68 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:24719885/1‑68 (MQ=255)
                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:33818588/1‑68 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:26299323/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:4450746/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:5369383/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:6212445/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:11760219/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:32374780/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                    >  1:9378806/1‑67 (MQ=255)
                                                  ctgcGCCTACGTCGGCAGGCGCTAAATACCGGTAATCAGGTTCCGGAATTTATCCTGCTTAACtctc                    >  1:27226229/1‑67 (MQ=255)
                                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                   <  1:30365560/67‑1 (MQ=255)
                                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                  >  1:21516314/1‑68 (MQ=255)
                                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt               >  1:23683015/1‑68 (MQ=255)
                                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt               >  1:19370414/1‑68 (MQ=255)
                                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt               >  1:12414456/1‑68 (MQ=255)
                                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt               >  1:27742712/1‑68 (MQ=255)
                                                       ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg              >  1:14456052/1‑68 (MQ=255)
                                                        cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg              <  1:26572011/67‑1 (MQ=255)
                                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc            <  1:6568937/68‑1 (MQ=255)
                                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc            <  1:8403013/68‑1 (MQ=255)
                                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa           <  1:20708331/68‑1 (MQ=255)
                                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa           <  1:10480971/68‑1 (MQ=255)
                                                            gTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAcc         >  1:16848533/1‑68 (MQ=255)
                                                              cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg       <  1:18932460/68‑1 (MQ=255)
                                                               tggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:23963657/67‑1 (MQ=255)
                                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt      <  1:3553948/68‑1 (MQ=255)
                                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATAGCACCAGt      <  1:10752064/68‑1 (MQ=255)
                                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg   >  1:33446931/1‑68 (MQ=255)
                                                                   |                                                                   
GTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGT  >  CP009273/2762226‑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: