Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,688 C→T 5.3% T370T (ACC→ACT icd → e14 prophage; isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,191,6880CT 5.3% 1547.9 / 51.3 456T370T (ACC→ACTicde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (244/185);  new base T (9/15);  total (255/201)
Fisher's exact test for biased strand distribution p-value = 8.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 471 total aligned reads displayed.

ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACT  >  CP009273/1191624‑1191755
                                                                |                                                                   
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                                  >  1:20894054/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCg                                                                    <  1:1746082/66‑1 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                  >  1:29805437/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                  >  1:25227110/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                  >  1:7607622/1‑68 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                               >  1:10299906/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                               >  1:18521032/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                               >  1:16972381/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                               >  1:24617304/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                               >  1:13681774/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                               >  1:26897790/1‑66 (MQ=255)
         aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                                         <  1:4128324/68‑1 (MQ=255)
         aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                                         <  1:31746515/68‑1 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        >  1:2644430/1‑68 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        >  1:33957271/1‑68 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        >  1:8144084/1‑68 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        >  1:26180782/1‑68 (MQ=255)
          tCCTGGCTCTATAATTCTCTCCGCTGAGATGATGCTGCGCCACATGCGTTGGACCGAAGCGGCTGACt                                                        >  1:17879533/1‑68 (MQ=255)
            cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctg                                                      >  1:645500/1‑65 (MQ=21)
            cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                      >  1:20735903/1‑68 (MQ=255)
                ctctATTATTATCTCCGCGGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCGGACGTAATTg                                                  >  1:3846878/1‑68 (MQ=255)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                  <  1:15131119/67‑1 (MQ=255)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                 <  1:22621628/68‑1 (MQ=255)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                 <  1:2005992/68‑1 (MQ=255)
                 tctATTATGCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                  <  1:32952793/67‑1 (MQ=255)
                      tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                            >  1:4680002/1‑68 (MQ=255)
                      tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                            >  1:32615077/1‑68 (MQ=255)
                      tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                            <  1:3625921/68‑1 (MQ=255)
                      tattCTCTCCGCTGAGATGATGCTGCGCCACATGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:28012101/68‑1 (MQ=255)
                       attCTCTCCGCTGTGATGATGCTGCGCCACTTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:27256027/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGTTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:15266343/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                           <  1:21266942/68‑1 (MQ=14)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           >  1:32302476/1‑68 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:8989162/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:23817710/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:7816895/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGATACTTTGGTTGGACCGAAGAGGCTTACTTAATTGTTATAgg                                           <  1:26684265/68‑1 (MQ=255)
                        ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                          <  1:28888567/68‑1 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                         >  1:19791000/1‑68 (MQ=11)
                         tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                         >  1:4346311/1‑68 (MQ=11)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         <  1:10015693/68‑1 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:28481017/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:2311947/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:571433/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:12494373/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:30695674/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:7133584/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACATAATTGTTAAAGGTa                                         >  1:8526600/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGCCTGACTTAATTGTTAAAGGTa                                         >  1:14517723/1‑68 (MQ=255)
                          ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                                        <  1:5911804/68‑1 (MQ=255)
                           tctcCGCTGAGATGATGCTTCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:15863665/1‑68 (MQ=255)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:22153790/1‑68 (MQ=255)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:29686656/1‑68 (MQ=255)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:13047277/1‑68 (MQ=255)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:24298256/1‑68 (MQ=255)
                                gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                  <  1:33050204/68‑1 (MQ=255)
                                     gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                             >  1:11618673/1‑68 (MQ=255)
                                      atgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                             >  1:14180861/1‑67 (MQ=255)
                                       tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc                           <  1:11303568/68‑1 (MQ=255)
                                          tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:3159500/1‑68 (MQ=255)
                                          tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGAAATCAAc                        >  1:20485067/1‑68 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                       <  1:94100/68‑1 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                       <  1:23997076/68‑1 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:6775188/1‑67 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:34362465/1‑67 (MQ=255)
                                                  cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                <  1:23390327/68‑1 (MQ=255)
                                                   acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc               <  1:10589887/68‑1 (MQ=255)
                                                     aTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGt             >  1:14652357/1‑68 (MQ=255)
                                                     aTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGt             >  1:30946809/1‑68 (MQ=255)
                                                     aTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGt             >  1:21953936/1‑68 (MQ=255)
                                                     aTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGt             >  1:2593547/1‑68 (MQ=255)
                                                        ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa           >  1:16459852/1‑67 (MQ=255)
                                                        ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa           >  1:17573723/1‑67 (MQ=255)
                                                        ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa           >  1:18248580/1‑67 (MQ=255)
                                                        ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:27671547/68‑1 (MQ=255)
                                                        ggTAGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa           >  1:16151264/1‑67 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:30532534/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:32757457/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:22879326/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:12046548/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:21598777/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:19374059/67‑1 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:15535686/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:33800231/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:5191918/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:25354800/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:25580334/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:12763654/1‑68 (MQ=255)
                                                           tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        <  1:13309655/67‑1 (MQ=255)
                                                           tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        <  1:33688859/67‑1 (MQ=255)
                                                           tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        <  1:33363265/67‑1 (MQ=255)
                                                           tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTa       >  1:12202589/1‑68 (MQ=255)
                                                           tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTa       >  1:31297577/1‑68 (MQ=255)
                                                            ggACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt      <  1:113735/68‑1 (MQ=255)
                                                            ggACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt      <  1:10841498/68‑1 (MQ=255)
                                                              aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg     >  1:29495169/1‑67 (MQ=255)
                                                              aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg     >  1:1891447/1‑67 (MQ=255)
                                                               ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg     >  1:33110743/1‑66 (MQ=255)
                                                               ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGa    <  1:27852966/67‑1 (MQ=255)
                                                                cGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACt  <  1:5458296/68‑1 (MQ=255)
                                                                |                                                                   
ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACT  >  CP009273/1191624‑1191755

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: