Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,314 G→A 10.3% Q97Q (CAG→CAA) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3140GA10.3% 1805.3 / 157.2 606Q97Q (CAG→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (293/248);  new base A (55/7);  total (350/256)
Fisher's exact test for biased strand distribution p-value = 4.36e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 618 total aligned reads displayed.

CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381
                                                                   |                                                                   
cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                     >  1:7870995/1‑68 (MQ=255)
   gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                                                   >  1:19822110/1‑67 (MQ=255)
    tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                                  >  1:21605357/1‑67 (MQ=255)
        cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                             >  1:11894630/1‑68 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:26396327/68‑1 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:25997970/68‑1 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:30780246/68‑1 (MQ=255)
         gtGAGCATACAAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:11604433/68‑1 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                        >  1:16796289/1‑68 (MQ=255)
              cATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                        >  1:23235968/1‑67 (MQ=255)
               aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       <  1:20947479/67‑1 (MQ=255)
                 aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:28690453/1‑65 (MQ=255)
                 aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    <  1:31173214/68‑1 (MQ=255)
                  cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                                   <  1:34441504/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:2518434/1‑66 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:13479558/68‑1 (MQ=255)
                    aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                                 >  1:10303445/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:10034764/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:2065771/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:32385254/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:27074642/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:12696702/1‑68 (MQ=255)
                      acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                               >  1:71282/2‑68 (MQ=255)
                      acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                               >  1:26717805/2‑68 (MQ=255)
                      acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                               >  1:4017149/2‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:29343515/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:13127551/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:25073811/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:8126553/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:14206347/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:16362126/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:1059315/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:19468235/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa                                              >  1:15632952/1‑68 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATGCTGCTTa                                               >  1:33993666/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:18678621/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:18470683/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:22740350/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:28941115/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:24078713/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTa                                               >  1:31324739/1‑67 (MQ=255)
                       cATATGCAGCGCCTACGTCGGGAGGCGCAAAATACCGGTAAACACGTTCCGCAAATTATACTGCATa                                               >  1:8476866/1‑67 (MQ=21)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:15175767/1‑66 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:30268990/1‑66 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             <  1:14583762/68‑1 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                            <  1:10850017/68‑1 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:626367/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:22281078/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:21373711/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:5406577/1‑67 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:29620160/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:12185400/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:18059418/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:27623004/1‑68 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:8759827/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:7503084/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:33541761/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:25591834/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:11295086/1‑67 (MQ=255)
                              tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                       >  1:2991751/1‑68 (MQ=255)
                              tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                       >  1:1762603/1‑68 (MQ=255)
                                 gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                    >  1:28437337/1‑68 (MQ=255)
                                 gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAACTATCCTGCTTAACTCTCACGa                                     <  1:24897543/67‑1 (MQ=255)
                                  ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   >  1:645015/1‑68 (MQ=255)
                                  ccTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   >  1:9649409/1‑68 (MQ=255)
                                   cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   >  1:33335086/1‑67 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:28257133/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:4451475/68‑1 (MQ=255)
                                     aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                <  1:27899363/68‑1 (MQ=255)
                                     aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                <  1:23057966/68‑1 (MQ=255)
                                      cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                >  1:31480127/1‑67 (MQ=255)
                                        tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa                             >  1:24646520/1‑68 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:22022041/68‑1 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:9305150/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:11094361/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:32720947/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:15542069/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:20128413/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:6003980/68‑1 (MQ=255)
                                            gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                         >  1:1911712/1‑68 (MQ=255)
                                             aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                        >  1:26847937/1‑68 (MQ=255)
                                             aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                        >  1:17239796/1‑68 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:9117156/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:26645071/66‑1 (MQ=255)
                                                    aaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGa                 <  1:30550299/68‑1 (MQ=255)
                                                      ttACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt                >  1:23659705/1‑67 (MQ=255)
                                                        aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct             <  1:3530276/68‑1 (MQ=255)
                                                        aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct             <  1:21034899/68‑1 (MQ=255)
                                                           ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc          >  1:17085037/1‑68 (MQ=255)
                                                              aaaCAGGTTCCGTAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:31911466/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:26287798/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:14797512/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:13746733/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:524200/68‑1 (MQ=255)
                                                              aaaCAGGTTACGGAATTTATGCTGCTAAACTCTCACGATGGCACCAGTTCGTTCCAGATGCTGCCggt       <  1:34655972/68‑2 (MQ=18)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCCCGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     >  1:504614/1‑68 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     >  1:16616531/1‑68 (MQ=255)
                                                                  aGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATg   >  1:6883102/1‑68 (MQ=255)
                                                                  aGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATg   >  1:20334776/1‑68 (MQ=255)
                                                                   ggTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCTCCAGTTCGTACCAGATGCTGCCGGGCATGt  >  1:178723/1‑68 (MQ=255)
                                                                   |                                                                   
CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381

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: