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→GGAyfaDtransposase_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

CAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAG  >  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)
                                  cAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGc                                                                                                         <  2:24746/149‑1 (MQ=255)
                                   aaTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGcc                                                                                                        <  1:266356/149‑1 (MQ=255)
                                       cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGACCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTa                                                                                                    <  1:24813/149‑1 (MQ=255)
                                       cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTa                                                                                                    >  2:966550/1‑149 (MQ=255)
                                        ggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaa                                                                                                   >  2:380126/1‑149 (MQ=255)
                                        ggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGAGATTaa                                                                                                   >  2:622551/1‑149 (MQ=255)
                                          cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAATCCCATTaaa                                                                                                  >  1:256102/1‑149 (MQ=255)
                                                  cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGCCTATTGCTGCACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATTCATGAACAAGAAAATAAAAATGCATTg                                                                                         >  1:127750/1‑149 (MQ=25)
                                                  cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAATTGCTTTgc                                                                                        >  1:809195/1‑149 (MQ=255)
                                                                   aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                        <  1:593899/149‑1 (MQ=255)
                                                                       aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGTTCATCAAATTGGCTGGCAGGAAGGTATGCATGAACACGCCATTAAAATTTCTTTGCGCATGCTGGAGCAGGGCttt                                                                    <  2:883120/149‑1 (MQ=255)
                                                                       aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCCTCAAATAGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCttt                                                                    <  2:252588/149‑1 (MQ=255)
                                                                        ggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTg                                                                   >  1:469523/1‑149 (MQ=255)
                                                                           gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                                <  2:916186/149‑1 (MQ=255)
                                                                                   aaCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAtt                                                        >  2:857403/1‑149 (MQ=255)
                                                                                    aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCATTGATCGTGACATtg                                                       >  1:453160/1‑149 (MQ=255)
                                                                                          accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                                 <  1:380126/149‑1 (MQ=255)
                                                                                           ccacTACAAAAGGCGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                                <  1:338311/149‑1 (MQ=255)
                                                                                            cacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGACATTAAAAATGCTATGCGAATGCTGGACCTCGTCTTTGAACGTGACATTGAGCTGGCg                                               >  1:903386/1‑149 (MQ=21)
                                                                                               tACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGCCa                                            >  2:680331/1‑149 (MQ=255)
                                                                                                 cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCa                                                                                                                               <  1:811986/65‑1 (MQ=14)
                                                                                                 cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCa                                                                                                                               >  2:811986/1‑65 (MQ=14)
                                                                                                          agaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAATGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCAc                                 <  1:563789/149‑1 (MQ=255)
                                                                                                              tAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGAt                             <  2:246544/149‑1 (MQ=255)
                                                                                                               aaTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATg                            <  2:202322/149‑1 (MQ=255)
                                                                                                                 tGACTATTGCTGACCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTCTGCCTGACCCAGCCATTAACATTGCTATGCGCATGCTGGAGCAGGGCTTTGAACTTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCt                          <  1:698465/149‑1 (MQ=37)
                                                                                                                 tGACTATGTATGCACGGCCCCGGCCGGCAGGGCCTCCACTTGGCGGGCAGGAAGGTATGCATGAACAAGCCATTAACATTGCTTTGCGCATGCTGTAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCt                          <  2:1021642/149‑1 (MQ=255)
                                                                                                                        tGCTGACCGGCTACGGCATGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCTTGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCATGGCTTTGAACGTGAGATTGTGCTGGCGACAACCAAACTCACTGATGCTGATATTc                   <  1:847363/149‑1 (MQ=38)
                                                                                                                            gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa               <  1:857403/149‑1 (MQ=255)
                                                                                                                            gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGGATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa               <  1:858202/149‑1 (MQ=255)
                                                                                                                                  cTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCa         <  2:203291/149‑1 (MQ=255)
                                                                                                                                    aCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAtt       <  1:966550/149‑1 (MQ=255)
                                                                                                                                        cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCGTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaaa   <  1:756415/149‑1 (MQ=255)
                                                                                                                                        cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAg  <  2:809195/149‑1 (MQ=255)
                                                                                                                                         aGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAg  <  1:680331/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAG  >  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: