Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,582,843 A→G V85A (GTT→GCT)  ydeS ← fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,582,8430AG93.8% 95.3 / ‑0.3 32V85A (GTT→GCT) ydeSfimbrial protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (13/17);  total (14/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 = 8.61e-01

AAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGA  >  NZ_CP009273/1582707‑1582978
                                                                                                                                        |                                                                                                                                       
aaGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGcc                                                                                                                                       <  1:161894/139‑1 (MQ=255)
aaGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGcc                                                                                                                                       <  2:1260982/139‑1 (MQ=255)
      aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGt                                                                                                                                  <  1:265523/138‑1 (MQ=255)
             ttCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCtt                                                                                                                          >  1:901715/1‑139 (MQ=255)
                 ccGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTtgct                                                                                                                      >  1:631271/1‑139 (MQ=255)
                    gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCa                                                                                                                   <  2:474083/139‑1 (MQ=255)
                      cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCAtt                                                                                                                 >  1:663769/1‑139 (MQ=255)
                                  tCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTaa                                                                                                     >  2:379849/1‑139 (MQ=255)
                                   ccAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAAt                                                                                                    >  2:230908/1‑139 (MQ=255)
                                                               cGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTa                                                                        <  2:484175/139‑1 (MQ=255)
                                                                gTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTAc                                                                       <  2:901715/139‑1 (MQ=255)
                                                                 tATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGt                                                                                                                                  >  2:70682/1‑79 (MQ=255)
                                                                 tATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGt                                                                                                                                  <  1:70682/79‑1 (MQ=255)
                                                                              aTCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGaaa                                                          <  2:232877/138‑1 (MQ=255)
                                                                                        aacaaAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCa                                                <  1:379849/138‑1 (MQ=255)
                                                                                             aGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTGGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGttt                                          <  2:921697/139‑1 (MQ=255)
                                                                                                       tCTGCATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTg                                     >  2:655668/1‑134 (MQ=255)
                                                                                                           cATCCTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATc                               <  2:976100/136‑1 (MQ=255)
                                                                                                              ccTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTaa                                                                                                     <  2:21081/63‑1 (MQ=255)
                                                                                                              ccTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTaa                                                                                                     >  1:21081/1‑63 (MQ=255)
                                                                                                               cTCAATACCGTTAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGt                         <  2:631271/138‑1 (MQ=255)
                                                                                                                         ttAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGc               >  2:99424/1‑138 (MQ=255)
                                                                                                                         ttAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGc               >  2:246553/1‑138 (MQ=255)
                                                                                                                         ttAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGc               >  1:796212/1‑138 (MQ=255)
                                                                                                                         ttAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGc               >  1:206808/1‑138 (MQ=255)
                                                                                                                          tAAATGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCg              >  1:150918/1‑138 (MQ=255)
                                                                                                                            aaTGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGaa            <  2:48639/138‑1 (MQ=255)
                                                                                                                            aaTGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGaa            <  2:206808/138‑1 (MQ=255)
                                                                                                                            aaTGCAACCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGaa            <  1:1025206/138‑1 (MQ=255)
                                                                                                                                 aaCCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACtt                                                                              <  1:928214/67‑1 (MQ=255)
                                                                                                                                 aaCCCGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACtt                                                                              >  2:928214/1‑67 (MQ=255)
                                                                                                                                     cGAGCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATgaga  <  1:728849/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGA  >  NZ_CP009273/1582707‑1582978

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