Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,314 G→A 6.2% 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,3140GA 6.2% 1941.1 / 84.7 598Q97Q (CAG→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (310/251);  new base A (24/13);  total (334/264)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 615 total aligned reads displayed.

CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381
                                                                   |                                                                   
        cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                             >  1:26104866/1‑68 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:28298447/68‑1 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:27627530/68‑1 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                        >  1:34338447/1‑68 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                        >  1:12077266/1‑68 (MQ=255)
               aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       <  1:5820774/67‑1 (MQ=255)
               aTACTAAGCATATGCTGCGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       <  1:7035763/67‑1 (MQ=255)
                tACTAAGCATATGCTGCGCCTACGTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:4359927/1‑66 (MQ=255)
                 aCTACGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    <  1:3984851/68‑1 (MQ=255)
                  cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:10054902/1‑67 (MQ=255)
                  cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:30072947/1‑67 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:18276414/1‑66 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:10464644/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:2949069/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:21641196/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:15728144/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:14790800/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:10968466/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:11697434/68‑1 (MQ=255)
                    aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                  <  1:12825980/64‑1 (MQ=255)
                     aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                                >  1:3718783/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:9657380/1‑68 (MQ=255)
                      acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                               <  1:15986431/67‑1 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:1976419/1‑68 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:5119104/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:28536935/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:30909816/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:18494544/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:11214640/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:10696727/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:16161631/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:25765953/1‑67 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:22893261/1‑66 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:20387605/1‑66 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             <  1:15596844/68‑1 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             <  1:13354688/68‑1 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             <  1:15177738/68‑1 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             <  1:23173907/68‑1 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                            <  1:9856465/68‑1 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             >  1:24901544/1‑67 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             >  1:25179442/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            >  1:4234062/1‑50 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:33331234/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:31333213/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:20950094/1‑67 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGTTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:28793012/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:1649207/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:14265670/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:33492374/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:30780470/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:29185305/1‑68 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:14515966/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:6212445/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:32374780/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:1380331/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:23905936/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCATGCTTAACtctc                                         >  1:17744953/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGCAGGCGCTAAATACCGGTAATCAGGTTCCGGAATTTATCCTGCTTAACtctc                                         >  1:27226229/1‑67 (MQ=255)
                              tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                                        <  1:5657826/67‑1 (MQ=255)
                              tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                       >  1:7863666/1‑68 (MQ=255)
                              tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                       >  1:21516314/1‑68 (MQ=255)
                               gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                      <  1:26476535/68‑1 (MQ=255)
                                 gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                    >  1:23683015/1‑68 (MQ=255)
                                 gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                    >  1:12414456/1‑68 (MQ=255)
                                  ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   >  1:22444155/1‑68 (MQ=255)
                                   cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg                                  <  1:34016511/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:4932764/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:27436333/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:33744171/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:6568937/68‑1 (MQ=255)
                                     aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                <  1:26370762/68‑1 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:3193314/68‑1 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:32981577/68‑1 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:18932460/68‑1 (MQ=255)
                                          ggtaGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:2428171/65‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:4550388/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:11511646/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:34466250/68‑1 (MQ=255)
                                            gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                         >  1:16617579/1‑68 (MQ=255)
                                            gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                         >  1:30392864/1‑68 (MQ=255)
                                             aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                        >  1:8589226/1‑68 (MQ=255)
                                              ggggCAAATTACCGGTAAACAGGTTCCTGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       >  1:19965857/1‑68 (MQ=255)
                                                  cAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTAcc                    <  1:32676190/67‑1 (MQ=255)
                                                     aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCCCGATGGCACCAGGTAGTACCa                   >  1:21971185/1‑65 (MQ=21)
                                                     aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt                <  1:12648202/68‑1 (MQ=255)
                                                       tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTATCTCTCACGATGGCACCAGTTCGTACCAGAt                >  1:21243731/1‑66 (MQ=255)
                                                       tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt                >  1:13061071/1‑66 (MQ=255)
                                                        aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct             <  1:17453563/68‑1 (MQ=255)
                                                        aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct             <  1:19151095/68‑1 (MQ=255)
                                                           ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc          <  1:19370726/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:26957369/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:29788793/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:7317445/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:7652890/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:24494281/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:9231471/68‑1 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCTGATGCTGCCGGGCa     >  1:27835895/1‑68 (MQ=255)
                                                                  aGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATg   >  1:6889362/1‑68 (MQ=255)
                                                                   ggTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt  >  1:26806292/1‑68 (MQ=255)
                                                                   ggTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt  >  1:1688092/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: