Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 155,463 A→T 100% *247K (TAA→AAA)  ecpD ← predicted periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048155,4630AT100.0% 204.2 / NA 67*247K (TAA→AAA) ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (37/30);  total (37/30)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCA  >  AP009048/155396‑155527
                                                                   |                                                                
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                               >  1:855836/1‑71 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                             <  1:4312206/71‑1 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                             <  1:436102/71‑1 (MQ=255)
    atAATTTTTAGTATATTCTATAGTCACTCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                            <  1:3884737/70‑1 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                          <  1:3399208/71‑1 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                          <  1:2476327/71‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                     >  1:4301477/1‑71 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                     >  1:4296712/1‑71 (MQ=255)
           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                    >  1:3258360/1‑71 (MQ=255)
              gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                  <  1:1071299/70‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                              >  1:1122141/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                               <  1:2664992/70‑1 (MQ=255)
                    tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                           >  1:2042693/1‑71 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                <  1:3618663/65‑1 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                          >  1:1757686/1‑71 (MQ=255)
                      tataGTCACGCTATGCTTCCTGAGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                         >  1:868231/1‑71 (MQ=255)
                       ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                        >  1:3723726/1‑71 (MQ=255)
                       ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                        >  1:3718843/1‑71 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCtt                                                 >  1:1936563/1‑56 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                   >  1:961145/1‑70 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                   >  1:3288745/1‑70 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                   >  1:1445224/1‑70 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                   <  1:3774077/70‑1 (MQ=255)
                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                 <  1:1713433/71‑1 (MQ=255)
                               gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                 <  1:4255164/70‑1 (MQ=255)
                                 tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCtt                                                 >  1:3582009/1‑52 (MQ=255)
                                 tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                <  1:2146594/53‑1 (MQ=255)
                                 tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                               >  1:2010979/1‑70 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                             >  1:1125650/1‑71 (MQ=255)
                                   tGCTTCCTGCGGATTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                              <  1:1442479/69‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                          <  1:1391160/71‑1 (MQ=255)
                                       tCCTGCGGAGTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                        >  1:2305251/1‑71 (MQ=255)
                                       tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                            <  1:2220210/51‑1 (MQ=255)
                                       tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                        >  1:2733243/1‑71 (MQ=255)
                                       tCCTGCGGAATTATATCCTTGACTGATTTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTAAt                        >  1:4535761/1‑71 (MQ=255)
                                         cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                    >  1:287045/1‑41 (MQ=255)
                                          tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                         <  1:380520/51‑1 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                            >  1:4109431/1‑62 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                <  1:4567077/71‑1 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                <  1:1679128/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa               <  1:1465146/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa               <  1:3707923/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                >  1:3041243/1‑70 (MQ=255)
                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                >  1:842511/1‑36 (MQ=255)
                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                         <  1:2624627/59‑1 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                             >  1:1166308/1‑54 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                           >  1:4208459/1‑56 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg            <  1:2679727/70‑1 (MQ=255)
                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                            <  1:2240856/36‑1 (MQ=255)
                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                <  1:1430545/64‑1 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt       <  1:3726580/71‑1 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt       <  1:2947965/71‑1 (MQ=255)
                                                        cTTGCCAGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                 >  1:600946/1‑45 (MQ=38)
                                                         ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt       >  1:1936274/1‑70 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt     >  1:989552/1‑71 (MQ=255)
                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg    >  1:4046422/1‑71 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                 >  1:1679775/1‑41 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt     <  1:2676857/69‑1 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg    <  1:2495208/70‑1 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg    <  1:4000317/70‑1 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg    <  1:3128786/70‑1 (MQ=255)
                                                              tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa              >  1:3970029/1‑58 (MQ=255)
                                                              tGATTTCAGCCTGGAATTACCTTCAATTGCGCCACCAAAGTTATTAATGGCGTAAAAATGCACTTTTg    >  1:1854048/1‑68 (MQ=255)
                                                               gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                           >  1:3185408/1‑44 (MQ=255)
                                                               gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCa  >  1:1277616/1‑69 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                        >  1:4108358/1‑46 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt     >  1:1370398/1‑65 (MQ=255)
                                                                   |                                                                
GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCA  >  AP009048/155396‑155527

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