Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,191,432 C→T 100% G72S (GGC→AGC)  yehC ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,191,4320CT96.1% 233.7 / 0.4 76G72S (GGC→AGC) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (29/44);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191335‑2191528
                                                                                                 |                                                                                                
aTTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtat                                                                                                >  1:1231779/1‑100 (MQ=255)
 tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  2:1593164/100‑1 (MQ=255)
 tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  1:428420/100‑1 (MQ=255)
 tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATAACAGGAACCTGAATTTTTTCTCGTGGTAATGACGTATCGCtatc                                                                                               >  2:1097081/1‑100 (MQ=255)
  ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              >  2:626617/1‑100 (MQ=255)
  ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              <  2:2181431/100‑1 (MQ=255)
    ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            <  1:538278/100‑1 (MQ=255)
      aCTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCaa                                                                                          <  1:1279180/100‑1 (MQ=255)
      aCTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCaa                                                                                          >  1:339719/1‑100 (MQ=255)
        ttGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                        <  1:1386478/100‑1 (MQ=255)
         tgcttcCCGGAATTTGCCCCTATTTTTGCCATTGGTGGTGTTAACATGAAAGCAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCGTCAATcc                                                                                       <  2:641181/100‑1 (MQ=255)
           ctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCAc                                                                                     <  2:191697/100‑1 (MQ=255)
            tgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACg                                                                                    <  1:2076893/100‑1 (MQ=255)
             gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  1:2126664/1‑100 (MQ=255)
               ccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCt                                                                                 >  1:2249069/1‑100 (MQ=255)
                 ggAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgc                                                                               >  1:1825452/1‑100 (MQ=255)
                  gAATTTGCCCCTATTTTTGCCACTGGTGGCGTGAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgca                                                                              >  1:2279409/1‑100 (MQ=255)
                  gAATTTGCCCCTATTTTTGCCACTGGTGGCGTAAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgca                                                                              >  1:2278615/1‑100 (MQ=255)
                    aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  1:48788/100‑1 (MQ=255)
                    aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgca                                                                              <  1:206493/98‑1 (MQ=255)
                     tttGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagca                                                                           <  1:20287/100‑1 (MQ=255)
                       tggtgcTAGTTTTGCGACTGGTGGCTTTAACATGAAAGGAACCTGAAGTTTTTCTGGTGGTAATGACGTATAGCTATTATCAATCCACGCCTGCagcagc                                                                         <  2:1037062/95‑1 (MQ=255)
                        gCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCg                                                                        >  2:241290/1‑100 (MQ=255)
                         ccccTATTTTTGTCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       <  1:273436/100‑1 (MQ=255)
                           ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                     <  1:957227/100‑1 (MQ=255)
                            cTATTTTTGCCACTGTTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaaga                                                                    >  2:1124431/1‑100 (MQ=255)
                            cTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaaga                                                                    <  1:1415364/100‑1 (MQ=255)
                            cTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                     <  1:188964/99‑1 (MQ=255)
                             tATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                     <  1:818123/98‑1 (MQ=255)
                               tttttGCAACTGGTTGCGTTAACATGAAAGTAACCTGAATTTTTTGTGGTGGTAATGACGTATCGCTATCATCAATCCACGTGTGCAGCAGCGAAGAAc                                                                  <  2:682591/99‑1 (MQ=255)
                                ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                                <  1:267892/100‑1 (MQ=255)
                                  ttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGtt                                                              <  2:775836/100‑1 (MQ=255)
                                  ttGCCACTGGTCGCGTTAACATGAAAGGCACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCCATCCACGCCTGCAGCAGCGAAGAACGGtt                                                              <  2:1795088/100‑1 (MQ=255)
                                   tGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGttt                                                             <  1:1864437/100‑1 (MQ=255)
                                    gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                            >  2:1980048/1‑100 (MQ=255)
                                    gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                            >  2:360208/1‑100 (MQ=255)
                                       aCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCt                                                         <  1:1643052/100‑1 (MQ=255)
                                         tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCt                                                         <  2:955272/98‑1 (MQ=255)
                                         tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTg                                                        <  1:743698/99‑1 (MQ=255)
                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                       >  1:1349919/1‑99 (MQ=255)
                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      <  2:966446/100‑1 (MQ=255)
                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      >  1:2256577/1‑100 (MQ=255)
                                           gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                     >  2:1991049/1‑100 (MQ=255)
                                            tggCGTTAACATGAAAGCAACCTGAATTTTTTCTGGTGGTAATGCCGTATCGCTATCATCAATCCACGCCTGCAGCCGCGAAGAACGGTTTCCCTGGTtc                                                    >  2:344548/1‑100 (MQ=255)
                                               cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  2:361463/1‑100 (MQ=255)
                                               cGTTAACATGAAAGGAACCTGAATTATTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  2:1150282/1‑100 (MQ=255)
                                                  tAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAAc                                              >  2:1315776/1‑100 (MQ=255)
                                                   aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             <  2:2256551/100‑1 (MQ=255)
                                                     cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAAc                                              >  1:124768/1‑97 (MQ=255)
                                                        gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                        <  2:752450/100‑1 (MQ=255)
                                                        gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                        <  1:361466/100‑1 (MQ=255)
                                                         aaaGTATCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCa                                       <  1:663266/100‑1 (MQ=255)
                                                           agtaaCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATc                                     <  2:993366/97‑1 (MQ=255)
                                                           aGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATc                                     <  2:17327/100‑1 (MQ=255)
                                                                  tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              >  2:1242129/1‑100 (MQ=255)
                                                                     aTTTTTTTTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  1:789915/1‑100 (MQ=255)
                                                                     aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTa                            <  2:745528/99‑1 (MQ=255)
                                                                     aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  2:172268/1‑100 (MQ=255)
                                                                       tttttCTGGTGGTAATGACGTATCGCTATCATCATTCCACGCCTGCAGCAGCGAAGCACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                         <  2:1031947/100‑1 (MQ=255)
                                                                       tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  1:1237216/1‑98 (MQ=255)
                                                                        ttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtttt                        >  1:486127/1‑100 (MQ=255)
                                                                         tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                       <  2:1825431/100‑1 (MQ=255)
                                                                          ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                      <  2:2154128/100‑1 (MQ=255)
                                                                           tCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                     <  1:2125316/100‑1 (MQ=255)
                                                                            cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGc                    <  1:270064/100‑1 (MQ=255)
                                                                              ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  >  2:1113792/1‑100 (MQ=255)
                                                                               gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTACACCATCACTTCTTTATTTTCTGCCgg                 <  2:673218/100‑1 (MQ=255)
                                                                                tggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCggg                <  2:2377081/100‑1 (MQ=255)
                                                                                  gtAATGACGTATCGCTATCATCACTCCTCGCCTGCAGCATCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTa              <  2:1654305/100‑1 (MQ=255)
                                                                                  gtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTa              <  1:510663/100‑1 (MQ=255)
                                                                                  gtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAATTTCTGCCGGGTa              <  1:1609354/100‑1 (MQ=255)
                                                                                   tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa             >  2:698779/1‑100 (MQ=255)
                                                                                   tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTa              <  2:2250641/99‑1 (MQ=255)
                                                                                    aaTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaaa            <  1:1124441/100‑1 (MQ=255)
                                                                                    aaTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa             >  2:691734/1‑99 (MQ=255)
                                                                                     aTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat           >  1:2137873/1‑100 (MQ=255)
                                                                                      tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata          <  2:2328567/100‑1 (MQ=255)
                                                                                      tGACGTATCGCTATCATCAATCAACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata          >  2:1651152/1‑100 (MQ=255)
                                                                                          gTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgc      >  2:437434/1‑100 (MQ=255)
                                                                                          gTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgc      <  1:1242142/100‑1 (MQ=255)
                                                                                             tCGCTATGATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  2:1282126/100‑1 (MQ=255)
                                                                                              cGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGt  >  2:2381666/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191335‑2191528

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