| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 990,149 | C→G | 33.0% | S14C (TCC→TGC) | oppD1 → | oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 990,149 | 0 | C | G | 33.0% | 36.6 / 25.1 | 31 | S14C (TCC→TGC) | oppD1 | oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein | 
| Reads supporting (aligned to +/- strand): ref base C (9/11); new base G (5/5); total (14/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 = 9.67e-01 | |||||||||||
TTTTACTTATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAATTACAAAATTA  >  CP000730/990036‑990282                                                                                                                 |                                                                                                                                     ttttACTTATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTGGATATTACAGGAGGGGAAGAgc                                                                                                             >  2:22887/1‑140 (MQ=255)   tACTTATTTAGTGAAAGAATACGAGATGCAATGGAACCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTgcagg                                                                                                          <  1:11532/140‑1 (MQ=255)    acctATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCACATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTgcaggc                                                                                                         <  2:316772/137‑1 (MQ=255)      ttATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTgcaggcag                                                                                                       <  1:198413/140‑1 (MQ=255)         tttAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGAAGAGGAAGTGCAGGCAGTga                                                                                                    >  2:185805/1‑140 (MQ=255)              tGATGGATTACGTTATTCATATGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGTAAGTGCAGGCAGTGAGAGGc                                                                                               <  2:70906/140‑1 (MQ=255)              tGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGc                                                                                               <  2:48038/140‑1 (MQ=255)                      tACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGAttt                                                                                       <  2:248503/140‑1 (MQ=255)                                 tttGATCTGAAAATGCGTAAATAAAAAGGGGGCATAGCAGATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACa                                                                            <  1:34620/140‑1 (MQ=255)                                 ttggATCCGATAATGCGTAAATAGAGAGGGGGTATTGCATATGACTGAATGAATATTAGAAGTAAAGGATTTGCATGTTTCCTTTGATATTACAGCAGGGGGAGTGCAGGCAGTGAGAGGCGTAGTTTTTTATTTGAACa                                                                            <  2:170438/137‑1 (MQ=255)                                  ttGATTCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACaa                                                                           <  1:322610/140‑1 (MQ=255)                                      tCCTAAAATATGTAAATAAAAAGTGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTACTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAggg                                                                       <  1:138899/140‑1 (MQ=255)                                            aaTGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGTTATTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACa                                                                 <  1:245164/140‑1 (MQ=255)                                             aTGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACGAAGGGGACACAt                                                                >  2:222050/1‑140 (MQ=255)                                             aTGCGTAAATAAAAAGGGGGCATAGCATATGACAGAAAGAAGATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGACGTGCAGGCAGTGAGAGGCGTAGATTTTTAGTTGAACAAAGGGGAAACAt                                                                >  2:12861/1‑140 (MQ=255)                                              tGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACAtt                                                               >  1:298793/1‑140 (MQ=255)                                                cGTAAAAAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTgg                                                             <  2:215178/140‑1 (MQ=255)                                                agaaaaaaaaaAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTgg                                                             <  2:63157/137‑1 (MQ=255)                                                  taaataaaAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCa                                                           <  1:43180/140‑1 (MQ=255)                                                         aaaGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAAttgttg                                                    >  1:336754/1‑140 (MQ=255)                                                                        tatGACTGAAATAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAAtcaggttcag                                     <  2:315482/140‑1 (MQ=255)                                                                           taCTGTAAGAATATTAGAGGTAAATGATGTGCATGTTTCCTTTGATATTACAGCAGGGGAAGGGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGCATCAGGTTCAGGTa                                  <  2:265530/139‑1 (MQ=255)                                                                              tGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAAt                               >  1:306211/1‑140 (MQ=255)                                                                               gAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATc                              >  2:295592/1‑140 (MQ=255)                                                                               gAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTAAGGTAAATc                              >  2:295525/1‑140 (MQ=255)                                                                                     atatTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTaac                        <  1:14715/140‑1 (MQ=255)                                                                                        ttAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTaacaac                     >  1:58468/1‑140 (MQ=255)                                                                                             aGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAAAAACAAAAg                >  1:128748/1‑140 (MQ=255)                                                                                                 aaTGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGGAACATTGGCAATTGATGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAAt            >  1:68506/1‑140 (MQ=255)                                                                                                     aTTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAATTACa        >  1:180373/1‑140 (MQ=255)                                                                                                           aTGTTTCCTTTGATATTACAGCAGCGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGGAACAACAAAAGAAATTACAAAatta  >  1:226145/1‑140 (MQ=255)                                                                                                                 |                                                                                                                                     TTTTACTTATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAATTACAAAATTA  >  CP000730/990036‑990282 | 
| 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: ‑ |