Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,045,793 C→T 9.1% L286L (CTG→TTG)  yeeJ → putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,045,7930CT 9.1% 115.6 / 3.8 44L286L (CTG→TTG) yeeJputative adhesin
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base T (2/2);  total (20/24)
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.49e-01

TTAGGTTGGCGTCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAG  >  NC_000913/2045646‑2045940
                                                                                                                                                   |                                                                                                                                                   
ttAGGTTGGCGTCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACGATCTAAAATTAAGCAGGAACGGCTATTTGCGACt                                                                                                                                                    >  2:206187/1‑149 (MQ=255)
    gTTGGCGTCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGATTGAcc                                                                                                                                                >  1:1007445/1‑149 (MQ=255)
    gTTGGCGTCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGATTGAcc                                                                                                                                                >  1:701100/1‑149 (MQ=255)
       ggcggcATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAAc                                                                                                                                             <  2:413180/149‑1 (MQ=255)
          gTCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTgg                                                                                                                                          <  2:344905/149‑1 (MQ=255)
           tCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTgcc                                                                                                                                         >  2:1000057/1‑147 (MQ=255)
             aTTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGgcgc                                                                                                                                       <  2:381854/149‑1 (MQ=255)
             aTTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGgcgc                                                                                                                                       >  1:627422/1‑149 (MQ=255)
               ttCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGAGACTGACAAACTGgcgcag                                                                                                                                     >  2:82252/1‑149 (MQ=255)
                           tcgATGTCGGTCAACACCTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACt                                                                                                                         <  2:968266/147‑1 (MQ=255)
                              aTGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGa                                                                                                                      <  1:848091/149‑1 (MQ=255)
                              aTGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTAATGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTTAACTGGa                                                                                                                      <  2:553804/149‑1 (MQ=255)
                                      cATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAATTTAAGCAGTAACGGCTTTTTGCGACTTACCAACTGGCGCAGCGAACCTGAACTGGACAACGagt                                                                                                              >  2:211731/1‑147 (MQ=255)
                                        tCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTAt                                                                                                            <  1:77131/149‑1 (MQ=255)
                                        tCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTAt                                                                                                            <  1:517598/149‑1 (MQ=255)
                                          aaCTTCTTTTTCGCCCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGCCCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGa                                                                                                          >  1:94584/1‑149 (MQ=255)
                                           aCTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGaa                                                                                                         <  1:556471/149‑1 (MQ=255)
                                                               cTTAGCCGTTACCACTCCCGCGCCGGCATTCGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCAGCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTg                                                                                     <  2:459539/149‑1 (MQ=255)
                                                                ttAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTgg                                                                                    <  2:219756/149‑1 (MQ=255)
                                                                 tAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGCCTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGGCCGACTGGCGCAGCGCACCTGAAGTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTggg                                                                                   >  2:205235/1‑149 (MQ=255)
                                                                           cACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGCCTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACgcgc                                                                         >  2:432672/1‑149 (MQ=255)
                                                                           cACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGCCCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGCAGCACGCCCGGCCAATGGCTGGGATTTACgcgc                                                                         >  2:664611/1‑149 (MQ=255)
                                                                            aCTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCa                                                                        <  1:1000057/149‑1 (MQ=255)
                                                                                ccGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGCACGCGCAGaaa                                                                    <  1:814337/149‑1 (MQ=255)
                                                                                     ccGGCATTTGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATATTCGACTGACCAACTGGCGCAGCGCACATTAACTGGACAACGATTATGAAGCACGCCCGGCCAATTGCTGGGATGTACGCGCAGAAAGCTgg                                                               <  1:763693/149‑1 (MQ=255)
                                                                                     ccGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTgg                                                               <  1:834470/149‑1 (MQ=255)
                                                                                          aTTGGCGCGGAGTACTGGCGCGACTATCTAAAATTGAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTAcc                                                          >  1:305175/1‑149 (MQ=255)
                                                                                          aTTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTAcc                                                          >  1:636874/1‑149 (MQ=255)
                                                                                            tGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGATTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCg                                                        <  2:1007445/149‑1 (MQ=255)
                                                                                            tGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGATTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCg                                                        <  2:701100/149‑1 (MQ=255)
                                                                                              gcgcGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCATTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGcc                                                      <  1:798914/149‑1 (MQ=255)
                                                                                              gcgcGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGcc                                                      <  1:86964/149‑1 (MQ=255)
                                                                                              gcgcGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGcc                                                      >  1:327331/1‑149 (MQ=255)
                                                                                                   gAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGcc                                                 <  2:586656/149‑1 (MQ=255)
                                                                                                     gTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGCCAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGc                                               >  2:281009/1‑149 (MQ=255)
                                                                                                     gTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAACGCTGGCTACCCGCCTGGCCGc                                               >  2:659105/1‑149 (MQ=255)
                                                                                                       aCTGGCGCGACTATGTAAAATTAAGCAGTATCGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGAGTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCAc                                             <  2:330338/149‑1 (MQ=255)
                                                                                                         tGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCt                                           >  1:996764/1‑149 (MQ=255)
                                                                                                              cgACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCg                                      <  1:163491/149‑1 (MQ=255)
                                                                                                                    tCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAAc                                >  2:980539/1‑149 (MQ=255)
                                                                                                                                      cGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatt              <  2:362779/149‑1 (MQ=255)
                                                                                                                                       ggCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta             >  1:491900/1‑149 (MQ=255)
                                                                                                                                                 gACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCCCGCCCGGCCAATGGCTGGGATGTCCGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCGTGGCGGTAAACTGGTCTATGCACAGTATTATGGCGCTGat   >  1:578149/1‑148 (MQ=255)
                                                                                                                                                 gACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGaa   <  1:971857/149‑1 (MQ=255)
                                                                                                                                                  ccTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGTATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAg  <  1:211731/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTAGGTTGGCGTCATTTCACTCCCACATGGATGTCGGGCATCAACTTCTTTTTCGACCACGATCTTAGCCGTTACCACTCCCGCGCCGGCATTGGCGCGGAGTACTGGCGCGACTATCTAAAATTAAGCAGTAACGGCTATTTGCGACTGACCAACTGGCGCAGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAG  >  NC_000913/2045646‑2045940

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