Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 985,937 G→C 10.2% L233F (TTG→TTC)  pilF → type IV pili biogenesis protein PilF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947985,9370GC10.2% 125.6 / 7.3 49L233F (TTG→TTC) pilFtype IV pili biogenesis protein PilF
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base C (2/3);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

AGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACGG  >  NC_002947/985792‑986083
                                                                                                                                                 |                                                                                                                                                  
aGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTggg                                                                                                                                                  <  4:162456/148‑1 (MQ=255)
aGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTgg                                                                                                                                                   <  7:181901/147‑1 (MQ=255)
 gTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCACGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGc                                                                                                                                                 <  2:101149/148‑1 (MQ=255)
  ttGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCa                                                                                                                                               <  3:167455/149‑1 (MQ=255)
  ttGTCCTACGAAAACAGGCATTATGTGCCGGCCCAGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCCAATTCGGCCa                                                                                                                                               <  1:225721/149‑1 (MQ=255)
       cTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCaa                                                                                                                                           <  3:42680/148‑1 (MQ=255)
       cTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCCAATTCGGCCAgtgt                                                                                                                                           <  8:87605/148‑4 (MQ=255)
          cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAtt                                                                                                                                         <  1:293393/147‑1 (MQ=255)
          cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAtt                                                                                                                                         <  3:259899/147‑1 (MQ=255)
               aCAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGa                                                                                                                                  <  6:186874/149‑1 (MQ=255)
                 aGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACt                                                                                                                                >  8:268000/1‑149 (MQ=255)
                    cATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                          >  1:24301/1‑136 (MQ=255)
                    cATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                          <  2:24301/136‑1 (MQ=255)
                      ttATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATc                                                                                                                           >  8:42887/1‑149 (MQ=255)
                      ttATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAAATGGGCCAGCAATTACAACGACTTTATc                                                                                                                           <  4:330483/149‑1 (MQ=255)
                                     ggACTACTACGATCGTTTCAGCCAGTTGAGTGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGaa                                                                                                             >  2:49280/1‑148 (MQ=255)
                                        ctactaCGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATc                                                                                                         <  7:42887/149‑1 (MQ=255)
                                           ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                       >  3:231633/1‑148 (MQ=255)
                                           ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                       <  7:268000/148‑1 (MQ=255)
                                           ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                       <  2:67497/148‑1 (MQ=255)
                                           ctacCATCGTTTCAGCCAGTGGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGATTTGGGCCAGCAATTACATCGACTTTATCCCGGTACGCCGGAATATc                                                                                                         <  8:283944/146‑1 (MQ=255)
                                            tacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTCCTGATGGGCAGCCGCCTTGCCAGGCTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAAAATCAGc                                                                                                      >  5:22355/1‑148 (MQ=255)
                                             acGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCaa                                                                                                    <  7:122075/149‑1 (MQ=255)
                                                                gAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGa                                                                                         >  6:106784/1‑141 (MQ=255)
                                                                gAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGa                                                                                         <  5:106784/141‑1 (MQ=255)
                                                                    gACCACGACGCCCGTAGCCTGCTTCTGGGCAGCCTCCTTTCCTGGGTGTTCGACGAGCATGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCTGTACGCCGGAATATGAGCAATACCTGTCGGAGCAACGATGAAc                                                                             <  3:190447/149‑1 (MQ=255)
                                                                    gACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAc                                                                             >  4:153360/1‑149 (MQ=255)
                                                                     aCCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTTCAACGACTTTATCCCGGTACGCCGGAGTATCAGCCATACCTGTCGGAGCAACGATGAAcg                                                                            >  8:87841/1‑149 (MQ=255)
                                                                      ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAcgc                                                                           >  3:254678/1‑149 (MQ=255)
                                                                       cACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAcgcc                                                                          >  5:220453/1‑149 (MQ=255)
                                                                         cgacgCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAcgccgc                                                                        <  4:254678/149‑1 (MQ=255)
                                                                              cccGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGTATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATc                                                                   <  7:226290/149‑1 (MQ=255)
                                                                               ccGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATcc                                                                  >  5:153917/1‑149 (MQ=255)
                                                                                    gCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCTGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGGAACGATGAACGCCGCGCATCCCGaa                                                              <  6:22355/148‑1 (MQ=255)
                                                                                       tgctgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtag                                                          >  6:8763/1‑149 (MQ=255)
                                                                                            ctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagcag                                                       <  7:62374/147‑1 (MQ=255)
                                                                                                              aGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAAcccc                                   >  4:324536/1‑149 (MQ=255)
                                                                                                                 gtgtTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGt                                <  3:153360/149‑1 (MQ=255)
                                                                                                                       gacgaGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCcgcg                                                                       >  1:105338/1‑104 (MQ=255)
                                                                                                                       gacgaGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCcgcg                                                                       <  2:105338/104‑1 (MQ=255)
                                                                                                                             cAGGGCACACTGGCCGAATTGGGCCAGCAATGACAACGACTTTATCCCGGTACGCCGGAATATCAGCACTACCTGTCGGCGCAACGCTGAACGCCGCGCATCCCGAAGTCGCAGTAGCGCCTGGCCCGCACCCCGGTGAGCGTTTgcgc                    >  3:172497/1‑149 (MQ=255)
                                                                                                                                    tgcTGGCCCAATTCGGCCAGCAATTACAACGACTTTATCACGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGccc             >  4:129958/3‑149 (MQ=255)
                                                                                                                                      cTGGCCCAATTCGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           <  4:61125/117‑1 (MQ=255)
                                                                                                                                      cTGGCCCAATTCGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           >  3:61125/1‑117 (MQ=255)
                                                                                                                                         gCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGt           >  8:184192/1‑146 (MQ=255)
                                                                                                                                          ccGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTgag        >  5:316850/1‑148 (MQ=255)
                                                                                                                                            gAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGaaa     >  7:273899/1‑149 (MQ=255)
                                                                                                                                             aaTTGGGCCTGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAAc    >  2:256229/1‑149 (MQ=255)
                                                                                                                                             aaTTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAAc    >  4:51299/1‑149 (MQ=255)
                                                                                                                                             aaTTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAAc    <  6:316850/149‑1 (MQ=255)
                                                                                                                                               ttGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACgg  <  4:231633/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACGG  >  NC_002947/985792‑986083

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: ‑