Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,802,323 C→G 12.5% V171L (GTG→CTG)  pilQ ← type IV pili biogenesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,802,3230CG12.5% 136.0 / 10.2 48V171L (GTG→CTG) pilQtype IV pili biogenesis protein
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base G (3/3);  total (21/27)
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.76e-01

TCGACGAACAGCTCCTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCTTGGC  >  NC_002947/5802175‑5802467
                                                                                                                                                    |                                                                                                                                                
tCGACGAACAGCTCCTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGccac                                                                                                                                                  >  4:316563/1‑149 (MQ=255)
           cTCCTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCa                                                                                                                                       >  2:264977/1‑149 (MQ=255)
             ccTTGCCGGGCAGCAGTTTCTCGCCGTACAGCGGCACACCACAGCGTACCACCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATAGCTACCTTGCGCACCGGCAAATCCAACTGCGCCACCAGCTGGAGCAta                                                                                                                                     >  1:259042/1‑148 (MQ=255)
              cTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCAGCTGGTGGCGCagtt                                                                                                                                    <  7:215253/149‑4 (MQ=255)
                   cGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCAGCTGGTGGCGCagttg                                                                                                                                   >  2:182716/1‑141 (MQ=255)
                   cGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCAGCTGGTGGCGCagttg                                                                                                                                   <  1:182716/145‑5 (MQ=255)
                     ggCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAg                                                                                                                             <  4:346897/149‑1 (MQ=255)
                               tctcGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGc                                                                                                                    <  3:316563/148‑1 (MQ=255)
                                   tccGTACAGCGGCCCACCCCAGCGTACCCACAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTTCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCggg                                                                                                                 <  5:345808/146‑1 (MQ=255)
                                   gCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCggg                                                                                                                 <  3:118113/147‑1 (MQ=255)
                                   gCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGc                                                                                                                <  8:229947/148‑1 (MQ=255)
                                   gCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCa                                                                                                               <  4:79327/149‑1 (MQ=255)
                                   gCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCa                                                                                                               <  1:203554/149‑1 (MQ=255)
                                         cAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGca                                                                                                         >  7:354494/1‑149 (MQ=255)
                                                 cACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAAc                                                                                                 <  7:145772/149‑1 (MQ=255)
                                                 cACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAAc                                                                                                 <  1:264977/149‑1 (MQ=255)
                                                    cccAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACgtg                                                                                              >  4:359498/1‑149 (MQ=255)
                                                       aGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTg                                                                                           <  8:354494/149‑1 (MQ=255)
                                                              ccccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCAtctc                                                                                                                                   >  1:162910/1‑102 (MQ=255)
                                                              ccccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCAtctc                                                                                                                                   <  2:162910/102‑1 (MQ=255)
                                                               cccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCAGCTGGTGGCGCagtt                                                                                                                                    >  1:353952/1‑97 (MQ=255)
                                                               cccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCAGCTGGTGGCGCagtt                                                                                                                                    <  2:353952/100‑4 (MQ=255)
                                                                ccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCg                                                                                                                                  <  2:326225/101‑1 (MQ=255)
                                                                ccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCg                                                                                                                                  >  1:326225/1‑101 (MQ=255)
                                                                 cAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAATGTGTTGGTGCGCTCAt                                                                                 <  3:83512/149‑1 (MQ=255)
                                                                  aGGCTTTTCTCATAATCGCCGGTGGCCTCGACGATACTCGCCGCGTTCTCTACCTTGCTCACCGTCACATCCAACTGCGCCACCAGCTGGCGCATCTCGTCCAGGCGGTCAGAGGGCAGATGCACCAACAACGTGGTGTTGCGCTCATc                                                                                <  2:259042/149‑1 (MQ=255)
                                                                        ttCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGaca                                                                            <  3:359498/147‑1 (MQ=255)
                                                                          ctcATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTc                                                                         >  2:142088/1‑148 (MQ=255)
                                                                             aTAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAAGTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTTGTGCGCTCATCGAAACTCa                                                                        <  4:243476/146‑1 (MQ=255)
                                                                               aaTCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGc                                                                     >  5:26063/1‑147 (MQ=255)
                                                                                 tCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGTGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCaa                                                                       <  5:190172/143‑1 (MQ=255)
                                                                                 tCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGc                                                                 >  7:226784/1‑149 (MQ=255)
                                                                                 tCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACAATCaa                                                                       >  6:190172/1‑143 (MQ=255)
                                                                                           ccTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATg                                                       >  1:275874/1‑149 (MQ=255)
                                                                                            cTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTa                                                                   <  2:208669/136‑1 (MQ=255)
                                                                                             tCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAggg                                                          >  5:106765/1‑144 (MQ=255)
                                                                                             tCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAggg                                                          <  6:106765/144‑1 (MQ=255)
                                                                                                   aTACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGTATGCCGTCAcc                                               <  1:142088/149‑1 (MQ=255)
                                                                                                         gcCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAg                                         >  6:102747/1‑149 (MQ=255)
                                                                                                           cTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGc                                       >  1:46687/1‑149 (MQ=255)
                                                                                                               aTCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGc                                   >  3:83502/1‑149 (MQ=255)
                                                                                                                        tGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCcagca                              <  2:179774/145‑1 (MQ=255)
                                                                                                                        tGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCcagca                              >  1:179774/1‑145 (MQ=255)
                                                                                                                            gcACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTgc                      <  2:46687/149‑1 (MQ=255)
                                                                                                                             cACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTgcg                     >  7:141947/1‑149 (MQ=255)
                                                                                                                                          aaCTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCggcggc        >  3:69469/1‑149 (MQ=255)
                                                                                                                                             tGCGCCAGCTGGTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCtt     >  5:15968/1‑149 (MQ=255)
                                                                                                                                                gcCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCTTGGc  <  3:157319/149‑1 (MQ=255)
                                                                                                                                                gcCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCTTGGc  <  8:141947/149‑1 (MQ=255)
                                                                                                                                                   accaGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATcggcg          >  1:201330/1‑138 (MQ=255)
                                                                                                                                                   accaGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATcggcg          <  2:201330/138‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCGACGAACAGCTCCTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCTTGGC  >  NC_002947/5802175‑5802467

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