Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 309,486 C→T G181S (GGT→AGT)  ecpB ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913309,4860CT98.2% 169.1 / ‑4.1 55G181S (GGT→AGT) ecpBECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/25);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

AGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGT  >  NC_000913/309352‑309612
                                                                                                                                      |                                                                                                                              
aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTg                                                                                                                            >  1:139972/1‑139 (MQ=255)
            cACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAgg                                                                                                                <  1:1680937/139‑1 (MQ=255)
              cGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCa                                                                                                              <  1:1634980/139‑1 (MQ=255)
                           aTCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCg                                                                                                 >  2:1370290/1‑139 (MQ=255)
                                ccTGCCAAAGTGCGACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCg                                                                                                 >  2:188189/1‑134 (MQ=255)
                                ccTGCCAAAGTGCGACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCg                                                                                                 <  1:188189/134‑1 (MQ=255)
                                ccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGAt                                                                                            >  2:481032/1‑139 (MQ=255)
                                   gCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCg                                                                                         >  1:736740/1‑139 (MQ=255)
                                             aaCCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGc                                                                                >  1:1005632/1‑138 (MQ=255)
                                             aaCCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGc                                                                                <  2:1005632/138‑1 (MQ=255)
                                                 cGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTc                                                                           <  1:110357/139‑1 (MQ=255)
                                                    ccTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTg                                                                        >  1:711306/1‑139 (MQ=255)
                                                     cTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGt                                                                       <  2:1606225/139‑1 (MQ=255)
                                                              ccGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTggcg                                                                                                                         >  2:1309583/1‑80 (MQ=255)
                                                              ccGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTggcg                                                                                                                         <  1:1309583/80‑1 (MQ=255)
                                                               cGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCa                                                                                       <  2:88675/113‑1 (MQ=255)
                                                               cGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCa                                                                                       >  1:88675/1‑113 (MQ=255)
                                                                  aGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCt                                                                         >  1:194850/1‑124 (MQ=255)
                                                                  aGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCt                                                                         <  2:194850/124‑1 (MQ=255)
                                                                    tgtCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCa                                                                               <  2:1300904/116‑1 (MQ=255)
                                                                    tgtCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCa                                                                               >  1:1300904/1‑116 (MQ=255)
                                                                      tCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCa                                                                                                              >  2:1645342/1‑83 (MQ=255)
                                                                      tCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCa                                                                                                              <  1:1645342/83‑1 (MQ=255)
                                                                        cACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTg                                                    >  2:108476/1‑139 (MQ=255)
                                                                         aCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGc                                                                                                               <  1:56711/79‑1 (MQ=255)
                                                                         aCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGc                                                                                                               >  2:56711/1‑79 (MQ=255)
                                                                                tAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCa                                                                                       >  2:1468479/1‑96 (MQ=255)
                                                                                tAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCa                                                                                       <  1:1468479/96‑1 (MQ=255)
                                                                                        cgcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCg                                                                                                       <  1:825663/72‑1 (MQ=255)
                                                                                        cgcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCg                                                                                                       >  2:825663/1‑72 (MQ=255)
                                                                                        cgcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGc                                                                                <  2:190635/95‑1 (MQ=255)
                                                                                        cgcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGc                                                                                >  1:190635/1‑95 (MQ=255)
                                                                                         gcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTggcg                                                                                                                         <  1:513880/53‑1 (MQ=255)
                                                                                         gcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTggcg                                                                                                                         >  2:513880/1‑53 (MQ=255)
                                                                                         gcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCAt                                                                              >  1:1022689/1‑96 (MQ=255)
                                                                                         gcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCAt                                                                              <  2:1022689/96‑1 (MQ=255)
                                                                                         gcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTc                                                                           >  2:563034/1‑99 (MQ=255)
                                                                                         gcgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTc                                                                           <  1:563034/99‑1 (MQ=255)
                                                                                          cgACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAg                                  <  2:971663/139‑1 (MQ=255)
                                                                                               tGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGc                             <  1:481032/139‑1 (MQ=255)
                                                                                                  ccGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGCTCCTCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTg                          <  2:260862/139‑1 (MQ=255)
                                                                                                   cGGCATCAGGTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGa                         >  2:1487304/1‑139 (MQ=255)
                                                                                                           ggTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGa                                           >  2:1271092/1‑113 (MQ=255)
                                                                                                           ggTAGTAATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGa                                           <  1:1271092/113‑1 (MQ=255)
                                                                                                                taATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCg                                                                                                 >  2:449524/1‑54 (MQ=255)
                                                                                                                taATTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCg                                                                                                 <  1:449524/54‑1 (MQ=255)
                                                                                                                 aaTTCTCTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCgacgac           <  2:550777/139‑1 (MQ=255)
                                                                                                                    tctcTTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAg        >  2:1345361/1‑139 (MQ=255)
                                                                                                                     ctctTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAgg                                                                                              >  1:54569/1‑52 (MQ=255)
                                                                                                                     ctctTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAgg                                                                                              <  2:54569/52‑1 (MQ=255)
                                                                                                                     ctctTTACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGa       >  2:875098/1‑139 (MQ=255)
                                                                                                                          tACACTCCTTACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGt  >  2:430424/1‑139 (MQ=255)
                                                                                                                                  ttACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGt                                                                                   <  2:39657/50‑1 (MQ=255)
                                                                                                                                  ttACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGt                                                                                   >  1:39657/1‑50 (MQ=255)
                                                                                                                                   tACTGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCg                                                                                 >  1:1352379/1‑51 (MQ=255)
                                                                                                                                      |                                                                                                                              
AGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGT  >  NC_000913/309352‑309612

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