Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,308 T→G 6.4% G95G (GGT→GGG)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3080TG 6.4% 2110.9 / 85.8 595G95G (GGT→GGG) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (328/227);  new base G (26/13);  total (355/240)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 602 total aligned reads displayed.

GTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCC  >  CP009273/2762241‑2762373
                                                                   |                                                                 
 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                  >  1:31843501/1‑68 (MQ=255)
 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                  >  1:3862680/1‑68 (MQ=255)
 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                  >  1:28249535/1‑68 (MQ=255)
  gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                 >  1:19621949/1‑68 (MQ=255)
    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                                                <  1:27849462/67‑1 (MQ=255)
     ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTTCGTCGGGAGGGGCTAATTACCGGGAATCa                                                              >  1:12099810/1‑68 (MQ=25)
           cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                                        >  1:18458010/1‑68 (MQ=255)
               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                    <  1:15817757/68‑1 (MQ=255)
               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                    <  1:8880685/68‑1 (MQ=255)
               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                    <  1:27627530/68‑1 (MQ=255)
                tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                                   >  1:23328500/1‑68 (MQ=255)
                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                >  1:34338447/1‑68 (MQ=255)
                     aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                               <  1:5820774/67‑1 (MQ=255)
                      tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                               >  1:16528053/1‑66 (MQ=255)
                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                            <  1:3253376/68‑1 (MQ=255)
                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                            >  1:10054902/1‑67 (MQ=255)
                         tAAGCTTATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                          <  1:22183774/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                          <  1:20440319/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                          <  1:1289186/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                          <  1:2920553/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                          <  1:2949069/68‑1 (MQ=255)
                         tAAACATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                          <  1:12648643/68‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                          <  1:20230958/64‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                         >  1:4805744/4‑68 (MQ=255)
                          aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                         >  1:23774747/1‑68 (MQ=255)
                          aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                         >  1:2651166/1‑68 (MQ=255)
                           aGCATATGCTGCGCCTACGTCGGGAAGTGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                        >  1:14572712/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:16375650/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:9493130/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAAATACCGCTAAACAGCTTCCTGAAATTATCCTGCTTa                                       >  1:19129264/1‑68 (MQ=25)
                            acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                       >  1:22730937/2‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                      >  1:21078742/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                      >  1:11259431/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                      <  1:7816833/68‑1 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:10982897/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:5073050/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:13167253/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:25305116/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:18929525/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:18252402/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:6982336/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                       >  1:28536935/1‑67 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                     <  1:25970949/68‑1 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                     <  1:24761432/68‑1 (MQ=255)
                                atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   <  1:32811250/68‑1 (MQ=255)
                                atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   <  1:15629437/68‑1 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:4253471/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:29015933/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:15229174/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:24387524/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:31333213/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:9818991/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:6627186/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                   >  1:7097306/1‑67 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATc                                              >  1:30875966/1‑55 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:5679865/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:26918801/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:4191224/1‑68 (MQ=255)
                                  gctgcGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTGCCGGAAATTATCCTGCTTAACtctc                                 >  1:33377422/1‑68 (MQ=18)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTTACtctc                                 >  1:3642915/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:5369383/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:12441780/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:7428873/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:11760219/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:21332596/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                 >  1:1990542/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCATGCTTAACtctc                                 >  1:17744953/1‑67 (MQ=255)
                                    tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                               >  1:29813136/1‑68 (MQ=255)
                                    tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                               >  1:30775073/1‑68 (MQ=255)
                                      cgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                             >  1:6158715/1‑68 (MQ=255)
                                       gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                             <  1:33511931/67‑1 (MQ=255)
                                        ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                           >  1:21598623/1‑68 (MQ=255)
                                         cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg                          <  1:34016511/68‑1 (MQ=255)
                                          tACTTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                         <  1:1162213/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                         <  1:22928403/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                         <  1:27436333/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                         <  1:33744171/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                        <  1:10480971/68‑1 (MQ=255)
                                               cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                    <  1:8406347/68‑1 (MQ=255)
                                                tggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                   <  1:23963657/67‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                   <  1:34466250/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                   <  1:32285148/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                   <  1:3008385/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                   <  1:14049069/68‑1 (MQ=255)
                                                 ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTAAATCTCACGATGGCACCAGtt                  >  1:25077349/1‑68 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt               <  1:30554058/66‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt               <  1:13360132/66‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt               <  1:33122278/66‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt               <  1:14344333/66‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt               <  1:14821978/66‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa              <  1:443155/67‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa              <  1:16831762/67‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGGTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt               <  1:26434706/66‑1 (MQ=255)
                                                       gCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTAcc            >  1:21883385/1‑68 (MQ=255)
                                                              aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct     <  1:17453563/68‑1 (MQ=255)
                                                              aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt        >  1:16041067/1‑65 (MQ=255)
                                                                 ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc  <  1:19370726/68‑1 (MQ=255)
                                                                 ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc  >  1:26156959/1‑68 (MQ=255)
                                                                 ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc  <  1:27184351/68‑1 (MQ=255)
                                                                 ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc  >  1:10708959/1‑68 (MQ=255)
                                                                   |                                                                 
GTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCC  >  CP009273/2762241‑2762373

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: ‑