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,0680AG94.7% 49.9 / ‑6.3 19V45A (GTA→GCA) yadVputative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (8/10);  minor base T (1/0);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGTAATG  >  NC_000913/155823‑156314
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
tttCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGCCACGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTg                                                                                                                                                                                                                                                   >  2:297830/1‑251 (MQ=255)
           ggtttcggtagtaCTTCCAGTACGGTTAACCAGACCACGCTCTCTCTGTTTTTAGGCATTGAGGTGCTTGCTGTTTACATTAATTTGATTGTTTGCCCACGTTTGGCAGCAAGTCGCGATACTTGTGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGAGTTGTCATCGCCATTATCTAACCATCTCTGGACCAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACtttttt                                                                                                                                                                                                                                         <  2:316219/238‑1 (MQ=255)
                      aaCTTCCAGTACGTTAAAACAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATAGCAAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGTCCCTTTTTTGATCGCTTTtat                                                                                                                                                                                                                             >  2:336415/1‑251 (MQ=255)
                       aCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTtata                                                                                                                                                                                                                            <  1:243370/251‑1 (MQ=255)
                               tACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGATGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTCAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCCTCGCCAGTATCTAACCAAATCTTGACAAGCAACGGGGTATTCCCTTTATTTTCCAGACGTTCGTTGCCACTTTTTTCATCGCTTTTTTATAtttctc                                                                                                                                                                                                                    >  2:225/1‑247 (MQ=255)
                                aCGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTAcgc                                                                                                                                                                                                                    <  1:297830/250‑1 (MQ=255)
                                      aaCCAGAACCCGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGCCTGTAATACTGCCAGGCTCAGTGTTGTCATCGCCAGTATCTAACCAGCTCTGGACAAGCAACGGGTTATTCCCGTTATTGTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTAcc                                                                                                                                                                                                             <  1:44691/251‑1 (MQ=255)
                                                                                                            cACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTa                                                                                                                                       <  1:336415/251‑1 (MQ=255)
                                                                                                                       cacagtcgttagacgtgaggtggaGCAGTAAAATGGACTTTAATACTGTTAGTCTCAGCGTTGTCAACGCCAGTATCTAACCAACTGTGGACAAGCAACTTGTTATGACCTTTATTTTCCAGACGTGCTTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGCAATGAAAATGTACGCAATTGACGATGATCTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTGTGGTGtt                                                                                                                            <  2:256005/228‑1 (MQ=255)
                                                                                                                            aCGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTaaaaa                                                                                                                       <  1:225/251‑1 (MQ=255)
                                                                                                                               cgATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAAc                                                                                                                     >  1:206674/1‑250 (MQ=255)
                                                                                                                                                    aaaagggctgttaattccgcccggCTCAGCTTTGTCATCGCCAGTATCTAACCACCTCTGGACAAGCAACGGGTTATTCCCTTTTTTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTTAGCGAGTACACGAAATGACAATGTCCGCAATAGACGATTAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGtt                                                                                               <  2:322297/229‑1 (MQ=255)
                                                                                                                                                    aaaaGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATTCACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTTAAAAACATAAACCATCCATGATAGAGtt                                                                                               >  2:83133/1‑251 (MQ=255)
                                                                                                                                                                aaTACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAg                                                                                    <  1:299298/250‑1 (MQ=255)
                                                                                                                                                                   aCTGCCAGGCTCAGCGTTGTCAACGCCAGTATCTAACCAACACTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGa                                                                                 <  1:83133/250‑1 (MQ=255)
                                                                                                                                                                         aGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAatat                                                                          >  1:235296/1‑251 (MQ=255)
                                                                                                                                                                                                 aTCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGt                                                  >  1:319905/1‑251 (MQ=255)
                                                                                                                                                                                                           ctctGGACAAGAAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACAATCCCTGATAGCGTTAATTAAAACAGGGCATAATATAAAATAGAAATTATTATTCCCTGTATattcattc                                        >  2:150535/1‑251 (MQ=255)
                                                                                                                                                                                                                                ttCCCTTTATATTCCAGACGTGAGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGATTGTTTGGTGTTAAAAAAAATAAACCATCCCCGATAGCGATAATTAAATCAGGACAAAATATAAAACAAGATTTATTATTCCATGTATATTCATTCAATCAAATTAACTGGCGAAt                    >  2:305017/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 acgtgcgtTGACTCTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCAGCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGTAATg  <  1:39520/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    tacgtTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaa                                                                                >  1:55927/1‑170 (MQ=255)
                                                                                                                                                                                                                                                    tacggTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaa                                                                                <  2:55927/170‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
TTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGTAATG  >  NC_000913/155823‑156314

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