Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,545,735 G→T 9.5% G214W (GGG→TGG)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,545,7350GT 9.5% 150.9 / 3.6 48G214W (GGG→TGG) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base T (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 = 1.00e+00

GCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTC  >  NC_000913/4545598‑4545871
                                                                                                                                         |                                                                                                                                        
gcgcgcgTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTggg                                                                                                                                        <  1:166994/140‑1 (MQ=255)
  gcgcgTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGtt                                                                                                                                      <  2:522381/140‑1 (MQ=255)
  gcgcgTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGtt                                                                                                                                      <  2:840416/140‑1 (MQ=255)
  gcgcgTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGtt                                                                                                                                      <  1:522601/140‑1 (MQ=255)
   cgcgTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTa                                                                                                                                     >  1:1109059/1‑140 (MQ=255)
   cgcgTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTa                                                                                                                                     >  1:26066/1‑140 (MQ=255)
   cgcgTGGATATATTCCACCTGAGTTAAGGGATCCCGTTAATAATGCCGGATTGCTCAATTAAAATTTCATCGGAAATAGTGTACAGAATCGGATAGGGTGAAACATCCATTAAGCATATTAAAACCTACAGATTGGGATa                                                                                                                                     >  1:1066984/1‑140 (MQ=255)
                    cctGAGTTATGTGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATCTTGGTGCGTGTCg                                                                                                                    >  2:369917/1‑140 (MQ=255)
                    cctGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCg                                                                                                                    >  1:476885/1‑140 (MQ=255)
                    cctGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCg                                                                                                                    >  1:269657/1‑140 (MQ=255)
                            aTGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTAcgcg                                                                                                            >  1:956758/1‑140 (MQ=255)
                                  tCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATa                                                                                                      >  1:860241/1‑140 (MQ=255)
                                    ccGGTATTAATGCCGTATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTTGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATacc                                                                                                    >  2:159710/1‑140 (MQ=255)
                                    ccGGTATTAATGCCGTATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTTGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATacc                                                                                                    >  2:159513/1‑140 (MQ=255)
                                       gTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATaccacc                                                                                                 >  2:1001772/1‑140 (MQ=255)
                                        tATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTTGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCt                                                                                                <  2:606575/140‑1 (MQ=255)
                                        tATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCt                                                                                                <  2:674423/140‑1 (MQ=255)
                                        tATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCt                                                                                                <  2:435365/140‑1 (MQ=255)
                                        tATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCt                                                                                                <  2:608032/140‑1 (MQ=255)
                                        tATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCt                                                                                                <  2:217421/140‑1 (MQ=255)
                                             aTGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGt                                                                                           >  2:526089/1‑140 (MQ=255)
                                             aTGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGt                                                                                           >  2:142620/1‑140 (MQ=255)
                                                ccGGATTGCACAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTtat                                                                                        <  2:39648/140‑1 (MQ=255)
                                                                       aGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACAACCTGGAGTTTTAACAGTAGCGACAGATCATCAgg                                                                 >  2:226041/1‑140 (MQ=255)
                                                                       aGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGGGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAgg                                                                 >  2:225021/1‑140 (MQ=255)
                                                                                       aGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATgg                                                 <  1:1061007/140‑1 (MQ=255)
                                                                                       aGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATgg                                                 <  2:300461/140‑1 (MQ=255)
                                                                                       aGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATgg                                                 <  1:1060749/140‑1 (MQ=255)
                                                                                       aGAATCGGATTGGGGGTAACAGCAATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATgg                                                 <  1:1061393/140‑1 (MQ=255)
                                                                                          aTCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcag                                              <  2:574368/140‑1 (MQ=255)
                                                                                          aTCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcag                                              <  2:366139/140‑1 (MQ=255)
                                                                                           tCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcagc                                             <  2:818198/140‑1 (MQ=255)
                                                                                           tCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcagc                                             >  1:211486/1‑140 (MQ=255)
                                                                                           tCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcagc                                             <  1:494725/140‑1 (MQ=255)
                                                                                            cGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcagca                                            >  1:248844/1‑140 (MQ=255)
                                                                                                   ggggTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAAt                                     >  2:552021/1‑140 (MQ=255)
                                                                                                      gTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATAcc                                  <  1:527326/140‑1 (MQ=255)
                                                                                                               aTTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAg                         >  1:619877/1‑140 (MQ=255)
                                                                                                               aTTATGCATATTGAAACCTCCAGAGTGCGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCCCCTGGAGTTATAACAGTAGCGACTGATCATCAGGTAGCAACACTAAATGGCAGCATTTCAATACCGGGCTTGAg                         <  2:135853/140‑1 (MQ=255)
                                                                                                               aTTATGCATAATTAAAGCAACAGACTGGGTTAAATTTTGGTTCGTGGCGTTTAAGCTAAAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAATTATAAATGGCTGCATATCAATTCCTGGCTTGat                         <  2:787806/140‑2 (MQ=255)
                                                                                                                 tATGCATATTTAAACCTACAGAGTTGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACTGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCg                       <  1:159710/140‑1 (MQ=255)
                                                                                                                 tATGCATATTTAAACCTACAGAGTTGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCg                       <  1:159513/140‑1 (MQ=255)
                                                                                                                 tATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCg                       <  2:487989/140‑1 (MQ=255)
                                                                                                                 tATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCg                       <  2:487924/140‑1 (MQ=255)
                                                                                                                 tATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCg                       <  1:282001/140‑1 (MQ=255)
                                                                                                                     cATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGAc                   >  1:998539/1‑140 (MQ=255)
                                                                                                                           tAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACataata             <  2:998539/140‑1 (MQ=255)
                                                                                                                                  aCAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTAc      >  1:999470/1‑140 (MQ=255)
                                                                                                                                      agTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAAGGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTc  <  1:737128/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTC  >  NC_000913/4545598‑4545871

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