Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 156,068 A→G V45A (GTA→GCA)  yadV ← putative fimbrial chaperone YadV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913156,0680AG97.1% 98.1 / ‑4.2 34V45A (GTA→GCA) yadVputative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (13/20);  total (13/21)
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.44e-01

TTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGT  >  NC_000913/155824‑156310
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ttCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTGGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTg                                                                                                                                                                                                                                               <  2:22953/250‑1 (MQ=255)
            ttttggatgAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTACTTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGAGGCGTAGCAGTAAACGGGACTGTAATACTGCCAGGCTCCGCGTTGTCATCGCCAGTATCTAACCAACTCTGGCCCAGCAACGGGTTATTCCCTTTATTTTCCATACGTGCGTTGACCCTTTTTTgtt                                                                                                                                                                                                                                  >  2:332642/9‑249 (MQ=255)
                     aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTtat                                                                                                                                                                                                                         <  2:140175/251‑1 (MQ=255)
                      aCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTtata                                                                                                                                                                                                                        <  2:364668/251‑1 (MQ=255)
                        ttCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTtatatt                                                                                                                                                                                                                      >  1:313401/1‑250 (MQ=255)
                            aGTACGTTAAATCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTGGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTAc                                                                                                                                                                                                                  <  1:278867/251‑1 (MQ=255)
                                    aaaCCAGAACACGCTCGCTCTGTCGTTAGTCAGTGAAGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACTTTTGGCATCAATTCGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACGGCCATGCTCAGCGTTGTCATCGCCAGTATCTGACCAACTCTTGACAATCATCGGGTTATTCCCATTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTa                                                                                                                                                                                                           <  2:112039/250‑1 (MQ=255)
                                              acGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGAc                                                                                                                                                                                                >  2:244537/1‑251 (MQ=255)
                                                                      gAGGTGCTTGCTGTGTACATTAATTTCATTTTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGAATGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGGATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTATTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGTCGATGAAc                                                                                                                                                                         <  2:322385/250‑1 (MQ=255)
                                                                                        aTTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTAc                                                                                                                                                      <  2:265321/251‑1 (MQ=255)
                                                                                         ttAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACtt                                                                                                                                                                                                                                         >  2:91073/1‑167 (MQ=255)
                                                                                         ttAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACtt                                                                                                                                                                                                                                         <  1:91073/167‑1 (MQ=255)
                                                                                                        gCCCACGTTTGGCATCAATACGCGATGCTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCGAACCAACTCTGGACAAGCAACGGGTTATTCCCTTGATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgt                                                                                                                                       <  2:335671/250‑1 (MQ=255)
                                                                                                             cGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGCCTGTACTACTGCCAGGATAAGCGCTGTCATCGCAAGTATCTAACCCACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCCGACGTGCGTTGACACTTTTTTGATCGATTTTATATATTACGCGAGTACCCGAAATGCCAATGTCCGCCATAGACGATGAACTAAAAGCCATACAGGTTAAGAAGCATAAAGCTGTTGTATg                                                                                                                                 >  1:481694/1‑251 (MQ=255)
                                                                                                                                          cgTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCt                                                                                                    >  1:142254/1‑251 (MQ=255)
                                                                                                                                                 gTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATag                                                                                               >  1:354235/1‑249 (MQ=255)
                                                                                                                                                   aaaaGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGGTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGtt                                                                                           <  2:354235/251‑1 (MQ=255)
                                                                                                                                                   aaaaGGGACTGTAATACTGCCAGGCTCAGCGTTGTCACCGCCAGTATCTACCCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAACACATAAACCATCCCTGATAAAGtt                                                                                           >  2:69288/1‑251 (MQ=255)
                                                                                                                                                             gTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACa                                                                                 <  1:152324/251‑1 (MQ=255)
                                                                                                                                                                atcctgccatgcgctgCGTTGTCATCGCTAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATTGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTTTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAggg                                                                              <  2:396026/236‑1 (MQ=255)
                                                                                                                                                                aTACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAGTTAAAACAggg                                                                              <  1:417447/251‑1 (MQ=255)
                                                                                                                                                                aTACCGCCAGGCTCAGCGTTGCCATCGCCAGAATCTAACCAACTCTGGACAAGCAACTGGTTATTCCCGTTCTTTTCCAGACGTGCGTTGATACTTTCTTGATCGCTGTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGTCGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGGATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAggt                                                                              <  2:116668/251‑2 (MQ=255)
                                                                                                                                                                  ccTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaa                                                                            <  2:142254/250‑1 (MQ=255)
                                                                                                                                                                    cgCCAGGCTCAGCGTTGTCATCGCCAGTAGCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaata                                                                          <  2:125734/250‑1 (MQ=255)
                                                                                                                                                                        aGGCTAAGCGGTGTCATCGCAAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAATTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGGACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAatat                                                                      <  2:313401/251‑1 (MQ=255)
                                                                                                                                                                             catggttttcatctctaGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGGTGACACGTTTTTGATCGCTTTTATATAGTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAATCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGGGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAAt                                                                 <  2:302556/235‑1 (MQ=255)
                                                                                                                                                                              aGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATa                                                                <  1:80646/251‑1 (MQ=255)
                                                                                                                                                                                                    aaCCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGta                                             <  1:244537/248‑1 (MQ=255)
                                                                                                                                                                                                                  aaGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAAtt                            >  1:23066/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaGCAACGGGTTATCCCTTTATTTTCCAGACGTGCGTTGACAATTTTTTGATCGATTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCCTACAGGTTACGAAGAATAAAGCTGTTGTATGTTTGGGGTTAAAAAACATAAAACATCCCTTATAAAGATAATTAAAAAAGGGAATAATAAAAAAAAGAAATTATTAATCCCttgtacttcattacattaattt                           >  2:252190/1‑230 (MQ=255)
                                                                                                                                                                                                                                         tttCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGAt      >  1:406878/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ccAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGtta                                                                                          <  1:244306/163‑1 (MQ=255)
                                                                                                                                                                                                                                            ccAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGtta                                                                                          >  2:244306/1‑163 (MQ=255)
                                                                                                                                                                                                                                            ccAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGt  >  1:299590/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
TTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGT  >  NC_000913/155824‑156310

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