Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571858 20,628:1 +G 100% coding (288/306 nt) pSAS24 → putative exported protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX57185820,6281.G100.0% 51.9 / NA 13K96K (AAA→AAGpSAS24putative exported protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (5/8);  total (5/8)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTATAAATTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAATATTCCTAGCTTTTTAAAATGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑G  >  BX571858/20532‑20652
                                                                                                 |                                                                                      
ttATAAATTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACaaaa                                                                                         <  2:323406/97‑1 (MQ=255)
ttATAAATTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACaaaa                                                                                         <  2:323428/97‑1 (MQ=255)
    aaaTTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACaaaagata                                                                                     <  2:153659/97‑5 (MQ=255)
                      tATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAAt                                                                   >  2:376810/1‑97 (MQ=255)
                       aTGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATg                                                                  <  1:273505/97‑1 (MQ=255)
                             cTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtga                                                            <  2:354609‑M2/97‑6 (MQ=255)
                                     gAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagt                                                    >  1:181728‑M2/1‑84 (MQ=255)
                                        tatGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgat                                                 >  2:102622‑M2/1‑81 (MQ=255)
                                                   aaaTCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaa                                      >  2:295215‑M2/1‑70 (MQ=255)
                                                      tCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAACGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgttataacttgaaaatg                                   <  1:124456‑M2/97‑31 (MQ=255)
                                                      tCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatg                                   <  1:364723‑M2/97‑31 (MQ=255)
                                                             gggAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgt                            <  2:346342‑M2/97‑38 (MQ=255)
                                                                    atGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaa                     <  2:8308‑M2/97‑45 (MQ=255)
                                                                                 aGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaagggaataaaaaca        <  2:335280‑M2/97‑58 (MQ=255)
                                                                                    aTACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacattt     <  1:87050‑M2/97‑61 (MQ=255)
                                                                                       cATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacattttag  >  1:215448‑M2/1‑34 (MQ=255)
                                                                                                 |                                                                                      
TTATAAATTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAATATTCCTAGCTTTTTAAAATGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑G  >  BX571858/20532‑20652

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