Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,543,600 A→G 6.4% A162A (GCA→GCG fimA → major type 1 subunit fimbrin (pilin)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,543,6000AG 6.4% 186.7 / 2.9 62A162A (GCA→GCGfimAmajor type 1 subunit fimbrin (pilin)
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base G (2/2);  total (34/28)
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.01e-01

GCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTA  >  NC_000913/4543457‑4543745
                                                                                                                                               |                                                                                                                                                 
gCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCgg                                                                                                                                              <  8:25268/149‑1 (MQ=255)
    cGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGTCTGAACGGGTGGTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCTGGCGCGTTATTTTGCAACCgtgg                                                                                                                                            <  8:95117/147‑1 (MQ=255)
      ggTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAAc                                                                                                                                                 >  6:30637/1‑140 (MQ=255)
      ggTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAAc                                                                                                                                                 <  5:30637/140‑1 (MQ=255)
              aacaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccat                                                                                                                                 <  7:49242‑M2/148‑6 (MQ=255)
               acaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                                                                               >  8:42‑M2/1‑142 (MQ=255)
                caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                                                                              <  8:353000‑M2/149‑9 (MQ=255)
                caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                                                                              >  3:101357‑M2/1‑141 (MQ=255)
                caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                                                                              <  8:83955‑M2/149‑9 (MQ=255)
                     gTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACCGCCCTGAATAACGGAACCACTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCACcccataagcgct                                                                                                                          <  8:151335‑M2/148‑13 (MQ=255)
                            tgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccataagcgctaacttaag                                                                                                                  <  5:12693‑M2/149‑21 (MQ=255)
                               aGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                                                                               <  2:130642‑M2/149‑24 (MQ=255)
                                    cTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaa                                                                                                           >  8:259667‑M2/1‑121 (MQ=255)
                                     tGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCCTTATTTTGCAACCGGGGCCCCAAcccataatcgctaacttaagggatgaac                                                                                                          <  2:197016‑M2/148‑29 (MQ=255)
                                       gACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaacca                                                                                                        >  8:32897‑M2/1‑118 (MQ=255)
                                        aCAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatc                                                                                                      >  1:303033‑M2/1‑117 (MQ=255)
                                             aCGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAATCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagggctaacttatgggttgaaccatctgaa                                                                                                  >  1:197016‑M2/1‑112 (MQ=255)
                                                ggTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCATAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaat                                                                                              >  4:9822‑M2/1‑109 (MQ=255)
                                                 gTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatg                                                                                             >  2:13707‑M2/1‑108 (MQ=255)
                                                 gTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccacctgaagaatg                                                                                             >  2:233330‑M2/1‑108 (MQ=255)
                                                    ctgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcga                                                                                          <  8:88276‑M2/149‑45 (MQ=255)
                                                     tgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgac                                                                                         <  7:212683‑M2/149‑46 (MQ=255)
                                                       cgcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgc                                                                                       >  4:107486‑M2/1‑102 (MQ=255)
                                                         cTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcct                                                                                     >  1:378656‑M2/1‑100 (MQ=255)
                                                              gCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCGTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtg                                                                                <  7:318445‑M2/149‑55 (MQ=255)
                                                                  gATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcct                                                                             >  3:86546‑M2/1‑91 (MQ=255)
                                                                    tGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgt                                                                          >  7:31864‑M2/1‑89 (MQ=255)
                                                                       tGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtt                                                                         <  7:46037‑M2/147‑62 (MQ=255)
                                                                          gACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaacaatctgaagaatgcgacgcctcggtgcctcgttaagac                                                                    <  7:153604‑M2/149‑67 (MQ=255)
                                                                             aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgat                                                                 <  1:209270‑M2/149‑70 (MQ=255)
                                                                              tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatg                                                                <  6:361133‑M2/149‑71 (MQ=255)
                                                                              tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatg                                                                <  4:6416‑M2/149‑71 (MQ=255)
                                                                                   ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcg                                                           <  8:23289‑M2/149‑76 (MQ=255)
                                                                                        aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGTCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttc                                                      >  1:284359‑M2/1‑69 (MQ=255)
                                                                                        aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttc                                                      <  6:199447‑M2/149‑81 (MQ=255)
                                                                                        aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatc                                                                                                      <  4:174112‑M2/101‑33 (MQ=255)
                                                                                            aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttc                                                   <  2:145621‑M2/148‑84 (MQ=255)
                                                                                                   aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcg                                            >  7:83955‑M2/1‑58 (MQ=255)
                                                                                                    aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtt                                          <  5:339140‑M2/149‑93 (MQ=255)
                                                                                                         ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgta                                     >  6:12693‑M2/1‑52 (MQ=255)
                                                                                                         ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgta                                     >  8:46037‑M2/1‑52 (MQ=255)
                                                                                                         ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttg                                       <  7:205748‑M2/147‑96 (MQ=255)
                                                                                                         ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttg                                       <  2:280947‑M2/147‑96 (MQ=255)
                                                                                                          gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtag                                    >  8:49242‑M2/1‑51 (MQ=255)
                                                                                                               aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgt                               >  7:23289‑M2/1‑46 (MQ=255)
                                                                                                                   ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcag                            >  6:339140‑M2/1‑42 (MQ=255)
                                                                                                                    cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaat                                                >  4:200903‑M2/1‑41 (MQ=255)
                                                                                                                    cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaat                                                <  3:200903‑M2/127‑87 (MQ=255)
                                                                                                                      ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtttt                                        <  6:56547‑M2/133‑95 (MQ=255)
                                                                                                                      ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtttt                                        >  5:56547‑M2/1‑39 (MQ=255)
                                                                                                                        ccTTTCCAGTCGCGTTATTTTGCAACCGGGGCCGCAAcccataagggctatgttaagggtttaaccatctgaagaatgagacggctgggtgggtcgttaagacgatgcatcgagttcttcaattgcgttttgtaggctgtcagggata                       >  8:318445‑M2/1‑37 (MQ=255)
                                                                                                                        ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggata                       >  2:209270‑M2/1‑37 (MQ=255)
                                                                                                                             ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtc                  >  7:353000‑M2/1‑32 (MQ=255)
                                                                                                                                ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccac              >  7:88276‑M2/1‑29 (MQ=255)
                                                                                                                                 gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtttaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcatttgcgttttgtaggctgtcagggatactgtcccac              >  4:147192‑M2/1‑28 (MQ=255)
                                                                                                                                  cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacg             <  5:57090‑M2/148‑122 (MQ=255)
                                                                                                                                    cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgcagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaat          >  6:180784‑M2/1‑25 (MQ=255)
                                                                                                                                    cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttctccaattgcgttttgtaggctgtcagggatactgtcccacgaat          >  3:6416‑M2/1‑25 (MQ=255)
                                                                                                                                     gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatg         >  1:43950‑M2/1‑24 (MQ=255)
                                                                                                                                      ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatgg        >  1:145621‑M2/1‑23 (MQ=255)
                                                                                                                                        aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggc       >  8:17943‑M2/1‑21 (MQ=255)
                                                                                                                                             gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacct  >  1:130642‑M2/1‑16 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTA  >  NC_000913/4543457‑4543745

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