Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 930,950 A→C 7.8% Q757H (CAA→CAC ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273930,9500AC 7.8% 1007.8 / 5.8 415Q757H (CAA→CACftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base A (141/206);  new base C (31/4);  total (203/212)
Fisher's exact test for biased strand distribution p-value = 3.27e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 434 total aligned reads displayed.

TTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATC  >  CP009273/930883‑931014
                                                                   |                                                                
ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGcaacaa                                                                  <  1:25990797/68‑1 (MQ=255)
 tgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGcaacaa                                                                  >  1:2481493/1‑67 (MQ=255)
       aCCTGTACAGCAGCCGCAACAACCGGTTGCCCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGtt                                                            >  1:38113076/1‑67 (MQ=255)
       aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCCGCCGGCACATCCAGTTc                                                           >  1:37491828/1‑68 (MQ=18)
       aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGtt                                                            >  1:29197652/1‑67 (MQ=255)
       aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGtt                                                            >  1:27849716/1‑67 (MQ=255)
       aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGtt                                                            >  1:27242328/1‑67 (MQ=255)
        ccTGTACAGCAGCCGCAACAACCGGTTGCCCCGCAGCAGCAATATCAGCAGCCGCAACCACCAGTTcc                                                          >  1:36347512/1‑68 (MQ=17)
        ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTcc                                                          >  1:22904102/1‑68 (MQ=255)
        ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTcc                                                          >  1:38468818/1‑68 (MQ=255)
        ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTcc                                                          >  1:20063455/1‑68 (MQ=255)
        ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTcc                                                          >  1:1008659/1‑68 (MQ=255)
         cTGTACAGCAGCCGCAACACCCGGTTGCACCGCAGCAGCATTATCAGCAGCCTCAACAGCCAGTTcc                                                          >  1:22499171/1‑67 (MQ=14)
         cTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTccg                                                         <  1:14321240/68‑1 (MQ=255)
         cTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTcc                                                          >  1:1800525/1‑67 (MQ=255)
            tACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTccgcc                                                       >  1:23265582/1‑67 (MQ=255)
            tACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTccgcc                                                       >  1:30420186/1‑67 (MQ=255)
             aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTGcgccgc                                                     >  1:13198907/1‑68 (MQ=255)
             aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCcgccgc                                                     >  1:26139959/1‑68 (MQ=255)
             aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCcgccgc                                                     >  1:34834687/1‑68 (MQ=255)
             aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCcgccgc                                                     >  1:3578146/1‑68 (MQ=255)
             aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCcgccgc                                                     >  1:21252469/1‑68 (MQ=255)
               agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAGCCAGTTCCgccgca                                                    >  1:20485127/1‑67 (MQ=255)
               agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTCCgccgca                                                    >  1:15795505/1‑67 (MQ=255)
               agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgca                                                    >  1:1749527/1‑67 (MQ=255)
                gcagcCGCAACAACCGGTTGCCCCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTCCgccgcag                                                   >  1:29356471/1‑67 (MQ=14)
                gcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAGCCAGTTGCgccgcag                                                   >  1:27038458/1‑67 (MQ=255)
                  agcCGCTCCAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACTACCAGTTCCgccgcagca                                                 <  1:22027001/67‑2 (MQ=17)
                  agcCGCAACAACCGGTTGCCCCGCAGCAGCAATATCAGCAGCCGCAACATCCAGTTCCgccgcagcc                                                 >  1:11583278/1‑67 (MQ=17)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCCGCCGCAACAGCCAGTTCCgccgcagcc                                                 >  1:6345097/1‑67 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACATCCAGTTCCgccgcagcc                                                 >  1:34712505/1‑67 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagcc                                                 >  1:33964758/1‑67 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagcc                                                 <  1:13491760/67‑1 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagca                                                 >  1:33571820/1‑66 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagc                                                  <  1:455115/66‑1 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagc                                                  <  1:11275134/66‑1 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcag                                                   <  1:18396018/65‑1 (MQ=255)
                  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgca                                                    <  1:35674858/64‑1 (MQ=255)
                     cGCAACAACCGGTTTCACCTGAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:36954406/68‑1 (MQ=17)
                     cGCAACAACCGGTTGCGCCGCAGCAGCAATATCAGCAGCCGCAACAGCCAGTTGCGccgcagccgcag                                             >  1:25494019/1‑68 (MQ=11)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagccgca                                              <  1:6020622/67‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagccgca                                              <  1:16809404/67‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagccgca                                              <  1:31216597/67‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcagccg                                                <  1:24732070/65‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:3028738/68‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:10339533/68‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:14004381/68‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:27994594/68‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:36481642/68‑1 (MQ=255)
                     cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                             <  1:12759321/68‑1 (MQ=255)
                      gCAACAACCGGTTGCACCGCAGCAGCAATATCTGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGt                                            <  1:10992224/68‑1 (MQ=255)
                      gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGt                                            <  1:21486847/68‑1 (MQ=255)
                      gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGt                                            <  1:20293264/68‑1 (MQ=255)
                      gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGt                                            <  1:791805/68‑1 (MQ=255)
                        aacaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTa                                           <  1:23843424/67‑1 (MQ=255)
                         acaacCGGTTGCACCGCAGCAGCAATATCACCAGCCGCAACACCCAGTTCCGCCCCACCCGCAGTAta                                         >  1:39059379/1‑67 (MQ=17)
                          caacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATca                                        >  1:29778383/1‑68 (MQ=255)
                           aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTGCGCCGCAGCCGCAGTATca                                        >  1:27538781/1‑67 (MQ=255)
                              cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagcag                                    <  1:14651377/68‑1 (MQ=255)
                              cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagcag                                    <  1:3330925/68‑1 (MQ=255)
                              cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagcag                                    <  1:38865479/68‑1 (MQ=255)
                              cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagcag                                    <  1:31122745/68‑1 (MQ=255)
                              cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagcag                                    <  1:10647946/68‑1 (MQ=255)
                              cGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcagca                                     >  1:7530829/1‑67 (MQ=255)
                                gTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGcc                                  >  1:32192086/1‑68 (MQ=255)
                                   gCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCcaca                               <  1:2417060/68‑1 (MQ=255)
                                   gCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCcaca                               >  1:14961222/1‑68 (MQ=255)
                                     aCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCACCAGCCTCAAc                             <  1:9199925‑M2/68‑2 (MQ=255)
                                      ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTGCGCCGCAGCCGCAGTATCAGCAGCCTCAAcc                            <  1:29617292‑M2/68‑3 (MQ=255)
                                           gcagcaATATCAGCAGCCGCAACAGCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCCCAAcagccgg                       <  1:33436356‑M2/68‑8 (MQ=255)
                                             agcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAAcagccggtt                     <  1:19843345‑M2/68‑10 (MQ=255)
                                                aaTATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                   >  1:19591712/1‑67 (MQ=255)
                                                aaTATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                   >  1:21567538/1‑67 (MQ=255)
                                                    tCAGCAGCCGCAACACCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCAGGTTGcgccgc              >  1:6900058/1‑68 (MQ=255)
                                                     cggcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCCGTTGcgccgc              >  1:5750474/3‑67 (MQ=255)
                                                     cagcagCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccac              >  1:5324244/1‑65 (MQ=255)
                                                      agcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCa             >  1:33833450/1‑67 (MQ=255)
                                                      agcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCa             >  1:35266246/1‑67 (MQ=255)
                                                         agcCGGACCACCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc          <  1:7214349/67‑1 (MQ=18)
                                                         agcCGCAACAACCAGTTCCGCCGCAGGCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGc           <  1:12134921/66‑1 (MQ=255)
                                                         agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc          >  1:9758382/1‑67 (MQ=255)
                                                         agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc          >  1:34385934/1‑67 (MQ=255)
                                                         agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc          >  1:28585737/1‑67 (MQ=255)
                                                         agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGc           <  1:19275943/66‑1 (MQ=255)
                                                         agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGc           <  1:25674919/66‑1 (MQ=255)
                                                         agcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAAAAGCCGGTTGCGCCGCAGc           <  1:28317653/66‑1 (MQ=255)
                                                          gcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCCGCAGCCACAACAGCCGGTTGCGCCGCAGcccc        >  1:421030/1‑66 (MQ=255)
                                                            ctcaaccaCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCcaca       <  1:16187638/65‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCGGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa      <  1:8349708/68‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCcaca       <  1:32318515/67‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCcaca       <  1:28876808/67‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCca         <  1:34028196/65‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCca         <  1:32620645/65‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa      <  1:5038548/68‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa      <  1:17233557/68‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa      <  1:31599946/68‑1 (MQ=255)
                                                            cGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa      <  1:1594835/68‑1 (MQ=255)
                                                            cGCAACAACCAGTCCCGCCGCAGCTGCAGTATCAGCACCCACAACAGCCGGTTGCGCCGCAGCcaca       <  1:8672213/67‑1 (MQ=18)
                                                             gCAACAACGAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAat     <  1:18865131/68‑1 (MQ=255)
                                                               aacaacCATTGCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAatat   <  1:15631129/68‑1 (MQ=255)
                                                                   |                                                                
TTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATC  >  CP009273/930883‑931014

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: