Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,676 C→T 7.9% 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,6760CT 7.9% 1340.4 / 77.0 417H366H (CAC→CATicde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (220/162);  new base T (13/20);  total (234/183)
Fisher's exact test for biased strand distribution p-value = 4.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 459 total aligned reads displayed.

AATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCG  >  CP009273/1191609‑1191743
                                                                   |                                                                   
aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                     <  1:6746874/68‑1 (MQ=255)
aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                     >  1:24498995/1‑68 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATgg                                                                 <  1:21773493/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:20436409/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:22858351/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:18563491/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:26940456/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCACCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:24663897/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCGCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:9548241/68‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGACTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:25992972/68‑1 (MQ=255)
    tGCCGGGCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                 <  1:28102550/68‑1 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                >  1:4710808/1‑68 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                >  1:13236713/1‑68 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                >  1:9256980/1‑68 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                <  1:823793/68‑1 (MQ=255)
     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTATCTCCGCTGAGATGATGCTGCGCCACATggg                                                                >  1:33962347/1‑68 (MQ=255)
          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTgg                                                            <  1:27864569/67‑1 (MQ=255)
          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTgg                                                            >  1:4473785/1‑67 (MQ=255)
           cAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGa                                                           <  1:14432119/67‑1 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                      >  1:4118461/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                      >  1:18196830/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                      >  1:17473973/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                      >  1:25227110/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                      >  1:11319004/1‑68 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                      >  1:33797685/1‑68 (MQ=255)
                  aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                   <  1:16186003/68‑1 (MQ=255)
                   aGTAAATCCTGGCTGTATTATGCTCTCCGCTGAGATGATGCTGCGCGACATGGGTTGGACCGAAGCgg                                                  <  1:29795638/68‑1 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGTGTTGGACCGAAGCg                                                   >  1:13644388/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                   >  1:30604907/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                   >  1:20939785/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                   >  1:10299906/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                   >  1:10905427/1‑66 (MQ=255)
                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                   >  1:16972381/1‑66 (MQ=255)
                        aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                             <  1:5329204/68‑1 (MQ=255)
                        aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                             <  1:13032429/68‑1 (MQ=255)
                        aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                             <  1:27586934/68‑1 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            >  1:5614039/1‑68 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            >  1:33652178/1‑68 (MQ=255)
                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                            >  1:32402164/1‑68 (MQ=255)
                          ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACtt                                           >  1:25307760/1‑68 (MQ=255)
                               ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                      >  1:20569103/1‑68 (MQ=255)
                               ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                      >  1:3846866/1‑68 (MQ=255)
                               ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                      >  1:29363886/1‑68 (MQ=255)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                      <  1:25408822/67‑1 (MQ=17)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                     <  1:30883177/68‑1 (MQ=17)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                     <  1:31666927/68‑1 (MQ=255)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                     <  1:2005992/68‑1 (MQ=255)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                     <  1:32939563/68‑1 (MQ=255)
                                tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                     >  1:12138104/1‑68 (MQ=255)
                                      attCTCTCCGCTGTGATGATGCTGCGCCACTTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:27256027/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGTTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:15266343/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTGAAgg                               <  1:22561381/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:11433081/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               >  1:26220482/1‑68 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:6985655/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:433046/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:12630381/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:7343784/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:19038237/68‑1 (MQ=255)
                                      attCTCTCCGCTGAGATGATGCTGCGATACTTTGGTTGGACCGAAGAGGCTTACTTAATTGTTATAgg                               <  1:26684265/68‑1 (MQ=255)
                                      aatctgtcCGCGGAGATGATGCTGCCCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                               <  1:29677926/66‑1 (MQ=21)
                                       ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                              <  1:28888567/68‑1 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:6459927/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:5969669/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:33029885/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:7704826/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             <  1:10015693/68‑1 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:3106245/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:23650731/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:15487359/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:14776920/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:28481017/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                             >  1:14544608/1‑68 (MQ=255)
                                        tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACATAATTGTTAAAGGTa                             >  1:8526600/1‑68 (MQ=255)
                                         ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                            >  1:34344135/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTTCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:15863665/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:7886629/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:1669112/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:22153790/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:23192782/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:25852656/1‑68 (MQ=255)
                                          tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                           >  1:13775508/1‑68 (MQ=255)
                                               gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                      <  1:3614021/68‑1 (MQ=255)
                                                 tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                     >  1:12356747/1‑67 (MQ=255)
                                                    gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                 >  1:30286851/1‑68 (MQ=255)
                                                    gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCa                  <  1:11654162/67‑1 (MQ=255)
                                                          gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg           <  1:23997076/68‑1 (MQ=255)
                                                          gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc            >  1:1787136/1‑67 (MQ=255)
                                                          gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc            >  1:5870343/1‑67 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc          >  1:26451009/1‑68 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc          >  1:21970295/1‑68 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc          >  1:886859/1‑68 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc          >  1:3223415/1‑68 (MQ=255)
                                                            tgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgcg         >  1:33067118/1‑68 (MQ=255)
                                                               gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa      >  1:19478624/1‑68 (MQ=255)
                                                                 cacaTGGGTTGGACCGATGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa     <  1:33217968/67‑1 (MQ=255)
                                                                 cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa     <  1:32645653/67‑1 (MQ=255)
                                                                 cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc    <  1:21404824/68‑1 (MQ=255)
                                                                 cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc    <  1:19853378/68‑1 (MQ=255)
                                                                  acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc   <  1:10589887/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: