| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 2,709,297 | T→A | 28.0% | K544* (AAA→TAA) | USA300HOU_2562 ← | possible acetyltransferase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,709,297 | 0 | T | A | 28.0% | 36.6 / 16.1 | 25 | K544* (AAA→TAA) | USA300HOU_2562 | possible acetyltransferase | 
| Reads supporting (aligned to +/- strand): ref base T (6/12); new base A (2/5); total (8/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01 | |||||||||||
GCTTTACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTTAATTGATCTTTCGTAAATGCCCCATTTGTACCAA  >  CP000730/2709170‑2709416                                                                                                                               |                                                                                                                       gcTTTACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTAcc                                                                                                             <  2:328857/140‑1 (MQ=255)  tttACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCtt                                                                                                           <  2:71468/140‑1 (MQ=255)      ctacCTGCATATTATAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGATTTTCAGCTGACCCACCAGTTAATGTATTTATTCTACCTTCATa                                                                                                       >  2:73819/1‑140 (MQ=255)       tacCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATaa                                                                                                      >  2:240550/1‑140 (MQ=255)           tGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCACAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCACAAtctc                                                                                                  <  2:196380/140‑1 (MQ=255)                  tCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTAc                                                                                           <  1:37346/140‑1 (MQ=255)                    tAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTActc                                                                                         <  1:200119/140‑1 (MQ=255)                    tAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTActc                                                                                         <  1:254030/140‑1 (MQ=255)                         gAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATCATCTCTAGGTACTCTAATa                                                                                    <  1:213580/140‑1 (MQ=255)                                      tGCAAAGTATTCCGGATGACCTGCAGAAGCATTATACCAATCGACTGGATGTACGTTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAatata                                                                       >  2:162535/1‑140 (MQ=255)                                          aaGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTc                                                                   <  2:321078/140‑1 (MQ=255)                                          aaGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTc                                                                   <  2:97326/140‑1 (MQ=255)                                          aaGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTc                                                                   <  2:240710/140‑1 (MQ=255)                                            gTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTg                                                                 <  2:279920/140‑1 (MQ=255)                                              aTTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCt                                                               <  1:314681/140‑1 (MQ=255)                                               ttCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCtt                                                              <  2:22003/140‑1 (MQ=255)                                                        aCCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAAc                                                     <  2:196108/140‑1 (MQ=255)                                                             cAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATcc                                                >  2:52819/1‑140 (MQ=255)                                                                taGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAAt                                             <  1:103883/139‑1 (MQ=255)                                                                   cgttATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATTGt                                          <  2:5183/138‑1 (MQ=255)                                                                   cTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATATATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGt                                          <  2:195523/140‑1 (MQ=255)                                                                        taCCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACGCTAATAGATACTAAATATATGTCTGCTTTTCTAAAACTATCCACTAGTTCAtt                                     >  2:188335/1‑140 (MQ=255)                                                                                           cATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTTAATTGATCTTTCGTAAAt                  >  1:185439/1‑140 (MQ=255)                                                                                                  gcgcTTTTCAGCTGCCTCATGAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTAAATTGATCTTTCGTAAATGCCCCAt           >  1:215846/1‑140 (MQ=255)                                                                                                           aGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAATTCATTTAATTGATCTTTCGTAAATCCCCCATTTGTACCaa  >  1:23757/1‑140 (MQ=255)                                                                                                                               |                                                                                                                       GCTTTACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTTAATTGATCTTTCGTAAATGCCCCATTTGTACCAA  >  CP000730/2709170‑2709416 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |