Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,299 A→T 10.1% Q92H (CAA→CAT yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2990AT10.1% 1879.9 / 152.9 610Q92H (CAA→CATyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (329/219);  new base T (55/7);  total (384/226)
Fisher's exact test for biased strand distribution p-value = 2.88e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 638 total aligned reads displayed.

ACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGA  >  CP009273/2762233‑2762366
                                                                  |                                                                   
acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                                    >  1:26299306/1‑68 (MQ=255)
      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                               <  1:20817694/67‑1 (MQ=255)
      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                               <  1:28203190/67‑1 (MQ=255)
      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                                              >  1:15306778/1‑68 (MQ=255)
        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGGAAATTACCGGt                                                            >  1:11078389/1‑68 (MQ=255)
         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                           >  1:6746867/1‑68 (MQ=255)
           aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                          >  1:6522539/1‑67 (MQ=255)
            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                        >  1:19743859/1‑68 (MQ=255)
              cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                      >  1:6076309/1‑68 (MQ=255)
              cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                      >  1:12250712/1‑68 (MQ=255)
                       gtGAGCATACTAAGCATATGCTGCGCCTCCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                             <  1:21582794/68‑1 (MQ=255)
                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                             <  1:26396327/68‑1 (MQ=255)
                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                             <  1:4106869/68‑1 (MQ=255)
                        tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                            >  1:24341604/1‑68 (MQ=255)
                           gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAttt                                         >  1:8820723/1‑67 (MQ=255)
                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTAAACAGGTTCCGGAAATTATCCt                                     <  1:30015463/68‑1 (MQ=255)
                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                        >  1:31240921/1‑65 (MQ=255)
                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                     <  1:34325357/68‑1 (MQ=255)
                                cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                    <  1:937904/68‑1 (MQ=255)
                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:22161951/68‑1 (MQ=255)
                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:185853/68‑1 (MQ=255)
                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:28619170/68‑1 (MQ=255)
                                  aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                  >  1:10303445/1‑68 (MQ=255)
                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:34808626/1‑68 (MQ=255)
                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:27074642/1‑68 (MQ=255)
                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:12696702/1‑68 (MQ=255)
                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:10034764/1‑68 (MQ=255)
                                    acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCACGAAGTTCCTGAAATTATCCTGCTTa                                >  1:8212610/2‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:23690893/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:32056291/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:31771890/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:20124569/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:20514560/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:7658819/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:19468235/1‑68 (MQ=255)
                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:17193476/1‑68 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa                               >  1:33949800/1‑68 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:15553206/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:5405274/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:18470683/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:16690053/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:32431648/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:17849792/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:17254874/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:24825059/1‑67 (MQ=255)
                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:16921040/1‑67 (MQ=255)
                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:29114777/1‑66 (MQ=255)
                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:22625242/1‑66 (MQ=255)
                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                              <  1:4642141/68‑1 (MQ=255)
                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                              <  1:28884586/68‑1 (MQ=255)
                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                             <  1:10850017/68‑1 (MQ=255)
                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            >  1:22859665/1‑67 (MQ=255)
                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            >  1:21373711/1‑67 (MQ=255)
                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            >  1:13659459/1‑67 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:11793163/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:8671199/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:29620160/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:30667107/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:27623004/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:26733988/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:7208248/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:15996758/1‑68 (MQ=255)
                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:18059418/1‑68 (MQ=255)
                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcgc                          >  1:20975141/1‑65 (MQ=255)
                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:9800992/1‑67 (MQ=255)
                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:11295086/1‑67 (MQ=255)
                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:14234876/1‑67 (MQ=255)
                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:14669817/1‑67 (MQ=255)
                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCATAACtctc                          >  1:13044381/1‑67 (MQ=255)
                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                         <  1:1502479/67‑1 (MQ=255)
                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                        >  1:3030181/1‑68 (MQ=255)
                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                        >  1:19029407/1‑68 (MQ=255)
                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                        >  1:1762603/1‑68 (MQ=255)
                                            tgcGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                        >  1:25637414/1‑68 (MQ=255)
                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                      <  1:13816718/67‑1 (MQ=255)
                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                      <  1:25082965/67‑1 (MQ=255)
                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                     >  1:11542573/1‑68 (MQ=255)
                                                 cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                    >  1:33335086/1‑67 (MQ=255)
                                                 cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                    <  1:23180613/67‑1 (MQ=255)
                                                  tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                  <  1:16426284/68‑1 (MQ=255)
                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 <  1:27899363/68‑1 (MQ=255)
                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 <  1:21921880/68‑1 (MQ=255)
                                                    cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 >  1:2838349/1‑67 (MQ=255)
                                                      tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa              >  1:506796/1‑68 (MQ=255)
                                                      tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa              >  1:7970673/1‑68 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:10569488/68‑1 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:20179952/68‑1 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:29345025/68‑1 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:30947394/68‑1 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:32720947/68‑1 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:33528953/68‑1 (MQ=255)
                                                        gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:5713810/68‑1 (MQ=255)
                                                           aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg         >  1:31410024/1‑68 (MQ=255)
                                                           aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg         >  1:26847937/1‑68 (MQ=255)
                                                           aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg         >  1:3756529/1‑68 (MQ=255)
                                                              ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt        <  1:26079664/66‑1 (MQ=255)
                                                              ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt        <  1:26645071/66‑1 (MQ=255)
                                                              ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt        <  1:9117156/66‑1 (MQ=255)
                                                              ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa       <  1:6350563/67‑1 (MQ=255)
                                                              ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa       <  1:25451937/67‑1 (MQ=255)
                                                                  |                                                                   
ACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGA  >  CP009273/2762233‑2762366

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: