Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 985,932 G→C 10.4% E232Q (GAA→CAA)  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,9320GC10.4% 142.3 / 7.3 48E232Q (GAA→CAA) pilFtype IV pili biogenesis protein PilF
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base C (2/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 = 4.19e-01

AGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAA  >  NC_002947/985792‑986080
                                                                                                                                            |                                                                                                                                                    
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)
         aCGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAAtt                                                                                                                                                  <  6:105140/136‑1 (MQ=255)
         aCGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAAtt                                                                                                                                                  >  5:105140/1‑136 (MQ=255)
          cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAtt                                                                                                                                      <  3:259899/147‑1 (MQ=255)
          cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAtt                                                                                                                                      <  1:293393/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                                                                                                    <  2:67497/148‑1 (MQ=255)
                                           ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                    <  7:268000/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                                                                                      <  5:106784/141‑1 (MQ=255)
                                                                gAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGa                                                                                      >  6:106784/1‑141 (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                                        >  3:61125/1‑117 (MQ=255)
                                                                                                                                      cTGGCCCAATTCGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                        <  4:61125/117‑1 (MQ=255)
                                                                                                                                         gCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGt        >  8:184192/1‑146 (MQ=255)
                                                                                                                                          ccGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTgag     >  5:316850/1‑148 (MQ=255)
                                                                                                                                            gAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGaaa  >  7:273899/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
AGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAA  >  NC_002947/985792‑986080

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: