Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,688 C→T 22.2% 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,6880CT22.2% 1027.3 / 288.9 443T370T (ACC→ACTicde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (200/144);  new base T (48/50);  total (249/194)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 457 total aligned reads displayed.

ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACT  >  CP009273/1191624‑1191755
                                                                |                                                                    
aCAAAGTAAATCCTGGCTCTATTATTCTTTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                   >  1:27890539/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                                   >  1:10660043/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                                   >  1:21781125/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                   >  1:27634177/1‑68 (MQ=255)
aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATAATGCTGCGCCATATGGGTTGGACTGaa                                                                   >  1:30522597/1‑68 (MQ=37)
   aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCg                                                                <  1:3045746/68‑1 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                                >  1:34369573/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                                >  1:26562594/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                                >  1:9247988/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                                >  1:9123477/1‑66 (MQ=255)
     gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                                >  1:22593318/1‑66 (MQ=255)
        aaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGa                                                          <  1:14912263/68‑1 (MQ=255)
         aTCCTGTATCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                                         <  1:17593443/68‑1 (MQ=255)
         aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                                         <  1:36800342/68‑1 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                                          >  1:1319288/1‑66 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                                          >  1:32923301/1‑66 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        >  1:37421567/1‑68 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        <  1:5310973/68‑1 (MQ=255)
          tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                        >  1:7746350/1‑68 (MQ=255)
          tCCTGGCTCAATTATTCTCTCCGCTGAGATGATGCAGCGCCATATGGTTTGGACTGAAGCGGCTGAcc                                                        >  1:4285871/1‑67 (MQ=21)
           ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                                          >  1:35465346/1‑65 (MQ=255)
           ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACtt                                                       >  1:1609980/1‑68 (MQ=255)
            cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctg                                                      >  1:7244260/1‑65 (MQ=21)
            cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                      >  1:30320740/1‑68 (MQ=255)
                ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                  >  1:37726194/1‑68 (MQ=255)
                ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                  >  1:17055369/1‑68 (MQ=255)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                  <  1:23005192/67‑1 (MQ=17)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                 <  1:24712715/68‑1 (MQ=17)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                 <  1:11889752/68‑1 (MQ=255)
                 tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCTGACTTAATTGt                                                 <  1:37990376/68‑1 (MQ=255)
                 tctATTATTCTCACCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGCAGCGGCTGACTTAATTGt                                                 <  1:8258136/68‑1 (MQ=255)
                     ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                             <  1:10952025/68‑1 (MQ=255)
                     ttattCTCGCCGCTGAGATGATGCTGCGCCACATGGGTTCGACCGAAGCGGATGAATTAATTGTTaaa                                             >  1:27156090/1‑68 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                           <  1:13553184/68‑1 (MQ=14)
                       attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                           <  1:2359402/68‑1 (MQ=14)
                       attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                           <  1:27367304/68‑1 (MQ=14)
                       attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                           <  1:29530465/68‑1 (MQ=14)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:14151463/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:2133593/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:39174069/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:13497185/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:3944411/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:34953933/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:31366841/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:16843348/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:5398667/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:11441118/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:5849983/68‑1 (MQ=255)
                       attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                           <  1:6590777/68‑1 (MQ=255)
                        ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                          <  1:29984540/68‑1 (MQ=255)
                         tctctcCGCTGAGGTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:38299685/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                         >  1:28915819/1‑68 (MQ=11)
                         tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                         <  1:36034261/68‑1 (MQ=11)
                         tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                         >  1:39546708/1‑68 (MQ=11)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:10021107/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:25413835/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:24167966/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:4865150/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:1471927/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                         >  1:9651402/1‑68 (MQ=255)
                         tctctcCGCTGAGATGATGCTGCGACACATGGGTTGGACCGAAGCGTCTGACTTAATTGTTAAAGGTa                                         >  1:11123273/1‑68 (MQ=25)
                           tctcCGCTGAGATGGTGCTGCGCCACATGGGTGGGACCGAAGCGGCTGGGTTAATTGTTAAAGGTATg                                       >  1:28559491/1‑68 (MQ=21)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:19372987/1‑68 (MQ=255)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:37814002/1‑68 (MQ=255)
                           tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                       >  1:15114413/1‑68 (MQ=255)
                                gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                  <  1:33552556/68‑1 (MQ=255)
                                gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                                   >  1:1901114/1‑67 (MQ=255)
                                     gatgatGCTGCGCCACATGGTTTGGACGAAACGGCGTGCCTTAATGTTAAAGGTATTGAAGGCCCaa                             >  1:1631780/1‑67 (MQ=9)
                                     gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                             >  1:17960990/1‑68 (MQ=255)
                                      ctaatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                             >  1:20203228/4‑67 (MQ=255)
                                          tgctgcGCCACATGGGTTGGCCCGCAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:30858457/1‑68 (MQ=255)
                                          tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:6211303/1‑68 (MQ=255)
                                          tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:38527340/1‑68 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:23900226/1‑67 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:15723313/1‑67 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:22743775/1‑67 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                        >  1:569288/1‑67 (MQ=255)
                                           gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTCTTAAAGGTATGGAACGCGCAATCAAc                        >  1:33199377/1‑67 (MQ=255)
                                            ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                      >  1:28191067/1‑68 (MQ=255)
                                            ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                      >  1:6975861/1‑68 (MQ=255)
                                            ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                      >  1:33816752/1‑68 (MQ=255)
                                             tgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgcg                     >  1:12392528/1‑68 (MQ=255)
                                                gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa                  >  1:14039259/1‑68 (MQ=255)
                                                  cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                 <  1:10470622/67‑1 (MQ=255)
                                                  cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                <  1:25938527/68‑1 (MQ=255)
                                                  cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                <  1:31209339/68‑1 (MQ=255)
                                                   acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc               <  1:7490749/68‑1 (MQ=255)
                                                       gggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTa            >  1:14416158/1‑67 (MQ=255)
                                                        ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa           >  1:18007505/1‑67 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:21026645/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:23280836/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:34605464/67‑1 (MQ=255)
                                                         gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc          <  1:3196648/67‑1 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:36505013/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:35758967/1‑68 (MQ=255)
                                                          ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        >  1:1270190/1‑68 (MQ=255)
                                                           tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt        <  1:12286779/67‑1 (MQ=255)
                                                              aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg     >  1:19737133/1‑67 (MQ=255)
                                                              aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg     >  1:38897849/1‑67 (MQ=255)
                                                              aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg     >  1:8748822/1‑67 (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: