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% 234.9 / NA 77*247K (TAA→AAA) ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (40/37);  total (40/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  AP009048/155395‑155530
                                                                    |                                                                   
tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:3619228/71‑1 (MQ=255)
 gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:2603509/1‑71 (MQ=255)
  atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:2454417/1‑71 (MQ=255)
   tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                <  1:1735630/71‑1 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:2713475/1‑71 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                              >  1:4181576/1‑70 (MQ=255)
           tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                        >  1:4720342/1‑71 (MQ=255)
            ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:3488313/1‑71 (MQ=255)
            ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:2539098/1‑71 (MQ=255)
               gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                     <  1:5148682/70‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:3558132/1‑71 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:5274570/71‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:5250861/71‑1 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:2341397/1‑51 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:2684209/1‑71 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:2863302/1‑71 (MQ=255)
                       tataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                           >  1:2436367/1‑56 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                         >  1:2394295/1‑70 (MQ=255)
                            tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:3775053/1‑47 (MQ=255)
                            tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:577340/53‑1 (MQ=255)
                              aCGCTATGCTTCCTCCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:2156593/70‑1 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   <  1:2858375/68‑1 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:644598/1‑71 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                <  1:1173016/38‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                  <  1:3218119/52‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                               >  1:4287585/1‑71 (MQ=255)
                                     gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:664251/1‑52 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                              >  1:1007634/1‑38 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                      <  1:3220990/46‑1 (MQ=255)
                                       ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             <  1:445056/54‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:3145197/1‑36 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:5133697/1‑42 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            <  1:3254353/69‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                          >  1:222218/1‑71 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:2798870/39‑1 (MQ=255)
                                            gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          <  1:5170600/36‑1 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcc                                          >  1:3509645/1‑51 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:4702319/61‑1 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                      <  1:810039/71‑1 (MQ=255)
                                              ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:59977/54‑1 (MQ=255)
                                              ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                    >  1:1912152/1‑56 (MQ=255)
                                              ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                    >  1:5136727/1‑56 (MQ=255)
                                               gAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaa                    <  1:1773985/71‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:2399343/71‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaa                    >  1:3514838/1‑70 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:3143125/71‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   >  1:1142217/1‑70 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                            <  1:2061954/44‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa                 >  1:4741616/1‑71 (MQ=255)
                                                   tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa                 <  1:3689816/70‑1 (MQ=255)
                                                   tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1945107/68‑1 (MQ=255)
                                                   tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1951222/68‑1 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:225110/1‑37 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:325053/1‑37 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAAt                <  1:4762007/70‑1 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:4161112/1‑71 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                    >  1:3731518/1‑49 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            <  1:914243/57‑1 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               <  1:4158664/70‑1 (MQ=255)
                                                      atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:821079/52‑1 (MQ=255)
                                                      atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:2279649/66‑1 (MQ=255)
                                                        ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa                 <  1:3700940/65‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                             >  1:2781028/1‑52 (MQ=255)
                                                          ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                    <  1:3178076/44‑1 (MQ=255)
                                                          ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:769309/1‑45 (MQ=255)
                                                          ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:3841114/1‑45 (MQ=255)
                                                           tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc            >  1:2804666/1‑67 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGc                        <  1:4193189/54‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       >  1:3753358/1‑71 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       >  1:1482727/1‑71 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       >  1:3156161/1‑71 (MQ=255)
                                                             ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       <  1:416682/70‑1 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                      <  1:232771/54‑1 (MQ=255)
                                                               tGATTTCAGCCTGGTATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                         <  1:2563018/50‑1 (MQ=255)
                                                               tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1691217/45‑1 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                      >  1:4412292/1‑52 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:5221089/1‑71 (MQ=255)
                                                                    |                                                                   
TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  AP009048/155395‑155530

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: