Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,004,188 T→C L81P (CTG→CCG)  ycbF → putative fimbrial chaperone YcbF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,004,1880TC97.7% 143.9 / ‑4.2 44L81P (CTG→CCG) ycbFputative fimbrial chaperone YcbF
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/18);  total (26/18)
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.50e-01

CGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCAAA  >  NC_000913/1004059‑1004324
                                                                                                                                 |                                                                                                                                        
cGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCg                                                                                                                                >  2:81298/1‑140 (MQ=255)
cGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCg                                                                                                                                >  1:597037/1‑140 (MQ=255)
                                aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                                                                                                >  1:166325/1‑140 (MQ=255)
                                aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                                                                                                >  1:416990/1‑140 (MQ=255)
                                  tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                                                                                              >  1:94689/1‑140 (MQ=255)
                                     tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                                                                                           <  1:229870/140‑1 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  2:336574/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGGCAGGCACGTTgcgc                                                                                      >  1:421891/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  1:262725/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  1:808498/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  2:336908/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  1:399329/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  2:728530/1‑140 (MQ=255)
                                          gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                                                      >  2:477073/1‑140 (MQ=255)
                                           tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt                                                                                     >  1:645904/1‑140 (MQ=255)
                                                            ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGc                                                                    <  2:155690/140‑1 (MQ=255)
                                                              aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGa                                                                        >  1:795147/1‑134 (MQ=255)
                                                              aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGa                                                                        <  2:795147/134‑1 (MQ=255)
                                                                   tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACAt                                                               <  1:388276/138‑1 (MQ=255)
                                                                   tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACAt                                                               >  2:388276/1‑138 (MQ=255)
                                                                   tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCt                                                             <  1:355414/140‑1 (MQ=255)
                                                                   tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCt                                                             <  2:808498/140‑1 (MQ=255)
                                                                            ccGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTgg                                                    <  2:658649/140‑1 (MQ=255)
                                                                             cGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGa                                                   >  2:464409/1‑140 (MQ=255)
                                                                             cGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGa                                                   >  2:424285/1‑140 (MQ=255)
                                                                                    cgcgTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGaa                                            >  2:267360/1‑140 (MQ=255)
                                                                                     gcgTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGaaa                                           <  1:336908/140‑1 (MQ=255)
                                                                                      cgTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAAc                                          >  1:330398/1‑140 (MQ=255)
                                                                                      cgTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAAc                                          <  1:81298/140‑1 (MQ=255)
                                                                                          gggCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGaaa                                           >  2:423866/1‑135 (MQ=255)
                                                                                          gggCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGaaa                                           <  1:423866/135‑1 (MQ=255)
                                                                                          gggCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTa                                      >  2:784096/1‑140 (MQ=255)
                                                                                           ggCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTAt                                     <  2:597037/140‑1 (MQ=255)
                                                                                             ccGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTAttt                                   <  2:399329/140‑1 (MQ=255)
                                                                                                 gggggAAGCGTCGACAGTGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGt                               <  1:267360/140‑1 (MQ=255)
                                                                                                                  tGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTg              >  1:380675/1‑140 (MQ=255)
                                                                                                                  tGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTg              <  1:784096/140‑1 (MQ=255)
                                                                                                                  tGCCAGCACCATTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTg              <  2:330398/140‑1 (MQ=255)
                                                                                                                            aTTACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCa    >  2:709363/1‑140 (MQ=255)
                                                                                                                             ttACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCaa   <  1:424285/140‑1 (MQ=255)
                                                                                                                             ttACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCaa   <  1:328480/140‑1 (MQ=255)
                                                                                                                             ttACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCaa   <  2:380675/140‑1 (MQ=255)
                                                                                                                              tACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCaaa  >  1:384354/1‑140 (MQ=255)
                                                                                                                              tACCGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCaaa  >  2:322601/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                        
CGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTTGAGTTAAGCATTGCCAGCGTGCCATCCGGCAAA  >  NC_000913/1004059‑1004324

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