Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,216,702 T→C 47.5% Q677Q (CAA→CAG prc ← carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,216,7020TC47.5% ‑4.7 / 55.2 40Q677Q (CAA→CAGprccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (7/14);  new base C (6/13);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-02

CGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAT  >  Exported/1216641‑1216769
                                                             |                                                                   
cGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCcgcg                                                            >  1:940825/1‑71 (MQ=255)
 gCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCgg                                                           <  1:3590293/71‑1 (MQ=255)
     ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGt                                                         <  1:1536877/69‑1 (MQ=255)
     ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGt                                                         <  1:772533/69‑1 (MQ=255)
     ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTc                                                        <  1:2238637/70‑1 (MQ=255)
     ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTc                                                        <  1:1796495/70‑1 (MQ=255)
      tAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTg                                                      >  1:1481072/1‑71 (MQ=255)
        aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGc                                                    <  1:776386/71‑1 (MQ=255)
          aaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCtt                                                  >  1:651032/1‑71 (MQ=255)
            aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCtttt                                                <  1:43951/71‑1 (MQ=255)
            aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCtttt                                                <  1:2059974/71‑1 (MQ=255)
            aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCtttt                                                <  1:561759/71‑1 (MQ=255)
            aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCCGGCtttt                                                <  1:708869/71‑1 (MQ=255)
             tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCttttt                                               <  1:84240/71‑1 (MQ=255)
             tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCtttt                                                >  1:2922181/1‑70 (MQ=255)
               aGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCa                                             <  1:558096/71‑1 (MQ=255)
                   acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCt                                         <  1:2962992/71‑1 (MQ=255)
                   acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCt                                         >  1:1227180/1‑71 (MQ=255)
                     aaTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTc                                       >  1:2442573/1‑71 (MQ=255)
                      aTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTc                                                        <  1:3251056/53‑1 (MQ=255)
                       tttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc                                     >  1:1789675/1‑71 (MQ=255)
                        ttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCa                                             <  1:1622715/62‑1 (MQ=255)
                         aattGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCcg                                                              >  1:3532698/3‑44 (MQ=255)
                           ttGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTTTGGCTTTTTCAAGCTTCGCCa                                   <  1:1992505/69‑1 (MQ=255)
                              ttccTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATc                               <  1:1477415/68‑1 (MQ=255)
                                ccTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTc                                                        <  1:191076/43‑1 (MQ=255)
                                ccTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTc                                                        <  1:2703118/43‑1 (MQ=255)
                                    attgatATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTg                          <  1:1005589/69‑1 (MQ=255)
                                     ttgatATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                       >  1:1038577/1‑71 (MQ=255)
                                          atTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt                   <  1:236918/70‑1 (MQ=255)
                                                 gACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct           >  1:1087673/1‑71 (MQ=255)
                                                 gACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtc            <  1:2477758/70‑1 (MQ=255)
                                                  aCGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATc                               <  1:3252647/50‑1 (MQ=255)
                                                  aCGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtctc          <  1:601337/71‑1 (MQ=255)
                                                   cGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCa         <  1:2880456/71‑1 (MQ=255)
                                                    gggAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGa                             >  1:2639630/1‑50 (MQ=255)
                                                      gAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc       >  1:2484004/1‑70 (MQ=255)
                                                      gAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc       >  1:679612/1‑70 (MQ=255)
                                                       aGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATcc      <  1:2417351/70‑1 (MQ=255)
                                                          gggCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAt  <  1:2483511/71‑1 (MQ=255)
                                                             |                                                                   
CGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAT  >  Exported/1216641‑1216769

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