Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 309,358 C→T *223* (TAG→TAA)  ecpB ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913309,3580CT100.0% 127.3 / NA 40*223* (TAG→TAA) ecpBECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

TGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGC  >  NC_000913/309226‑309492
                                                                                                                                     |                                                                                                                                      
tgctgcGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTc                                                                                                                                  >  2:975995/1‑139 (MQ=255)
     cgcACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACggg                                                                                                                             >  2:491781/1‑139 (MQ=255)
         cTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATg                                                                                                                         <  2:795874/139‑1 (MQ=255)
                      cgcTGGCATCGGGCTGAACGAACAAAAAGCCCAGGCCGAGGGC‑AACTGGGCCTTCCGGCCCGGGGAAAACCCGCCTAAAGGGAATGTTGTTAATCCTCTGCACAACCTAT‑TAATTCAAGGGAATGGAATTAATACCCTg                                                                                                           >  1:35028/1‑139 (MQ=255)
                      cgcTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCt                                                                                                            >  1:938674/1‑139 (MQ=255)
                        cTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGc                                                                                                          >  1:147621/1‑139 (MQ=255)
                                     gaacgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAccc                                                                                             <  1:449270/139‑1 (MQ=255)
                                     gaacgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGCGCAAccc                                                                                             <  2:938674/139‑1 (MQ=255)
                                      aacgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCg                                                                                            >  1:493567/1‑139 (MQ=255)
                                        cgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTc                                                                                          <  2:49738/139‑1 (MQ=255)
                                                aaGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                               >  2:587488/1‑126 (MQ=255)
                                                aaGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                               <  1:587488/126‑1 (MQ=255)
                                                     cATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGc                                                                             >  1:495311/1‑139 (MQ=255)
                                                              ggCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCAc                                                                    <  1:695471/139‑1 (MQ=255)
                                                                                        agaACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCg                                          <  1:645963/139‑1 (MQ=255)
                                                                                          aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCa                                                                                                        >  1:120089/1‑75 (MQ=255)
                                                                                          aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCa                                                                                                        <  2:120089/75‑1 (MQ=255)
                                                                                          aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                                        >  1:530100/1‑139 (MQ=255)
                                                                                          aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                                        >  2:194070/1‑139 (MQ=255)
                                                                                          aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                                        >  2:34414/1‑139 (MQ=255)
                                                                                            ccGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATg                                                                                                                         >  1:121861/1‑56 (MQ=255)
                                                                                            ccGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATg                                                                                                                         <  2:121861/56‑1 (MQ=255)
                                                                                            ccGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCAt                                      <  2:493567/139‑1 (MQ=255)
                                                                                                cgtAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                                  <  2:495311/139‑1 (MQ=255)
                                                                                                   aaaGGCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtag                               <  1:194070/139‑1 (MQ=255)
                                                                                                       gCATTGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAt                           <  1:34414/139‑1 (MQ=255)
                                                                                                          ttGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtct                        >  2:858960/1‑139 (MQ=255)
                                                                                                           tGTCGTTATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtctc                       >  1:257904/1‑139 (MQ=255)
                                                                                                                ttATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTac                  <  1:858960/139‑1 (MQ=255)
                                                                                                                 tATCCATCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                 >  2:704618/1‑139 (MQ=255)
                                                                                                                       tCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCttt                                                                                      >  2:143262/1‑64 (MQ=255)
                                                                                                                       tCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCttt                                                                                      <  1:143262/64‑1 (MQ=255)
                                                                                                                       tCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                 >  2:891483/1‑133 (MQ=255)
                                                                                                                       tCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                 <  1:891483/133‑1 (MQ=255)
                                                                                                                       tCTGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt           >  2:844779/1‑139 (MQ=255)
                                                                                                                         tGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATc                                                                                <  2:334767/68‑1 (MQ=255)
                                                                                                                         tGCAC‑AGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATc                                                                                >  1:334767/1‑68 (MQ=255)
                                                                                                                           cac‑aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATc                                     <  2:214567/109‑1 (MQ=255)
                                                                                                                           cac‑aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATc                                     >  1:214567/1‑109 (MQ=255)
                                                                                                                             c‑aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                            >  1:294818/1‑100 (MQ=255)
                                                                                                                             c‑aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                            <  2:294818/100‑1 (MQ=255)
                                                                                                                                 ctatttatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggc  >  2:411165/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCAC‑AGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGC  >  NC_000913/309226‑309492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑