Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,309 A→C 10.6% K96Q (AAA→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,3090AC10.6% 1873.9 / 181.8 581K96Q (AAA→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (297/221);  new base C (55/7);  total (352/229)
Fisher's exact test for biased strand distribution p-value = 4.07e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 629 total aligned reads displayed.

TGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGG  >  CP009273/2762242‑2762376
                                                                   |                                                                   
tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                     >  1:6746867/1‑68 (MQ=255)
 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                    <  1:27934808/68‑1 (MQ=255)
 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                    >  1:15231855/1‑68 (MQ=255)
 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                    >  1:3170851/1‑68 (MQ=255)
   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                  >  1:5439384/1‑68 (MQ=255)
   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                  >  1:32765835/1‑68 (MQ=255)
   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                  >  1:3066911/1‑68 (MQ=255)
     cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                >  1:7870995/1‑68 (MQ=255)
     cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                >  1:31143727/1‑68 (MQ=255)
     cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                >  1:12250712/1‑68 (MQ=255)
      ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                               <  1:33395152/68‑1 (MQ=255)
             cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTATACAGGTTCCGGa                                                        >  1:30360564/1‑68 (MQ=255)
              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       <  1:25997970/68‑1 (MQ=255)
              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGCaa                                                       <  1:25088336/68‑1 (MQ=255)
              gtGAGCATACAAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       <  1:11604433/68‑1 (MQ=255)
               tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                                      >  1:24341604/1‑68 (MQ=255)
                     tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:30135815/1‑66 (MQ=255)
                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:28690453/1‑65 (MQ=255)
                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:31240921/1‑65 (MQ=255)
                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGAAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:26302821/1‑65 (MQ=255)
                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               >  1:18486239/1‑67 (MQ=255)
                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                              <  1:937904/68‑1 (MQ=255)
                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                              <  1:34441504/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:26856458/68‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                             <  1:2183009/64‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                             <  1:16922034/64‑1 (MQ=255)
                         aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                            >  1:10303445/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:12696702/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:10034764/1‑68 (MQ=255)
                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                          >  1:16625383/2‑68 (MQ=255)
                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                          <  1:19828719/67‑1 (MQ=255)
                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                          >  1:71282/2‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGGAATTATCCTGCTTaa                                         >  1:24820268/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:13127551/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:16362126/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:32056291/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:1059315/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:23690893/1‑68 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAATTTATCCTGCTTa                                          >  1:5253776/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATGCTGCTTa                                          >  1:33993666/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:27745126/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:32431648/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCCGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:15945767/1‑67 (MQ=255)
                            cATATGCTGCGCCAACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:20963869/1‑67 (MQ=255)
                            cATATGCAGCGCCTACGTCGGGAGGCGCAAAATACCGGTAAACACGTTCCGCAAATTATACTGCATa                                          >  1:8476866/1‑67 (MQ=21)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:29114777/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:18288468/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:21226390/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        <  1:28447898/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                       <  1:23101991/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                       <  1:10850017/68‑1 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:14940136/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:13659459/1‑67 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTGCGGAAATTATCCTGCTTAACtctc                                    >  1:19667873/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcac                                    >  1:17707702/1‑66 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:26733988/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:26649578/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:29672668/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:11110821/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTAACtctc                                    >  1:13787432/1‑68 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:8759827/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:7503084/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:9800992/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:21441790/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:14234876/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:14669817/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:22042334/1‑67 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                                   <  1:33629490/67‑1 (MQ=255)
                                   tgcGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                  >  1:25637414/1‑68 (MQ=255)
                                    gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                 <  1:23315365/68‑1 (MQ=255)
                                       ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                              >  1:645015/1‑68 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:4451475/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:20444958/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:171632/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:9160478/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:26454329/68‑1 (MQ=255)
                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           <  1:28275756/68‑1 (MQ=255)
                                           cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           >  1:11301362/1‑67 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:15542069/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:29345025/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:20690298/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:22744163/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTACCCTTCTTAACTCTCACGATGGCACCAGt                      <  1:24929233/68‑1 (MQ=255)
                                                 gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGtt                     <  1:8136359/67‑1 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                    <  1:20141790/67‑1 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                    <  1:4972355/67‑1 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                   >  1:18813903/1‑68 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                   >  1:31410024/1‑68 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:3990624/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:6011048/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:6324707/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:17251444/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                 <  1:25451937/67‑1 (MQ=255)
                                                             aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct        <  1:3530276/68‑1 (MQ=255)
                                                              ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg       <  1:27212296/68‑1 (MQ=255)
                                                              ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg       <  1:19205618/68‑1 (MQ=255)
                                                                ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc     <  1:34755340/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:22293120/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:17966945/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:11900425/68‑1 (MQ=255)
                                                                   |                                                                   
TGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGG  >  CP009273/2762242‑2762376

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: