Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 114,992 G→A 13.8% Q245* (CAG→TAG)  hofC ← assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913114,9920GA13.8% 78.2 / 4.9 29Q245* (CAG→TAG) hofCassembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base A (2/2);  total (13/16)
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.57e-01

CAGAGAGCCGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACG  >  NC_000913/114705‑115240
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
cAGAGAGCCGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTTAACCCATAATCGGGATGCGTAACAGCAATTTCTg                                                                                                                                                                                                                                                          >  1:166345/1‑288 (MQ=255)
       ccGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTGGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGACCCCATCATCGTGCTGCGTAACCGTCATTTCTGCCGccata                                                                                                                                                                                                                                                  >  2:175314/1‑285 (MQ=255)
        cGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCg                                                                                                                                                                                                                                                       <  1:20437/283‑1 (MQ=255)
                       ttCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAATAGCAATTTCTACCGCACTATAAGCCAGgtcggtcg                                                                                                                                                                                                                                  >  2:275290/1‑289 (MQ=255)
                            aCTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCg                                                                                                                                                                                                                                                       <  2:55178/263‑1 (MQ=255)
                            aCTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCg                                                                                                                                                                                                                                                       >  1:55178/1‑263 (MQ=255)
                                              gggCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAATAGCAATTTCTACCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCt                                                                                                                                                                                                            <  1:275290/288‑1 (MQ=255)
                                              gggCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCt                                                                                                                                                                                                            <  2:105331/288‑1 (MQ=255)
                                                    aaCTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTACCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCca                                                                                                                                                                                                     >  2:206777/1‑289 (MQ=255)
                                                                  tttttggggcagccggATGGGTTGACCGTTACTGATATCGTGCTGGATTGGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGaac                                                                                                                                                                                       <  2:166345/280‑1 (MQ=255)
                                                                         gccagccaGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGGGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACca                                                                                                                                                                                <  2:38598/289‑1 (MQ=255)
                                                                                                         gTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGc                                                                                                                                                 >  2:117905/1‑288 (MQ=255)
                                                                                                                                      aGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACca                                                                                                                                                                                <  1:300686/228‑1 (MQ=255)
                                                                                                                                      aGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACca                                                                                                                                                                                >  2:300686/1‑228 (MQ=255)
                                                                                                                                           gggCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCAt                                                                                                                                        <  1:68317/263‑1 (MQ=255)
                                                                                                                                           gggCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCAt                                                                                                                                        >  2:68317/1‑263 (MQ=255)
                                                                                                                                                      tttttCTCTGACGCGCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTTTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAg                                                                                                   <  1:175314/287‑1 (MQ=255)
                                                                                                                                                      tGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTACCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAg                                                                                                   <  1:206777/289‑1 (MQ=255)
                                                                                                                                                                   ctctctACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTatagata                                                                                       >  1:130971/1‑288 (MQ=255)
                                                                                                                                                                      tctACGCCCTGTAAAAAAGTAATGCCTGCACTTTTGGTCAGCGCCAAAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCgg                                                                                   <  2:202618/289‑1 (MQ=255)
                                                                                                                                                                          cGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCaaa                                                                               >  1:76201/1‑289 (MQ=255)
                                                                                                                                                                                gTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAg                                                                                                   <  1:316750/263‑1 (MQ=255)
                                                                                                                                                                                gTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAg                                                                                                   >  2:316750/1‑263 (MQ=255)
                                                                                                                                                                                           aTGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTg                                                                                                              >  1:65575/1‑241 (MQ=255)
                                                                                                                                                                                           aTGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTg                                                                                                              <  2:65575/241‑1 (MQ=255)
                                                                                                                                                                                                                        cGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGcc                                 <  1:117905/289‑1 (MQ=255)
                                                                                                                                                                                                                                   gCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGCAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCAACCACCACCATGATGGCCATCACAAAAat                      >  1:26183/1‑289 (MQ=255)
                                                                                                                                                                                                                                                       aTCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGAtaacgtaacg  <  2:76201/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aTAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCt                           <  2:194652/252‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aTAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCt                           >  1:194652/1‑252 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
CAGAGAGCCGGATGCCTCTCCTGTTCTCACTAATTGCAAACAGAGCGGGCTAAACTCACCGGTATTTTTTAGCGCCAGCCAGATGGGTTGACCGTTACTGATATCGTGCTGGATTTGTGTCAGAAGTTGCACCCAGTACGGGCAGCGCATTGTTTCTCTGACGCTCTCTACGCCCTGTAAAAAAGTAATGCCTGCACTTTGTGTCAGCGCCAGAATCGTAAAGATCTGCGTGAGTTTTTGTCCCCGCATCAGTGAACCCATAATCGGGATGCGTAACAGCAATTTCTGCCGCACTATAAGCCAGGTCGGTCGGCGCATCAGCAACTTATTGGCTATCGCCAGCAGAAAGCCGAACAACACCAGCAGCCAGCTCCATTCGCCACTAAAGTCTGCCAGCGTCATGATCCCCTGCGTTAGTGCCGGTAGTGGGGTGTTGAAGGTCTTATAGATAGCGGCAAACTCCGGCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACG  >  NC_000913/114705‑115240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: