Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,676 C→T 29.7% H366H (CAC→CAT 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,6760CT29.7% 731.3 / inf 419H366H (CAC→CATicde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (188/106);  new base T (57/67);  total (245/174)
Fisher's exact test for biased strand distribution p-value = 7.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 457 total aligned reads displayed.

AATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCG  >  CP009273/1191609‑1191743
                                                                   |                                                                    
aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                      <  1:35198739/68‑1 (MQ=255)
aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                      >  1:7715135/1‑68 (MQ=255)
aaTATGCCGGACAGGACAAAGTAAATCCTGGCTCGATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                      >  1:2443059/1‑68 (MQ=255)
   atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATg                                                                   <  1:36926574/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATgg                                                                  <  1:19066389/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATgg                                                                  <  1:12464971/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                  <  1:38170541/68‑1 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                 >  1:29794945/1‑68 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                 <  1:15002780/68‑1 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                 >  1:9354366/1‑68 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                 >  1:14390918/1‑68 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                 >  1:5903764/1‑68 (MQ=255)
          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTgg                                                             <  1:13496801/67‑1 (MQ=255)
          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTgg                                                             <  1:32771841/67‑1 (MQ=255)
          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTgg                                                             >  1:29981498/1‑67 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGActg                                                         <  1:2049202/66‑3 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                       >  1:12336752/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                       >  1:21781125/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                       >  1:10660043/1‑68 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCg                                                    >  1:971038/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCg                                                    >  1:5728690/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCg                                                    >  1:924513/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                    >  1:26038580/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                    >  1:19642257/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                    >  1:22543273/1‑66 (MQ=255)
                      aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGc                                                  <  1:16018566/66‑1 (MQ=255)
                      aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTg                                               <  1:38409950/68‑1 (MQ=255)
                      aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTg                                               <  1:14668687/68‑1 (MQ=255)
                        aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc                                             <  1:3562713/68‑1 (MQ=255)
                        aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTg                                               >  1:36418394/1‑66 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                              >  1:30788902/1‑66 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                              >  1:1319288/1‑66 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            >  1:4850554/1‑68 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            <  1:30462452/68‑1 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            >  1:2383377/1‑68 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            <  1:5310973/68‑1 (MQ=255)
                          ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACtt                                           >  1:35873234/1‑68 (MQ=255)
                           cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                          >  1:19389669/1‑68 (MQ=255)
                               ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                      >  1:37487994/1‑68 (MQ=17)
                               ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                      >  1:17055369/1‑68 (MQ=255)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                      <  1:24846185/67‑1 (MQ=17)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                     <  1:18623120/68‑1 (MQ=17)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                     <  1:33151601/68‑1 (MQ=17)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                     <  1:38769704/68‑1 (MQ=255)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                     <  1:11889752/68‑1 (MQ=255)
                                 ctATTATTCTCTGCGCTGAGATGATGCTGCGCCAGATGGGTTGGACCGAAGCGGCTTACTTAATTGt                                     <  1:23926510/67‑1 (MQ=37)
                                    ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                 <  1:10952025/68‑1 (MQ=255)
                                     tattCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                <  1:4331929/68‑1 (MQ=14)
                                     tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                <  1:1121125/68‑1 (MQ=255)
                                     tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                <  1:33967192/68‑1 (MQ=255)
                                     tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACAGAAGCGGCTGACTTAATTGTTAAAg                                >  1:2763074/1‑68 (MQ=255)
                                      tttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:23010134/67‑1 (MQ=255)
                                      attCTGTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:33376504/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                               <  1:28704712/68‑1 (MQ=14)
                                      attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                               <  1:30273296/68‑1 (MQ=14)
                                      attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                               <  1:14811002/68‑1 (MQ=14)
                                      attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                               <  1:13553184/68‑1 (MQ=14)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:39174069/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:3944411/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:20115036/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:37736472/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:23460855/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:25554328/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:33656338/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:32340086/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:31366841/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:11441118/68‑1 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:2917634/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:13808309/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:34749006/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:718044/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:39546708/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:36772529/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                             >  1:5474546/1‑68 (MQ=11)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:8946777/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:10021107/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:14163073/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:25413835/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:27938861/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:27715879/1‑68 (MQ=255)
                                             ccGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGaa                        >  1:30440419/1‑68 (MQ=255)
                                             ccGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGaa                        >  1:37908238/1‑68 (MQ=255)
                                               gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                      <  1:8631323/68‑1 (MQ=255)
                                               gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                      <  1:8189306/68‑1 (MQ=255)
                                                cTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                     >  1:7394266/1‑68 (MQ=255)
                                                 tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                     >  1:6809873/1‑67 (MQ=255)
                                                    gatgatGCTGCGCCACATGGTTTGGACGAAACGGCGTGCCTTAATGTTAAAGGTATTGAAGGCCCaa                 >  1:1631780/1‑67 (MQ=9)
                                                      tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc               >  1:17990475/1‑68 (MQ=255)
                                                      tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc               >  1:27270001/1‑68 (MQ=255)
                                                         tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc            >  1:34432320/1‑68 (MQ=255)
                                                         tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc            >  1:3104627/1‑68 (MQ=255)
                                                          gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc            >  1:15723313/1‑67 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc          >  1:6289246/1‑68 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc          >  1:16755970/1‑68 (MQ=255)
                                                            tgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgcg         <  1:17685873/68‑1 (MQ=255)
                                                               gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa      >  1:21192740/1‑68 (MQ=255)
                                                                 cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa     <  1:10470622/67‑1 (MQ=255)
                                                                 cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc    <  1:22879549/68‑1 (MQ=255)
                                                                 cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc    <  1:28340200/68‑1 (MQ=255)
                                                                  acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc   <  1:26690466/68‑1 (MQ=255)
                                                                   |                                                                    
AATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCG  >  CP009273/1191609‑1191743

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: