Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,703 A→G 14.3% L375L (TTA→TTG) ‡ 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,7030AG14.3% 1251.0 / 148.0 405L375L (TTA→TTG) ‡icde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base A (176/171);  new base G (17/41);  total (193/212)
Fisher's exact test for biased strand distribution p-value = 2.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

Only 100 of 408 total aligned reads displayed.

CTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATG  >  CP009273/1191636‑1191769
                                                                    |                                                                  
    ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                                >  1:30403762/1‑68 (MQ=17)
    ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                                >  1:37726194/1‑68 (MQ=255)
     tctATTATTCTCTCCGCTTAGATGATGCTGCGCGACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                               <  1:14399095/68‑1 (MQ=255)
     tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                                <  1:5502603/67‑1 (MQ=17)
     tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                                <  1:23005192/67‑1 (MQ=17)
     tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                               <  1:24712715/68‑1 (MQ=17)
     tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTg                                                                <  1:31072420/67‑1 (MQ=255)
     tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGt                                                               <  1:30498566/68‑1 (MQ=255)
     tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCTGACTTAATTGt                                                               <  1:37990376/68‑1 (MQ=255)
      ctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGtt                                                              <  1:37095369/68‑1 (MQ=17)
         ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                           <  1:10952025/68‑1 (MQ=255)
          tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACAGAAGCGGCTGACTTAATTGTTAAAg                                                          >  1:2763074/1‑68 (MQ=255)
           attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                         >  1:25975931/1‑68 (MQ=14)
           attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                         <  1:5827515/68‑1 (MQ=14)
           attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                         <  1:14811002/68‑1 (MQ=14)
           attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                         <  1:13553184/68‑1 (MQ=14)
           attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                         <  1:12187405/68‑1 (MQ=14)
           attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                         <  1:21149081/68‑1 (MQ=14)
           attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                         <  1:39308860/68‑1 (MQ=255)
           attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                         <  1:4359899/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:11441118/68‑1 (MQ=255)
           attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                         <  1:32340086/68‑1 (MQ=255)
           attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                         <  1:33082999/68‑1 (MQ=255)
            ttCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACATGATTGTTAAAgg                                                         <  1:17728617/67‑1 (MQ=14)
             tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                       >  1:39546708/1‑68 (MQ=11)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:25379979/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:38580483/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:38348217/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:24167966/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:31751165/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:10021107/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:14624980/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                       >  1:12607091/1‑68 (MQ=255)
             tctctcCGCTGAGATGATGCTGCGACACATGGGTTGGACCGAAGCGTCTGACTTAATTGTTAAAGGTa                                                       >  1:11123273/1‑68 (MQ=25)
               tctcCGCTGAGATGGTGCTGCGCCACATGGGTGGGACCGAAGCGGCTGGGTTAATTGTTAAAGGTATg                                                     >  1:28559491/1‑68 (MQ=21)
               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                     >  1:5183751/1‑68 (MQ=255)
               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                     >  1:19372987/1‑68 (MQ=255)
               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                     >  1:9131772/1‑68 (MQ=255)
               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                     >  1:15961524/1‑68 (MQ=255)
               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATg                                                     >  1:6999374/1‑68 (MQ=255)
               tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAGTGTTATAGGTATg                                                     >  1:15041939/1‑68 (MQ=255)
                    gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                                <  1:22835059/68‑1 (MQ=255)
                    gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                                <  1:8631323/68‑1 (MQ=255)
                    gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                                                 >  1:1901114/1‑67 (MQ=255)
                    gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                                                 >  1:31366794/1‑67 (MQ=255)
                      tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                                               >  1:6809873/1‑67 (MQ=255)
                         gatgatGCTGCGCCACATGGTTTGGACGAAACGGCGTGCCTTAATGTTAAAGGTATTGAAGGCCCaa                                           >  1:1631780/1‑67 (MQ=9)
                         gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                           >  1:3781492/1‑68 (MQ=255)
                         gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                           >  1:3355986/1‑68 (MQ=255)
                         gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                           >  1:29913868/1‑68 (MQ=255)
                         gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                           >  1:8226801/1‑68 (MQ=255)
                          ctaatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                           >  1:20203228/4‑67 (MQ=255)
                              tgctgcGCCACATGGGTTGGCCCGCAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                      >  1:30858457/1‑68 (MQ=255)
                              tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                      >  1:17970107/1‑68 (MQ=255)
                               gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                                     <  1:11873176/68‑1 (MQ=255)
                               gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                      >  1:23900226/1‑67 (MQ=255)
                               gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTCTTAAAGGTATGGAACGCGCAATCAAc                                      >  1:33199377/1‑67 (MQ=255)
                                ctgcGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                                    >  1:28238130/1‑68 (MQ=255)
                                ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                                    >  1:32594834/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)
                                    gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa                                >  1:5372768/1‑68 (MQ=255)
                                     ccACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa                                >  1:34602740/1‑67 (MQ=255)
                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                               <  1:10470622/67‑1 (MQ=255)
                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                               <  1:23464040/67‑1 (MQ=255)
                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                              <  1:10585105/68‑1 (MQ=255)
                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc                              <  1:16910588/68‑1 (MQ=255)
                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGaaaa                               <  1:34755527/67‑1 (MQ=255)
                                       acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc                             <  1:29284940/68‑1 (MQ=255)
                                             gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                        <  1:36739740/67‑1 (MQ=255)
                                             gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                        <  1:12980751/67‑1 (MQ=255)
                                              ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                      >  1:35758967/1‑68 (MQ=255)
                                              ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                      >  1:7623739/1‑68 (MQ=255)
                                               tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                      <  1:19882362/67‑1 (MQ=255)
                                               tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                      <  1:14147757/67‑1 (MQ=255)
                                               tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                      <  1:22644340/67‑1 (MQ=255)
                                               tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                      <  1:27216660/67‑1 (MQ=255)
                                                 gACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                   >  1:15349689/1‑68 (MQ=255)
                                                  aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                   >  1:38897849/1‑67 (MQ=255)
                                                  aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                   >  1:17134628/1‑67 (MQ=255)
                                                   ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                   >  1:34036057/1‑66 (MQ=255)
                                                   ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                   >  1:7468210/1‑66 (MQ=255)
                                                       aGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCg             >  1:18261988/1‑68 (MQ=255)
                                                        gCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGa            >  1:26335293/1‑68 (MQ=255)
                                                            cTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCg         <  1:8936783/67‑1 (MQ=255)
                                                            cTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCg         <  1:6403401/67‑1 (MQ=255)
                                                              tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGt        <  1:9331259/67‑1 (MQ=255)
                                                              tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc       <  1:27910282/68‑1 (MQ=255)
                                                              tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc       <  1:30018496/68‑1 (MQ=255)
                                                              tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAACTATGACTTCGAGCGTc       <  1:2897652/68‑1 (MQ=255)
                                                               gACTTAATTGTTAAAGGTATGGAAGGCGCAATCATCGCGAAAACCGTAACCTATGACTTCGAGCGTCt      <  1:36177709/68‑1 (MQ=255)
                                                               gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt      <  1:16738092/68‑1 (MQ=255)
                                                               gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt      <  1:7197586/68‑1 (MQ=255)
                                                               gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt      <  1:6271227/68‑1 (MQ=255)
                                                               gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt      <  1:26929008/68‑1 (MQ=255)
                                                                 cTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga    <  1:22094818/68‑1 (MQ=255)
                                                                    |                                                                  
CTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATG  >  CP009273/1191636‑1191769

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: