Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,802,325 A→T 12.5% L170Q (CTG→CAG)  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,3250AT12.5% 119.0 / 8.7 48L170Q (CTG→CAG) pilQtype IV pili biogenesis protein
Reads supporting (aligned to +/- strand):  ref base A (18/24);  new base T (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 = 5.51e-01

CTCCTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCTTGGC  >  NC_002947/5802186‑5802467
                                                                                                                                           |                                                                                                                                              
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                                                                                                                                    <  2:353952/100‑4 (MQ=255)
                                                    cccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCAGCTGGTGGCGCagtt                                                                                                                                    >  1:353952/1‑97 (MQ=255)
                                                     ccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCg                                                                                                                                  >  1:326225/1‑101 (MQ=255)
                                                     ccAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCg                                                                                                                                  <  2:326225/101‑1 (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          <  2:201330/138‑1 (MQ=255)
                                                                                                                                        accaGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATcggcg          >  1:201330/1‑138 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CTCCTTGCCGGGCCGCAGTTTCTCGCCGTACAGCGGCCCACCCCAGCGTACCCCCAGGCTTTTCTCATAATCGACGTTGGCCTCGACGATACGCGCCTCGATCGCTACCTGGCGCACCGGCACATCCAACTGCGCCACCAGCTGGCGCATCTCGGCCAGGCGGTCAGCGGGCAGATGCACCACCAACGTGTTGGTGCGCTCATCGACACTCAAGCTAGCCGGAGGGATGCCGTCACCCTCCAGGCCTGCCAGCAGCAACTGCGCCAGATCGGCGGCCTTGGC  >  NC_002947/5802186‑5802467

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