Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 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
*W3110S.gb155,4630AT98.7% 225.6 / ‑3.2 76coding (739/741 nt)ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (39/36);  minor base . (1/0);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155395‑155528
                                                                    |                                                                 
tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:1480890/71‑1 (MQ=255)
    gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 >  1:781906/1‑67 (MQ=255)
    gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                              >  1:2847944/1‑70 (MQ=255)
    gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             >  1:2098403/1‑71 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                          <  1:3490139/70‑1 (MQ=255)
             tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                     <  1:3692148/70‑1 (MQ=255)
             tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                     <  1:528330/70‑1 (MQ=255)
              aGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                   >  1:3291702/1‑71 (MQ=255)
                tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:3606583/71‑1 (MQ=255)
                tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:1661293/71‑1 (MQ=255)
                 atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             >  1:410989/1‑58 (MQ=255)
                 atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             >  1:3148668/1‑58 (MQ=255)
                 atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCtt                                                  >  1:3769565/1‑69 (MQ=255)
                 atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:1482272/1‑70 (MQ=255)
                  tatTCTATAGTCAGGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:382017/1‑71 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2717097/1‑71 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:2827456/70‑1 (MQ=255)
                     tatataGTCACGCTATGCTTCCAGCGGAATTATATACTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                            >  1:2773512/3‑71 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:558804/1‑67 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2541247/1‑67 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                           >  1:2817505/1‑71 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                            >  1:705156/1‑70 (MQ=255)
                       tataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                          >  1:4273714/1‑71 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:2778699/1‑43 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacta                                    >  1:3126936/1‑68 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    >  1:3360605/1‑70 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                >  1:816867/1‑70 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              >  1:3488844/1‑71 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              >  1:3101312/1‑71 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:1658831/35‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:3940305/35‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:4616656/35‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:4667382/35‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:3347954/52‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                               <  1:3561369/69‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTCAGCCTGGCATTa                                                      >  1:1281135/1‑45 (MQ=38)
                                      cTTCCTGCGTAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:4727710/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                              >  1:3134217/1‑36 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:4258540/1‑35 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:4670762/1‑35 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                     >  1:3128552/1‑45 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:3562205/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:2521319/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:4022425/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:60438/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:2740474/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:4017508/71‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:3630605/71‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:3534967/71‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:3352385/71‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                    >  1:462809/1‑43 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                    >  1:4019902/1‑43 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          <  1:89495/69‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          <  1:2570704/69‑1 (MQ=255)
                                           tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:2734329/65‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:259036/61‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:2123101/61‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                 <  1:3489169/71‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                <  1:1340135/71‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                 >  1:1214595/1‑70 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                 >  1:2108781/1‑70 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:2383946/1‑60 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              >  1:669720/1‑38 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              >  1:2527121/1‑38 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGc            <  1:3094415/71‑1 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              >  1:1524696/1‑51 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc          >  1:1152275/1‑71 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                             <  1:1300719/50‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt        <  1:10016/71‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGACACCAAAGTCAt                             <  1:2483526/50‑1 (MQ=255)
                                                           tGCCTGATTTCAGACTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt      >  1:2885294/1‑71 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:2736292/42‑1 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:1261993/42‑1 (MQ=255)
                                                               tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaa                  <  1:2154279/55‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  <  1:1592175/69‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  <  1:637084/69‑1 (MQ=255)
                                                                    |                                                                 
TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155395‑155528

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