Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,665 G→A 9.1% G254G (GGG→GGA)  yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6650GA 9.1% 120.3 / 2.2 55G254G (GGG→GGA) yfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (19/29);  new base A (2/3);  total (22/33)
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.00e+00

CAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAG  >  NC_000913/2357521‑2357806
                                                                                                                                                |                                                                                                                                              
cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATc                                                                                                                                            <  1:790724/149‑1 (MQ=255)
cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTAATGAGGAGATTGGCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAACGGCATc                                                                                                                                            <  1:797369/149‑1 (MQ=255)
  gccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGATCGGCTCCGGACGGAAGGGCaggaa                                                                                                                                          >  1:351737/1‑145 (MQ=255)
   ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGCTTGCCGAACGTTCACCACTACAAAAGTAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaat                                                                                                                                         >  2:756415/1‑148 (MQ=255)
     agcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAtt                                                                                                                                       <  1:941343/149‑1 (MQ=255)
       cTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCAACAAATtgg                                                                                                                                     >  1:670145/1‑149 (MQ=255)
         aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                                   <  2:727847/149‑1 (MQ=255)
         aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACATCCCGTTTCAGACGTTTTATTGAGGAGATTGCCTAACGTTCTCAACTAAAAAAGGAGAGATTAATGACTATTGCAGAACGGATACTGAATGAATGAAATCAAATtgaat                                                                                                                                   >  1:558107/1‑146 (MQ=255)
           aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctgg                                                                                                                                 <  2:657799/149‑1 (MQ=255)
             acacTGTTAATTTATCTGCTGCAATGCTGCGATACCTCCCGTTTCACCCTTTTTATTGAGGAGATTGCCGAACTTTCACCACTACAAAAGTAGAGATTAATGACTATTGCTGAACGGCTAAGGCAGGAAGGGCATCAAATTGGCTGGCa                                                                                                                               <  2:804714/149‑1 (MQ=37)
                  gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAg                                                                                                                          <  2:407169/149‑1 (MQ=255)
                   tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                         >  1:535028/1‑149 (MQ=255)
                   tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGAGTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgt                                                                                                                         >  1:203291/1‑148 (MQ=255)
                    cttccgaTCTGCTGCATTCCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGACCAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGt                                                                                                                        <  2:223667/143‑1 (MQ=25)
                     tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTa                                                                                                                       >  1:746231/1‑149 (MQ=255)
                     tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACAACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTa                                                                                                                       >  2:932668/1‑149 (MQ=255)
                        ttATCTGCTGCAATGCGGCGATACCTCACGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGc                                                                                                                    <  1:78337/149‑1 (MQ=255)
                          aTCTGCTGCAATGCGGCGATACCGCCCGTTTCACCCGTTTAATTTAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaatta                                                                                                                  <  2:746231/149‑6 (MQ=21)
                              gctgcAATGCGGCGATACCTCCCGTTTAACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGTGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                              >  1:883120/1‑149 (MQ=255)
                               ctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACa                                                                                                             >  2:177158/1‑149 (MQ=255)
                                  cAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑Gc                                                                                                         <  2:24746/149‑1 (MQ=255)
                                   aaTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑Gcc                                                                                                        <  1:266356/149‑1 (MQ=255)
                                       cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGACCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTa                                                                                                    <  1:24813/149‑1 (MQ=255)
                                       cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTa                                                                                                    >  2:966550/1‑149 (MQ=255)
                                        ggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTaa                                                                                                   >  2:380126/1‑149 (MQ=255)
                                        ggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GAGATTaa                                                                                                   >  2:622551/1‑149 (MQ=255)
                                          cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAATCCCATTaaa                                                                                                  >  1:256102/1‑149 (MQ=255)
                                                  cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGCCTATTGCTGCACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATTCATGAACAA‑GAAAATAAAAATGCATTg                                                                                         >  1:127750/1‑149 (MQ=25)
                                                  cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATT‑AAATTGCTTTgc                                                                                        >  1:809195/1‑149 (MQ=255)
                                                                   aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                        <  1:593899/149‑1 (MQ=255)
                                                                       aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGTTCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAC‑GCCATTAAAATTTCTTTGCGCATGCTGGAGCAGGGCttt                                                                    <  2:883120/149‑1 (MQ=255)
                                                                       aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCCTCAAATAGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCttt                                                                    <  2:252588/149‑1 (MQ=255)
                                                                        ggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTg                                                                   >  1:469523/1‑149 (MQ=255)
                                                                           gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                                <  2:916186/149‑1 (MQ=255)
                                                                                   aaCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAtt                                                        >  2:857403/1‑149 (MQ=255)
                                                                                    aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCATTGATCGTGACATtg                                                       >  1:453160/1‑149 (MQ=255)
                                                                                          accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                                 <  1:380126/149‑1 (MQ=255)
                                                                                           ccacTACAAAAGGCGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                                <  1:338311/149‑1 (MQ=255)
                                                                                            cacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAA‑GACATTAAAAATGCTATGCGAATGCTGGACCTCGTCTTTGAACGTGACATTGAGCTGGCg                                               >  1:903386/1‑149 (MQ=21)
                                                                                               tACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGCCa                                            >  2:680331/1‑149 (MQ=255)
                                                                                                 cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCa                                                                                                                               <  1:811986/65‑1 (MQ=14)
                                                                                                 cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCa                                                                                                                               >  2:811986/1‑65 (MQ=14)
                                                                                                          agaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAATGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCAc                                 <  1:563789/149‑1 (MQ=255)
                                                                                                              tAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGAt                             <  2:246544/149‑1 (MQ=255)
                                                                                                               aaTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGAGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATg                            <  2:202322/149‑1 (MQ=255)
                                                                                                                 tGACTATTGCTGACCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTCTGCCTGACCCA‑GCCATTAACATTGCTATGCGCATGCTGGAGCAGGGCTTTGAACTTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCt                          <  1:698465/149‑1 (MQ=37)
                                                                                                                 tGACTATGTATGCACGGCCCCGGCCGGCAGGGCCTCCACTTGGCGGGCAGGAAGGTATGCATGAACAA‑GCCATTAACATTGCTTTGCGCATGCTGTAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCt                          <  2:1021642/149‑1 (MQ=255)
                                                                                                                        tGCTGACCGGCTACGGCATGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCTTGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCATGGCTTTGAACGTGAGATTGTGCTGGCGACAACCAAACTCACTGATGCTGATATTc                   <  1:847363/149‑1 (MQ=38)
                                                                                                                            gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa               <  1:857403/149‑1 (MQ=255)
                                                                                                                            gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGGATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa               <  1:858202/149‑1 (MQ=255)
                                                                                                                                  cTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCa         <  2:203291/149‑1 (MQ=255)
                                                                                                                                    aCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAtt       <  1:966550/149‑1 (MQ=255)
                                                                                                                                        cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCGTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaaa   <  1:756415/149‑1 (MQ=255)
                                                                                                                                        cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATT‑AAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAg  <  2:809195/149‑1 (MQ=255)
                                                                                                                                         aGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAg  <  1:680331/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAA‑GCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAG  >  NC_000913/2357521‑2357806

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: ‑