Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,692,749 G→A 11.1% G125S (GGT→AGT)  PP_2360 → type I pili subunit CsuA/B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,692,7490GA11.1% 94.4 / 3.3 36G125S (GGT→AGT) PP_2360type I pili subunit CsuA/B
Reads supporting (aligned to +/- strand):  ref base G (17/15);  new base A (2/2);  total (19/17)
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.95e-01

GTGCACGCCTGGCGTGGCCATGAGCATGAGCCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACGGCCAGCGCCACCCAT  >  NC_002947/2692601‑2692897
                                                                                                                                                    |                                                                                                                                                    
gTGCACGCCTGGCGTGGCCATGAGGATGAGCCTGGACGGCGGGCAAAACTGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCGGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCg                                                                                                                                                      <  7:208882/149‑1 (MQ=255)
      gCCTGGCGTGGCCATGAGCATGAGCCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCaga                                                                                                                                                <  4:31278/149‑1 (MQ=255)
                   atgagcatgagcCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAg                                                                                                                                   <  3:199421/149‑1 (MQ=255)
                   atgagcatgagcCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAg                                                                                                                                   <  8:151471/149‑1 (MQ=255)
                   atgagcatgagcCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAg                                                                                                                                   >  5:240112/1‑149 (MQ=255)
                        catgagcCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTAtt                                                                                                                              >  7:14697/1‑149 (MQ=255)
                            agcCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTccga                                                                                                                          <  6:240112/149‑1 (MQ=255)
                                         ggCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATc                                                                                                             <  8:196758/149‑1 (MQ=255)
                                            aaaaCAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAg                                                                                                          >  8:253107/1‑149 (MQ=255)
                                               cgtgTGCCACTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTa                                                                                                       >  2:31333/4‑149 (MQ=255)
                                                     ccAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTccga                                                                                                                          <  3:60606/124‑1 (MQ=255)
                                                     ccAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTccga                                                                                                                          >  4:60606/1‑124 (MQ=255)
                                                      cAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCAGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGgtg                                                                                                 >  5:320891/1‑148 (MQ=255)
                                                      cAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCAGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGgtg                                                                                                 >  7:95785/1‑148 (MQ=255)
                                                              cacGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTAcgggcgg                                                                                        <  4:128758/149‑1 (MQ=255)
                                                               acGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACg                                                                                                                   <  1:60871/121‑1 (MQ=255)
                                                               acGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACg                                                                                                                   >  2:60871/1‑121 (MQ=255)
                                                                 gCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGTTCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACg                                                                                                                   <  5:41220/119‑1 (MQ=255)
                                                                 gCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGTTCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACg                                                                                                                   >  6:41220/1‑119 (MQ=255)
                                                                        gAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTATCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGa                                                                              >  6:124650/1‑149 (MQ=255)
                                                                                        aCGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAg                                                                                                          >  8:196362/1‑105 (MQ=255)
                                                                                        aCGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAg                                                                                                          <  7:196362/105‑1 (MQ=255)
                                                                                          gCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGCACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAg                                                            <  7:144562/149‑1 (MQ=255)
                                                                                                 cTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCAGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTa                                                     <  8:95785/149‑1 (MQ=255)
                                                                                                 cTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCAGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTa                                                     <  6:320891/149‑1 (MQ=255)
                                                                                                             ttGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATgtggt                                         >  4:256273/1‑149 (MQ=255)
                                                                                                                tACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGc                                       <  7:253107/148‑1 (MQ=255)
                                                                                                                      gACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGAGGTGCAAGTGAc                                >  2:223644/1‑149 (MQ=255)
                                                                                                                         gccgccTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGt                             >  2:219470/1‑149 (MQ=255)
                                                                                                                          ccgccTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGc                                       >  4:108188/1‑138 (MQ=255)
                                                                                                                          ccgccTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGc                                       <  3:108188/138‑1 (MQ=255)
                                                                                                                               tataGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTg                       >  5:125005/1‑149 (MQ=255)
                                                                                                                                    ccAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAAc                  >  8:42448/1‑149 (MQ=255)
                                                                                                                                     cAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACg                 <  6:125005/149‑1 (MQ=255)
                                                                                                                                       ggCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACGGc               <  8:14697/149‑1 (MQ=255)
                                                                                                                                         ccTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACGGc               >  8:88723/1‑147 (MQ=255)
                                                                                                                                                   cGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACGGCCAGCGCCACCCa   >  5:87434/1‑149 (MQ=255)
                                                                                                                                                    ggTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACGGCCAGCGCCACCCAt  <  3:256273/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGCACGCCTGGCGTGGCCATGAGCATGAGCCTGGACGGCGGGCAAAACAGTGCCAGTGGCACACGCAACCTGAAGCGCACCAGTGGCACGCAAGTGCTGGCCTACCAGTTGTACCGGGACGCCGCCTATAGCCAGGCCTTGGGGATCGGTCAGAGCGTGGCGGTGAGCTATTCCGACCCGACGGCGATCAAGCTACCGGTGTACGGGCGGGTTCAGCTGACCGGCACGCAGCCGGCAGGGACCTACACCGATGTGGTGCAAGTGACGGTGACCTGGTAACGGCCAGCGCCACCCAT  >  NC_002947/2692601‑2692897

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: