Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,701 T→C 14.0% L375L (TTA→CTA) ‡ 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,7010TC14.0% 1265.1 / 139.8 406L375L (TTA→CTA) ‡icde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base T (177/171);  new base C (17/40);  total (194/212)
Fisher's exact test for biased strand distribution p-value = 3.98e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 430 total aligned reads displayed.

TCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGA  >  CP009273/1191634‑1191767
                                                                   |                                                                  
tacTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                                    >  1:35709249/3‑68 (MQ=255)
tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                                                    >  1:34934119/1‑68 (MQ=255)
 ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACtt                                                                   >  1:13193690/1‑68 (MQ=255)
  cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                                  >  1:30320740/1‑68 (MQ=255)
      ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                              >  1:37487994/1‑68 (MQ=17)
      ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                              >  1:17055369/1‑68 (MQ=255)
       tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                             <  1:18623120/68‑1 (MQ=17)
       tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                             <  1:11889752/68‑1 (MQ=255)
       tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCTGACTTAATTGt                                                             <  1:37990376/68‑1 (MQ=255)
        ctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGtt                                                            <  1:5924508/68‑1 (MQ=255)
           ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                         <  1:10952025/68‑1 (MQ=255)
            tattCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaagg                                                        <  1:39182011/68‑3 (MQ=14)
            tattCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                        <  1:21308340/68‑1 (MQ=14)
            tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                                        <  1:1121125/68‑1 (MQ=255)
             tttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:23010134/67‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                       <  1:13553184/68‑1 (MQ=14)
             attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                       <  1:28704712/68‑1 (MQ=14)
             attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                       <  1:14811002/68‑1 (MQ=14)
             attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                       <  1:5827515/68‑1 (MQ=14)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:37736472/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:4359899/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:3944411/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:16843348/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:20115036/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:23460855/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:36153324/68‑1 (MQ=255)
             attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                       <  1:33082999/68‑1 (MQ=255)
              ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                                      <  1:29984540/68‑1 (MQ=255)
               tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                     >  1:32415333/1‑68 (MQ=11)
               tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                     >  1:39546708/1‑68 (MQ=11)
               tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                     >  1:5474546/1‑68 (MQ=11)
               tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                     >  1:718044/1‑68 (MQ=11)
               tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                     >  1:13808309/1‑68 (MQ=11)
               tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                     >  1:14331811/1‑68 (MQ=11)
               tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                     >  1:38580483/1‑68 (MQ=255)
               tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                     >  1:10021107/1‑68 (MQ=255)
               tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                     >  1:28108444/1‑68 (MQ=255)
               tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                     >  1:38348217/1‑68 (MQ=255)
               tctctcAGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                     >  1:15309661/1‑68 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                   >  1:5183751/1‑68 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                   >  1:9131772/1‑68 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                   >  1:11585272/1‑68 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                   >  1:15961524/1‑68 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                   >  1:27715879/1‑68 (MQ=255)
                      gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                              <  1:22835059/68‑1 (MQ=255)
                      gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                              <  1:8189306/68‑1 (MQ=255)
                       cTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                                             >  1:7394266/1‑68 (MQ=255)
                           gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                         >  1:3355986/1‑68 (MQ=255)
                           gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                         >  1:3781492/1‑68 (MQ=255)
                            ctaatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                         >  1:3657684/4‑67 (MQ=255)
                             tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc                                       >  1:27270001/1‑68 (MQ=255)
                             tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc                                       >  1:24450500/1‑68 (MQ=255)
                                tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                    >  1:3104627/1‑68 (MQ=255)
                                 gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                    >  1:30452292/1‑67 (MQ=255)
                                 gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTCTTAAAGGTATGGAACGCGCAATCAAc                                    >  1:33199377/1‑67 (MQ=255)
                                  ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                                  >  1:6975861/1‑68 (MQ=255)
                                  ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                                  >  1:33816752/1‑68 (MQ=255)
                                        caccTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                            <  1:6726838/68‑1 (MQ=255)
                                        cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                             <  1:10470622/67‑1 (MQ=255)
                                        cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                            <  1:10585105/68‑1 (MQ=255)
                                        cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                            <  1:25938527/68‑1 (MQ=255)
                                        cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                            <  1:28340200/68‑1 (MQ=255)
                                         acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                           <  1:1927856/68‑1 (MQ=255)
                                         acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                           <  1:31492616/68‑1 (MQ=255)
                                         acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                           <  1:36887771/68‑1 (MQ=255)
                                              ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                       >  1:5340584/1‑67 (MQ=255)
                                              ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                       >  1:18007505/1‑67 (MQ=255)
                                               gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                      <  1:34605464/67‑1 (MQ=255)
                                               gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                      <  1:12312232/67‑1 (MQ=255)
                                               gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                      <  1:15028745/67‑1 (MQ=255)
                                               gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                      <  1:21026645/67‑1 (MQ=255)
                                               gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                      <  1:31245908/67‑1 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCCACGCGAAAACCGTAACCt                    >  1:30762384/1‑68 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                    >  1:7623739/1‑68 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                    >  1:1270190/1‑68 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                    >  1:23829843/1‑68 (MQ=255)
                                                 tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                    <  1:14147757/67‑1 (MQ=255)
                                                 tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                    <  1:22644340/67‑1 (MQ=255)
                                                  ggACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                  <  1:1241007/68‑1 (MQ=255)
                                                    aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                 >  1:6329651/1‑67 (MQ=255)
                                                    aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                 >  1:38897849/1‑67 (MQ=255)
                                                    aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                 >  1:19737133/1‑67 (MQ=255)
                                                     ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                 >  1:17648366/1‑66 (MQ=255)
                                                     ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                 >  1:34036057/1‑66 (MQ=255)
                                                     ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGa                <  1:17912655/67‑1 (MQ=255)
                                                          gCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGa          >  1:29089719/1‑68 (MQ=255)
                                                            ggCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg         <  1:26427940/67‑1 (MQ=255)
                                                              cTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCg       <  1:16244537/67‑1 (MQ=255)
                                                              cTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCg       <  1:8936783/67‑1 (MQ=255)
                                                               tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGt      <  1:25256239/67‑1 (MQ=255)
                                                               tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGt      <  1:29811269/67‑1 (MQ=255)
                                                               tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGt      <  1:9331259/67‑1 (MQ=255)
                                                               tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc     <  1:24133080/68‑1 (MQ=255)
                                                               tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc     <  1:32639558/68‑1 (MQ=255)
                                                                gACTTAATTGTTAAAGGTATTGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt    <  1:8611353/68‑1 (MQ=255)
                                                                gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt    <  1:14557697/68‑1 (MQ=255)
                                                                gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt    <  1:13504481/68‑1 (MQ=255)
                                                                gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt    <  1:26974404/68‑1 (MQ=255)
                                                                  cTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga  <  1:22094818/68‑1 (MQ=255)
                                                                  cTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga  <  1:16593066/68‑1 (MQ=255)
                                                                   |                                                                  
TCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGA  >  CP009273/1191634‑1191767

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: