Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,288 C→A 6.8% R89R (CGG→AGG)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2880CA 6.8% 1931.7 / 93.6 555R89R (CGG→AGG) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (345/172);  new base A (25/13);  total (370/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

Only 100 of 566 total aligned reads displayed.

GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355
                                                                  |                                                                   
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:7785961/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:6404774/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:4283538/1‑68 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                              <  1:27907973/68‑1 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                               <  1:11673722/67‑1 (MQ=255)
          cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                            <  1:29641963/66‑1 (MQ=255)
          cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                            <  1:20787619/66‑1 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:31511389/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:8315227/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         <  1:8936032/68‑1 (MQ=255)
               gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                                      >  1:26126222/1‑67 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                    <  1:13117900/67‑1 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                                   >  1:281990/1‑68 (MQ=255)
                  cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                  <  1:16205691/68‑1 (MQ=255)
                   gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                 >  1:15993666/1‑68 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:29386044/1‑68 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:32183421/1‑68 (MQ=255)
                     gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:15814762/1‑68 (MQ=255)
                     gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:20177580/1‑68 (MQ=255)
                            gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                        <  1:34183946/68‑1 (MQ=255)
                             tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                        >  1:18730338/1‑67 (MQ=255)
                                 cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                   >  1:19081080/1‑68 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:28535901/68‑1 (MQ=255)
                                   tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                 >  1:30975204/1‑68 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                              >  1:26871120/1‑68 (MQ=255)
                                         tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:31835362/1‑66 (MQ=255)
                                         tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:4210264/1‑66 (MQ=255)
                                          aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          <  1:3253376/68‑1 (MQ=255)
                                          aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          <  1:25116067/68‑1 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:10054902/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:11371619/1‑67 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:23598357/1‑66 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:14790800/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:33697976/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:14182409/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:5937424/68‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:18228964/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:12825980/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                       >  1:4896892/4‑68 (MQ=255)
                                              aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                      >  1:125985/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:19887705/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:26497361/1‑68 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGGTTa                     >  1:13674921/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:27734103/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:12152470/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:6982336/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:2701474/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:17680959/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:26372856/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:5073050/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:21809673/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:13221960/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGGTCCGGAAATTATCCTGCTTa                     >  1:3722537/1‑67 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGTTAAACAGGTTCCGGAAATTATCCAGCTTa                     >  1:10852022/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:19507268/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:20387605/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:15558156/68‑1 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:27434957/68‑1 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:25970949/68‑1 (MQ=255)
                                                  tatGCTGCGCCTACGTCTGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                  <  1:3593854/68‑1 (MQ=255)
                                                  tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                  <  1:28607520/68‑1 (MQ=255)
                                                  tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   >  1:5740525/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:7097306/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:6627186/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:33420592/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:33331234/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:19209931/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:15229174/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTACCGGAAATTATCCTGCTTAActc                 >  1:24312170/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGCCAAATTACCGGAAAACAGGTTCCGGAAAATATCCTGCTTAActct                >  1:11824256/1‑68 (MQ=255)
                                                    tgctgcACCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:978422/1‑67 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcta               >  1:29143494/1‑67 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:24719885/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:21126739/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:18406827/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:22061508/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:16540241/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:1649207/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:30801940/1‑68 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:11121851/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:25442175/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:7428873/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:23938276/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:4450746/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:12441780/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:9378806/1‑67 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc             >  1:29813136/1‑68 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc             >  1:30136343/1‑68 (MQ=255)
                                                           ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         >  1:14456052/1‑68 (MQ=255)
                                                           ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         >  1:22444155/1‑68 (MQ=255)
                                                            cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         <  1:5455293/67‑1 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:21574119/68‑1 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:22928403/68‑1 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:3930969/68‑1 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:33905242/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:10480971/68‑1 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  <  1:17127602/68‑1 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  <  1:34404720/68‑1 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  <  1:30610741/68‑1 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  <  1:32981577/68‑1 (MQ=255)
                                                                  |                                                                   
GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355

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: