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% 160.9 / NA 38K96K (AAA→AAGpSAS24putative exported protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/22);  total (16/22)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAATATTCCTAGCTTTTTAAAATGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑G  >  BX571858/20538‑20652
                                                                                           |                                                                                                 
aTTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTc                                                                                             <  1:189731/98‑1 (MQ=255)
      ttttttAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCt                                                                                       <  1:127203/98‑1 (MQ=255)
             ggTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTaaa                                                                                <  2:166979/98‑1 (MQ=255)
               ttATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAAt                                                                              <  1:251758/98‑1 (MQ=255)
                 aTGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGa                                                                            >  1:141800/1‑98 (MQ=255)
                     tCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtg                                                                        <  1:91208‑M2/98‑5 (MQ=255)
                      ccTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtga                                                                       >  2:161036‑M2/1‑93 (MQ=255)
                         aCAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccg                                                                    <  2:6448‑M2/98‑9 (MQ=255)
                          cAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccga                                                                   <  2:219265‑M2/98‑10 (MQ=255)
                          cAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccga                                                                   <  1:85638‑M2/98‑10 (MQ=255)
                                aaTATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtga                                                             <  2:219211‑M2/98‑16 (MQ=255)
                                     gtgAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaa                                                        >  1:59388‑M2/1‑78 (MQ=255)
                                            tAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaa                                                 >  2:298675‑M2/1‑71 (MQ=255)
                                            tAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaa                                                 <  2:282287‑M2/98‑28 (MQ=255)
                                             aaaTCATAAAGGGAAGTCTGTGCGATCGATAGAATACATTACAAAAGATATTCCTGGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaa                                                <  2:118765‑M2/98‑29 (MQ=255)
                                              aaTCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaat                                               <  2:249191‑M2/98‑30 (MQ=255)
                                               aTCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatg                                              >  1:320384‑M2/1‑68 (MQ=255)
                                                 cATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtt                                            >  2:89425‑M2/1‑66 (MQ=255)
                                                   tAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttt                                          <  2:151337‑M2/98‑35 (MQ=255)
                                                       gggAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgta                                      >  2:237584‑M2/1‑60 (MQ=255)
                                                        ggAAATATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaa                                     >  2:171856‑M2/1‑59 (MQ=255)
                                                          aaaTATGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaa                                   <  2:319640‑M2/98‑42 (MQ=255)
                                                            atatGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaa                                 >  1:125589‑M2/1‑55 (MQ=255)
                                                            atatGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaa                                 <  1:3394‑M2/98‑44 (MQ=255)
                                                             tatGTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaa                                <  2:194765‑M2/98‑45 (MQ=255)
                                                             aatgTGCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaagcttgaaaatgtttttgtaaaaaag                                >  2:148205‑M2/2‑54 (MQ=255)
                                                                 tgCGTTCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggta                            <  1:296786‑M2/98‑49 (MQ=255)
                                                                     ttCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataa                        <  2:274857‑M2/98‑53 (MQ=255)
                                                                      tCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaa                       <  1:236021‑M2/98‑54 (MQ=255)
                                                                      tCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaa                       <  2:8349‑M2/98‑54 (MQ=255)
                                                                      tCGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaa                       <  2:206946‑M2/98‑54 (MQ=255)
                                                                       cGATAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaa                      >  1:157254‑M2/1‑44 (MQ=255)
                                                                          tAGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaaca                   >  2:130961‑M2/1‑41 (MQ=255)
                                                                           aGAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacat                  >  1:163013‑M2/1‑40 (MQ=255)
                                                                            gAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacatt                 >  1:205374‑M2/1‑39 (MQ=255)
                                                                            gAATACATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaaaattttaaaaaaaaaaggtaataaaaacatt                 >  2:121907‑M2/1‑39 (MQ=255)
                                                                                aCATTACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacattttag             >  2:141800‑M2/1‑35 (MQ=255)
                                                                                    tACAAAAGATATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacattttagacgg         <  2:69101‑M2/98‑68 (MQ=255)
                                                                                      aaaagatATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacattttagacggatg      <  1:205312‑M2/97‑71 (MQ=255)
                                                                                          aaTATTCCTAGCTTTTTAAAATGAggtgaccgatagtgataaaacttgaaaatgtttttgtaaaaaaaggtaataaaaacattttagacggatgtaac  >  2:166835‑M2/1‑24 (MQ=255)
                                                                                           |                                                                                                 
ATTGGATTTTTTAGGTTATGATCCTACAAAAGAATATGTGAAAGTAAATCATAAAGGGAAATATGTGCGTTCGATAGAATACATTACAAAAATATTCCTAGCTTTTTAAAATGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑G  >  BX571858/20538‑20652

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