Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,304 C→A 6.2% T94N (ACC→AAC)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3040CA 6.2% 2105.2 / 85.9 599T94N (ACC→AAC) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (335/225);  new base A (26/13);  total (361/238)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 612 total aligned reads displayed.

ACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTG  >  CP009273/2762239‑2762371
                                                                 |                                                                   
aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTGCGGCGGGAGGGGCAAATTAcc                                                                    <  1:33505395/67‑1 (MQ=255)
aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                                    <  1:13159758/67‑1 (MQ=255)
aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                                    <  1:34303730/67‑1 (MQ=255)
aCGTGACCCCCGTCGGCGTGAGCATACTAAGCATATGCTGCGCCTACGGCGGGAGGGGCAAATTAcc                                                                    <  1:12495836/67‑1 (MQ=25)
  gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                 >  1:30701170/1‑68 (MQ=255)
   tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                >  1:32183421/1‑68 (MQ=255)
    gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                               >  1:33343425/1‑68 (MQ=255)
    gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                               >  1:19621949/1‑68 (MQ=255)
    gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                               >  1:33741777/1‑68 (MQ=255)
     aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCAGTaa                                                               >  1:5820646/1‑67 (MQ=255)
      cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                                              <  1:22291424/67‑1 (MQ=255)
        cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAActg                                                           <  1:6066273/68‑3 (MQ=255)
         ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                          <  1:21406963/68‑1 (MQ=255)
           gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                        <  1:12028300/68‑1 (MQ=255)
             cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                                      >  1:18458010/1‑68 (MQ=255)
              gtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCg                                                     >  1:3592073/1‑68 (MQ=255)
                 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                  <  1:8880685/68‑1 (MQ=255)
                   gAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAt                                                <  1:11386292/68‑1 (MQ=255)
                     gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                              >  1:26871120/1‑68 (MQ=255)
                        tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                             >  1:31835362/1‑66 (MQ=255)
                         aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                          <  1:25116067/68‑1 (MQ=255)
                          cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                          >  1:10054902/1‑67 (MQ=255)
                           tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                        <  1:11697434/68‑1 (MQ=255)
                           tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                        <  1:31219621/68‑1 (MQ=255)
                            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                        <  1:17195611/64‑1 (MQ=255)
                            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                        <  1:20230958/64‑1 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:9493130/1‑68 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:15246836/1‑68 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAAATACCGCTAAACAGCTTCCTGAAATTATCCTGCTTa                                     >  1:19129264/1‑68 (MQ=25)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:13967956/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:26497361/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:423791/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    <  1:7816833/68‑1 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTTCTTa                                     >  1:28776239/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:6982336/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:18252402/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:4403429/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:18929525/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:32529350/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:15899599/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAATATCCTGCTTa                                     >  1:33971686/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATTATCCTGCTTa                                     >  1:4745800/1‑67 (MQ=255)
                                atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                   <  1:13354688/68‑1 (MQ=255)
                                atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                   <  1:27434957/68‑1 (MQ=255)
                                atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                   <  1:25970949/68‑1 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:9818991/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:7097306/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:21616138/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:4031490/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:33064135/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:16518271/1‑67 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:21126739/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:14265670/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               <  1:12292475/68‑1 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:29731348/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:30062130/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:30858016/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGATTAACtctc                               >  1:32772016/1‑68 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:12707709/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:11162926/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:25442175/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:14515966/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:22761565/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:1990542/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:7428873/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:16344092/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCTTAACtctc                               >  1:5065044/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCATGCTTAACtctc                               >  1:17744953/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAACTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:28118555/1‑67 (MQ=255)
                                      tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                              <  1:5657826/67‑1 (MQ=255)
                                      tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                              <  1:30365560/67‑1 (MQ=255)
                                         gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                          >  1:19370414/1‑68 (MQ=255)
                                            tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                       <  1:2407107/68‑1 (MQ=255)
                                            tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                       <  1:6568937/68‑1 (MQ=255)
                                             aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                      <  1:10480971/68‑1 (MQ=255)
                                             aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                      <  1:10772910/68‑1 (MQ=255)
                                             aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                      <  1:27800958/68‑1 (MQ=255)
                                             aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                      <  1:20708331/68‑1 (MQ=255)
                                                 cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                  <  1:8406347/68‑1 (MQ=255)
                                                 cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                  <  1:34404720/68‑1 (MQ=255)
                                                  gggTGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:24763238/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:5247479/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:2940509/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:2629487/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:27145915/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:27584460/68‑1 (MQ=255)
                                                    gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc               >  1:29162648/1‑68 (MQ=255)
                                                     aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc               <  1:1566629/67‑1 (MQ=255)
                                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt             <  1:24401234/66‑1 (MQ=255)
                                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt             <  1:3735765/66‑1 (MQ=255)
                                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt             <  1:14821978/66‑1 (MQ=255)
                                                        ggCAAATTACCGGTAAACAGGGTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt             <  1:26434706/66‑1 (MQ=255)
                                                             aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCCCGATGGCACCAGGTAGTACCa         >  1:21971185/1‑65 (MQ=21)
                                                                aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct   <  1:23795700/68‑1 (MQ=255)
                                                                aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct   <  1:28473005/68‑1 (MQ=255)
                                                                aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct   <  1:13002562/68‑1 (MQ=255)
                                                                aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt      >  1:16700302/1‑65 (MQ=255)
                                                                aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt      >  1:16041067/1‑65 (MQ=255)
                                                                 ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg  <  1:23380454/68‑1 (MQ=255)
                                                                 ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg  <  1:22948948/68‑1 (MQ=255)
                                                                 |                                                                   
ACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTG  >  CP009273/2762239‑2762371

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