Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,440,865 C→T 6.1% L247L (CTG→CTA) ‡ fabB ← 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,440,8650CT 6.1% 212.8 / 4.4 65L247L (CTG→CTA) ‡fabB3‑oxoacyl‑[acyl‑carrier‑protein] synthase I
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

TCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACA  >  NC_000913/2440720‑2440986
                                                                                                                                                 |                                                                                                                         
tCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAg                                                                                                                         <  2:410112/148‑1 (MQ=255)
    gTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGcgcg                                                                                                                     >  2:765627/1‑148 (MQ=255)
    gTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGcgcg                                                                                                                     <  1:765627/148‑1 (MQ=255)
                       aCGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCTGGCGCgtg                                                                                                                   <  3:1101588/131‑1 (MQ=255)
                       aCGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCTGGCGCgtg                                                                                                                   >  4:1101588/1‑131 (MQ=255)
                                    ccATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgcc                                                                                     >  1:1115337/1‑148 (MQ=255)
                                    ccATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgcc                                                                                     >  1:328393/1‑148 (MQ=255)
                                                 gCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTa                                                                                               >  4:867588/1‑125 (MQ=255)
                                                 gCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTa                                                                                               <  3:867588/125‑1 (MQ=255)
                                                                     aCGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCaa                                                                                                     <  4:1065257/99‑1 (MQ=255)
                                                                     aCGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCaa                                                                                                     >  3:1065257/1‑99 (MQ=255)
                                                                     aCGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGc                                                    >  2:764410/1‑148 (MQ=255)
                                                                        gAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAg                                                                                                             <  1:96225/88‑1 (MQ=255)
                                                                        gAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAg                                                                                                             >  2:96225/1‑88 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGcgcg                                                                                                                     <  4:148837/73‑1 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGcgcg                                                                                                                     >  3:148837/1‑73 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGctct                                                                                                         >  1:789961/1‑85 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGctct                                                                                                         <  2:789961/85‑1 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCa                                                            <  3:661578/130‑1 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCCCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCa                                                            >  4:661578/1‑130 (MQ=255)
                                                                                 tGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        <  2:328393/148‑1 (MQ=255)
                                                                                   tCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACggg                                      <  2:1115337/148‑1 (MQ=255)
                                                                                   tCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACggg                                      <  1:764410/148‑1 (MQ=255)
                                                                                       cACCATCAGAGGTTGCGCCGTAGCCAACGACTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACAACCATACCGCCGCCGCCAGCCATAACGAAACCGTCACGGTCAGCGTCGTAAGTACGGGAGGc                                  >  3:703330/1‑148 (MQ=255)
                                                                                         ccATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTa                                                                                               <  4:644770/85‑1 (MQ=255)
                                                                                         ccATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTa                                                                                               >  3:644770/1‑85 (MQ=255)
                                                                                         ccATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAAc                                                                 <  1:106098/115‑1 (MQ=255)
                                                                                         ccATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAAc                                                                 >  2:106098/1‑115 (MQ=255)
                                                                                             cAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCaa                                                                                                     >  1:841500/1‑75 (MQ=255)
                                                                                             cAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCaa                                                                                                     <  2:841500/75‑1 (MQ=255)
                                                                                             cAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATaa                                                                       <  2:579982/105‑1 (MQ=255)
                                                                                             cAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATaa                                                                       >  1:579982/1‑105 (MQ=255)
                                                                                                agtttGCGCCGTAGTCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGAGTTCGAGCTCTTCAACCACTACCATACCGCCGCTGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCgg                         <  4:703330/145‑1 (MQ=255)
                                                                                                agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGaa                                                                   <  3:168442/106‑1 (MQ=255)
                                                                                                agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGaa                                                                   >  4:168442/1‑106 (MQ=255)
                                                                                                agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        >  1:896571/1‑133 (MQ=255)
                                                                                                agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        <  2:896571/133‑1 (MQ=255)
                                                                                                                aaCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACggg                                      <  4:502276/119‑1 (MQ=255)
                                                                                                                aaCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACggg                                      >  3:502276/1‑119 (MQ=255)
                                                                                                                 aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        >  2:298510/1‑116 (MQ=255)
                                                                                                                 aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        <  1:1053281/116‑1 (MQ=255)
                                                                                                                 aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        >  2:1053281/1‑116 (MQ=255)
                                                                                                                 aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        <  1:298510/116‑1 (MQ=255)
                                                                                                                       tCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCgccgccgcc                                                                               <  3:191252/71‑1 (MQ=255)
                                                                                                                       tCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCgccgccgcc                                                                               >  4:191252/1‑71 (MQ=255)
                                                                                                                            aTAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        <  2:164136/105‑1 (MQ=255)
                                                                                                                            aTAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                        >  1:164136/1‑105 (MQ=255)
                                                                                                                             tAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAgg                                   >  4:959277/1‑109 (MQ=255)
                                                                                                                             tAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAgg                                   <  3:959277/109‑1 (MQ=255)
                                                                                                                               gATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACa  <  3:613017/140‑1 (MQ=255)
                                                                                                                               gATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACa  >  4:613017/1‑140 (MQ=255)
                                                                                                                                 tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGAc   >  3:803247/1‑137 (MQ=255)
                                                                                                                                 tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGAc   <  4:803247/137‑1 (MQ=255)
                                                                                                                                    gAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAg                                                     >  4:285742/1‑84 (MQ=255)
                                                                                                                                    gAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAg                                                     <  3:285742/84‑1 (MQ=255)
                                                                                                                                       caccacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGt                 >  4:449443/1‑117 (MQ=255)
                                                                                                                                       caccacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGt                 <  3:449443/117‑1 (MQ=255)
                                                                                                                                        accacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACTGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAttt          >  2:903663/1‑123 (MQ=255)
                                                                                                                                        accacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAttt          <  1:903663/123‑1 (MQ=255)
                                                                                                                                          cacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgtgt                    >  1:147208/1‑111 (MQ=255)
                                                                                                                                          cacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgtgt                    <  2:147208/111‑1 (MQ=255)
                                                                                                                                            cgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCg                  >  2:774629/1‑111 (MQ=255)
                                                                                                                                            cgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCg                  <  1:774629/111‑1 (MQ=255)
                                                                                                                                            cgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtag        >  3:845614/1‑121 (MQ=255)
                                                                                                                                            cgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGAGTCGTTGTATTtag        <  4:845614/121‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                         
TCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACA  >  NC_000913/2440720‑2440986

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: