Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 152,036 C→T V66I (GTA→ATA)  yadL ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913152,0360CT98.3% 180.0 / ‑6.6 58V66I (GTA→ATA) yadLputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/29);  total (28/30)
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

AGCAGTTCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAC‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAAGAGAGTAC  >  NC_000913/151905‑152162
                                                                                                                                     |                                                                                                                                
aGCAGTTCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCa                                                                                                                         <  2:488746/139‑1 (MQ=255)
     ttCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATc                                                                                                                    >  1:429404/1‑139 (MQ=255)
     ttCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATc                                                                                                                    <  2:1458177/139‑1 (MQ=255)
     ttCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATc                                                                                                                    >  2:1418838/1‑139 (MQ=255)
               tgtgGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          <  1:1418838/139‑1 (MQ=255)
                 tgGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAAt                                                                                                        >  2:1224008/1‑139 (MQ=255)
                      cGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCt                                                                                                   >  1:1084124/1‑139 (MQ=255)
                        cTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATAT‑‑ATCACCAAAATCAACAACAGAAATATCCt                                                                                                   >  2:1413298/1‑139 (MQ=255)
                              tAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCagag                                                                                           <  1:52284/139‑1 (MQ=255)
                               aaTTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATATATATCACCAAAATCAACAACAGAAATATCCTGATCaga                                                                                            <  1:1413298/139‑1 (MQ=255)
                                      tttGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTAc                                                                                   >  2:1347067/1‑139 (MQ=255)
                                            cTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGt                                                                             <  1:95946/139‑1 (MQ=255)
                                            cTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGt                                                                             <  1:677423/139‑1 (MQ=255)
                                                         gATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          >  2:610768/1‑97 (MQ=255)
                                                         gATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          <  1:610768/97‑1 (MQ=255)
                                                            aCCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAccc                                                             <  1:1224008/139‑1 (MQ=255)
                                                                 cacaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          >  1:1318296/1‑89 (MQ=255)
                                                                 cacaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          <  2:1318296/89‑1 (MQ=255)
                                                                 cacaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTc                                                        >  1:1378931/1‑139 (MQ=255)
                                                                  acaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCt                                                       >  2:1119949/1‑139 (MQ=255)
                                                                    aGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                     <  2:1461017/139‑1 (MQ=255)
                                                                    aGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                     <  1:1119949/139‑1 (MQ=255)
                                                                      tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCagag                                                                                           >  2:911077/1‑99 (MQ=255)
                                                                      tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCagag                                                                                           <  1:911077/99‑1 (MQ=255)
                                                                      tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCAttt                                                         >  2:1399921/1‑133 (MQ=255)
                                                                      tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCAttt                                                         <  1:1399921/133‑1 (MQ=255)
                                                                        tttGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          <  1:554864/82‑1 (MQ=255)
                                                                        tttGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          >  2:554864/1‑82 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          <  2:364297/74‑1 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          >  1:364297/1‑74 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCagag                                                                                           >  2:34921/1‑89 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCagag                                                                                           <  1:34921/89‑1 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTa                                                                <  1:1383324/116‑1 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTa                                                                >  2:1383324/1‑116 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAAc                                         >  1:1040805/1‑139 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAAc                                         >  1:1228056/1‑139 (MQ=255)
                                                                                tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAAc                                         >  1:1273923/1‑139 (MQ=255)
                                                                                   gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                      <  1:1347067/139‑1 (MQ=255)
                                                                                   gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                      <  1:485887/139‑1 (MQ=255)
                                                                                   gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                      <  2:429404/139‑1 (MQ=255)
                                                                                        aaTGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCAttt                                                                                     >  2:7265/1‑87 (MQ=255)
                                                                                        aaTGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCAttt                                                                                     <  1:7265/87‑1 (MQ=255)
                                                                                        aaTGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTc                                 <  1:233274/139‑1 (MQ=255)
                                                                                                    ttGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCa                                                                                              >  2:1272899/1‑66 (MQ=255)
                                                                                                    ttGGTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCa                                                                                              <  1:1272899/66‑1 (MQ=255)
                                                                                                      ggTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAccc                                                             >  1:852311/1‑97 (MQ=255)
                                                                                                      ggTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAccc                                                             <  2:852311/97‑1 (MQ=255)
                                                                                                      ggTCTTGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGtt                   >  1:47188/1‑139 (MQ=255)
                                                                                                          ttGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATc                                                                                                                    >  2:675047/1‑38 (MQ=255)
                                                                                                          ttGGCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATc                                                                                                                    <  1:675047/38‑1 (MQ=255)
                                                                                                            ggCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          <  1:11654/46‑1 (MQ=255)
                                                                                                            ggCATTGATTTCAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGaa                                                                                                          >  2:11654/1‑46 (MQ=255)
                                                                                                                   aTTTCAGAA‑‑ATATATAC‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAagaga      <  2:1040805/139‑1 (MQ=255)
                                                                                                                       cAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGc                                                                          >  1:729695/1‑67 (MQ=255)
                                                                                                                       cAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGc                                                                          <  2:729695/67‑1 (MQ=255)
                                                                                                                       cAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                     >  1:1485346/1‑88 (MQ=255)
                                                                                                                       cAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                     <  2:1485346/88‑1 (MQ=255)
                                                                                                                       cAGAA‑‑ATATATAT‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAAGAGAGTAc  >  2:1452373/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                
AGCAGTTCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAA‑‑ATATATAC‑‑ATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAAGAGAGTAC  >  NC_000913/151905‑152162

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