Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,191,432 C→T 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,4320CT98.3% 201.9 / ‑4.2 60G72S (GGC→AGC) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (32/27);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

ATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCG  >  NC_000913/2191335‑2191525
                                                                                                 |                                                                                             
aTTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtat                                                                                             <  1:579576/100‑1 (MQ=255)
 tttttACTTGCTGCCCGGAATTTGCCACTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                            <  2:611926/100‑1 (MQ=255)
    ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                         <  1:182459/100‑1 (MQ=255)
    ttACTTGCTGCCCGGAATTTGCCCCTATTCTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAACTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                         <  2:947080/100‑1 (MQ=255)
     tACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCAtcatca                                                                                        <  1:261156/100‑1 (MQ=255)
        ttGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                     <  2:573462/100‑1 (MQ=255)
         tgctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATcc                                                                                    >  1:471461/1‑100 (MQ=255)
           ctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCAc                                                                                  <  1:1360426/100‑1 (MQ=255)
             gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                >  1:1095786/1‑100 (MQ=255)
             gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                >  1:81097/1‑100 (MQ=255)
               ccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCt                                                                              >  2:389911/1‑100 (MQ=255)
                     tttGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagca                                                                        <  2:574131/100‑1 (MQ=255)
                      ttGCCCCTATTTTTGCTACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGcagcag                                                                       >  2:1113826/1‑100 (MQ=255)
                        gCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCg                                                                     >  1:1150557/1‑100 (MQ=255)
                         ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGTCGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                    <  1:154159/100‑1 (MQ=255)
                         ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                    >  1:24473/1‑100 (MQ=255)
                          cccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaa                                                                   >  1:1213542/1‑100 (MQ=255)
                          cccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaa                                                                   >  1:1404724/1‑100 (MQ=255)
                           ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                  >  1:222497/1‑100 (MQ=255)
                           ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGCATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                  >  1:1504986/1‑100 (MQ=255)
                            cTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaaga                                                                 <  2:1018352/100‑1 (MQ=255)
                              aattttGCCACTGATGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAAc                                                               >  1:1072202/3‑100 (MQ=255)
                                    gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                         <  1:154573/100‑1 (MQ=255)
                                    gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                         >  1:1019780/1‑100 (MQ=255)
                                     ccatgggtggCGTTAACATGAAAGGGACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                        <  1:1477736/95‑1 (MQ=255)
                                         tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                    <  2:964620/100‑1 (MQ=255)
                                         tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                    <  1:1297794/100‑1 (MQ=255)
                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                   >  1:1202507/1‑100 (MQ=255)
                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                   >  1:1202529/1‑100 (MQ=255)
                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                   >  2:156428/1‑100 (MQ=255)
                                           gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                  <  2:1202524/100‑1 (MQ=255)
                                           gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                  <  2:1202502/100‑1 (MQ=255)
                                           gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                  >  1:673317/1‑100 (MQ=255)
                                            tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                 <  2:1019774/100‑1 (MQ=255)
                                            tagCGGTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                 <  2:724702/98‑1 (MQ=255)
                                             ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGTTTTCCCTGGTtca                                                <  1:293757/100‑1 (MQ=255)
                                               cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCCATCCTCGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                              <  2:1127145/100‑1 (MQ=255)
                                               cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                              <  1:389911/100‑1 (MQ=255)
                                                       tGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGGCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAc                                      >  2:1109146/1‑100 (MQ=255)
                                                       ngAAAGGAACCTGAATTTTTTTTTTTTTTATTGACGTATCGCTATCATCAATCCACGCCTGTAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAc                                      >  2:1276056/2‑100 (MQ=255)
                                                        gacacgAACCTGAAGCTTGTCTGGTGGTAAAGAGGTAGCGCTACCATCCATCCTCGCCTTCAGCAGCGAAGAACGGCTGCCCTGGTGCATCAGCTGc                                        <  2:1507899/92‑1 (MQ=255)
                                                        gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATAAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                     >  2:850992/1‑100 (MQ=255)
                                                            ggAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCa                                 >  1:664039/1‑100 (MQ=255)
                                                                ccTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                             >  1:795700/1‑100 (MQ=255)
                                                                 cTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCGATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtct                            >  1:1084889/1‑100 (MQ=255)
                                                                  tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                           <  2:12081/100‑1 (MQ=255)
                                                                    aaTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTa                         <  2:81097/100‑1 (MQ=255)
                                                                    aaTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTa                         <  2:721553/100‑1 (MQ=255)
                                                                     aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                        >  1:1274619/1‑100 (MQ=255)
                                                                     aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                        >  1:1274220/1‑100 (MQ=255)
                                                                      ttttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtt                       <  1:163125/100‑1 (MQ=255)
                                                                       tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                      <  1:156428/100‑1 (MQ=255)
                                                                        ttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtttt                     >  2:381032/1‑100 (MQ=255)
                                                                         tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                    >  1:1190050/1‑100 (MQ=255)
                                                                          ttCTGGTGGTAATTACGTATCGCTAGCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTGTca                   <  2:760962/100‑2 (MQ=255)
                                                                             tggtggTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCtcc                 >  2:1262936/1‑97 (MQ=255)
                                                                              ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg               >  1:64116/1‑100 (MQ=255)
                                                                               gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTACCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg              >  1:715900/1‑100 (MQ=255)
                                                                                 ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCGTCACTTCTTTATTTTCTGCCGGGt            >  1:1098982/1‑100 (MQ=255)
                                                                                 ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt            <  2:1335209/100‑1 (MQ=255)
                                                                                     aTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat        >  1:594794/1‑100 (MQ=255)
                                                                                           tATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgcg  <  1:714391/100‑1 (MQ=255)
                                                                                                 |                                                                                             
ATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCG  >  NC_000913/2191335‑2191525

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: