| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 2,013,179 | A→T | 30.2% | Q19H (CAA→CAT) | yfhP → | possible metal‑dependent hydrolase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,013,179 | 0 | A | T | 30.2% | 40.8 / 21.7 | 33 | Q19H (CAA→CAT) | yfhP | possible metal‑dependent hydrolase | 
| Reads supporting (aligned to +/- strand): ref base A (14/9); new base T (4/6); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01 | |||||||||||
CATGGTTAGTATGTAAGTTAATTTAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTCGCATCATAGAGGTATCACGCA  >  CP000730/2013039‑2013313                                                                                                                                            |                                                                                                                                      cATGGTTAGTATGTAAGTTAATTTAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCaa                                                                                                                                        >  2:16757/1‑141 (MQ=255)         tATGTAAGTTAAATTAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCa                                                                                                                               <  1:152860/141‑1 (MQ=255)                ttaaaTTTAGGTATTTGCAAAATTGGATTATAATAGTATATATAATATTATGAACTGAGCGGTCGGTGATGGACACTGCAACACTTATCGCAATTGGGGTGGGCCTTACATCACTTGCAACTCATGATCCAGCAATTGCtt                                                                                                                        <  2:122256/138‑1 (MQ=255)                gTTAATTTAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCGACTCAAGATCCAGCAATGGCtt                                                                                                                        <  2:252009/141‑1 (MQ=255)                      ttAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGt                                                                                                                  >  2:8487/1‑141 (MQ=255)                      ttAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGt                                                                                                                  >  1:107280/1‑141 (MQ=255)                       tAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGtt                                                                                                                 >  2:127694/1‑141 (MQ=255)                        aGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAAATGATATAAACACGGCAAAACATATTTCAATTGGGGTAGAGCGTACAGCACTTGCAACTCAACATCCAGCAACAGGTTCTACGttt                                                                                                                >  1:241981/1‑141 (MQ=255)                             tttGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTGTACGTTTGGTGc                                                                                                           >  1:279753/1‑141 (MQ=255)                                  gAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGTTTGGTGCAACAg                                                                                                      <  2:100890/141‑1 (MQ=255)                                         gATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGTTTGGTGCAACACCTACAAc                                                                                               >  1:256289/1‑141 (MQ=255)                                         gATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAAc                                                                                               >  2:274030/1‑141 (MQ=255)                                                         ataataTTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTATGGAGCAACAGATACAACACTTGTcctgcgtaca                                                                               >  2:181203/1‑132 (MQ=255)                                                          taataTTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCAt                                                                              >  1:8981/1‑141 (MQ=255)                                                              tttattaAAGGAGTGAACGGATATGGACACTGCAACACATATCGCAATTGTGGTGGGCCTTACTGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAATCCTTATTGTTGGTTCATTAAt                                                                          <  2:236365/140‑1 (MQ=255)                                                              attatGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGCGGTGGGCCTTACGGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAAGCCTTATTGTTGGTTCATTAAt                                                                          <  2:236268/141‑1 (MQ=255)                                                                   gAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAAGTCCTg                                                                     <  1:346618/141‑1 (MQ=255)                                                                         aGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATgggg                                                               >  2:345295/1‑141 (MQ=255)                                                                          gTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGa                                                              <  1:47591/141‑1 (MQ=255)                                                                                 gATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTTATGGGGATACTGtt                                                       >  1:11073/1‑141 (MQ=255)                                                                                 gATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAACAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCCTTGGGTAATTAATTCATGCTCGGCATACTGtt                                                       >  2:60163/1‑141 (MQ=255)                                                                                   tatGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCt                                                     <  1:8487/141‑1 (MQ=255)                                                                                         cacTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCttaaatt                                               >  2:27724/1‑141 (MQ=255)                                                                                             gCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAg                                           >  1:258800/1‑141 (MQ=255)                                                                                                  acaTATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACaa                                      <  1:157717/141‑1 (MQ=255)                                                                                                      aaccataTTGGGTTGTGCCTTCTAGCACTTGCAACTCATCATCCAGCAATGGCTTCGACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCa                                  <  2:329471/135‑1 (MQ=255)                                                                                                        cGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAAc                                <  1:201662/141‑1 (MQ=255)                                                                                                                 ggTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTc                       >  1:277653/1‑141 (MQ=255)                                                                                                                 ggTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTc                       >  1:47728/1‑141 (MQ=255)                                                                                                                   tGGGCCTTACAGCACTTGCAACTCATGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTCGc                     <  1:345295/141‑1 (MQ=255)                                                                                                                   tGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTCGc                     >  1:337994/1‑141 (MQ=255)                                                                                                                         ttACAGCACTTGCAACTCATGATCCAGCAACGGCTTCTACGTTGGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTTAGGGGGACACTGTTCTTAAATTAAAGGACAATGCAACATATATTTCGCATCATa               <  1:330936/141‑1 (MQ=255)                                                                                                                                      aaCTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGTACAATGCAACATATATTTCGCATCATAGAGGTATCACGCa  >  2:145359/1‑141 (MQ=255)                                                                                                                                      aaCTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTCGCATCATAGAGGTATCACGCa  <  1:187645/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                      CATGGTTAGTATGTAAGTTAATTTAGTTATTTGCGAAATTGGATTATAATAGTATATATAATATTATGAAATGAGTGAACTGATATGGACACTGCAACACATATCGCAATTGGGGTGGGCCTTACAGCACTTGCAACTCAAGATCCAGCAATGGCTTCTACGTTTGGTGCAACAGCTACAACCCTTATCGTTGGTTCATTAATTCCTGATGGGGATACTGTTCTTAAATTAAAGGACAATGCAACATATATTTCGCATCATAGAGGTATCACGCA  >  CP000730/2013039‑2013313 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |