Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,665 G→A 12.3% 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,6650GA12.3% 98.4 / 4.2 50G254G (GGG→GGA) yfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (12/29);  new base A (3/3);  total (18/32)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGA  >  NC_000913/2357517‑2357812
                                                                                                                                                    |                                                                                                                                                   
aTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                                     <  2:20575/149‑1 (MQ=255)
  gACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACATCCCGTTTCACCCGGTTTATTGAGTAGATTGCCGAACGTTCACCACTACAAAAGAAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCa                                                                                                                                                   >  1:449848/1‑149 (MQ=39)
   aCAGCCAGCTACAAACACTGTTTATTTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCaa                                                                                                                                                  >  2:574673/1‑148 (MQ=255)
       ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATCCCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                              >  1:359984/1‑149 (MQ=255)
       ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                              >  1:20753/1‑149 (MQ=255)
       ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGATTCACCCGTTTTAATGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                              <  1:403998/149‑1 (MQ=255)
        catcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACACGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGATATATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAt                                                                                                                                             <  2:212887/146‑1 (MQ=255)
         agcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAtt                                                                                                                                            >  2:924520/1‑149 (MQ=255)
         agcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACACGTTTTATTGAGGAGATTGCCGAACGTACACCACTACAAAAGGAGAGATTAATGACTATTGATGAACGGCTACGGCAGGAAGGGCATCAAAtt                                                                                                                                            >  2:106401/1‑149 (MQ=255)
             acccccacTGGTTAATTCTCTGCTGCAATGCGGCGATACCTCCCGGTTCACCCGTGCTATTGCGGAGATTGCCGACCGTTCACCACTACAACAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                                        <  2:946901/144‑1 (MQ=37)
              cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGACTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                                       <  2:230959/149‑1 (MQ=255)
                     tttttAATTATCAGCTGCAATGCGGCGATACCTCGCGTTGCACCCGTTTTAATGACGAGATTGCCGTACGTGCACCACTTCAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                                                <  1:319180/147‑1 (MQ=25)
                      gTTTAATTATCTGGTGCTATGAGGCGATAACTACCGGTACACCCCGTTTATATAGGAGACTGCAGAACGTTAACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAg                                                                                                                               <  1:685612/149‑1 (MQ=17)
                       tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaaa                                                                                                                                              <  1:382678/133‑1 (MQ=255)
                       tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaaa                                                                                                                                              >  2:382678/1‑133 (MQ=255)
                        ttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGCGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                             <  2:973077/149‑1 (MQ=255)
                        ttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                             <  2:551425/149‑1 (MQ=255)
                              aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACCGCACTAAAGAAATCAAATTGCATGCCACGAAGTTaaaaaa                                                                                                                       >  2:89433/1‑144 (MQ=37)
                                   ctgcAATGCGGCGATACCTCCCGTTTCACCCGTTGTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGCTTTATGACGATTGCTGAACGGCGACGGCAGGAAGGGCAGCAAATTGGCTGGCAGGAAGGTATGCATGAACa                                                                                                                  >  2:584733/1‑149 (MQ=38)
                                         tGCGGCGATACCTCCCGTTTCACCCGTATTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATAGGCTGGCAGGAAGGTATGCATGAACAAGCCAt                                                                                                            >  1:695565/1‑149 (MQ=255)
                                         agcggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACAACAAAACGAGAGATTAATGACTAATGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAt                                                                                                            <  2:634575/148‑1 (MQ=255)
                                           cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTa                                                                                                          <  2:1031604/149‑1 (MQ=255)
                                              cgATACCTCCCGTTTCACACGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTTCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                                       <  1:721360/149‑1 (MQ=255)
                                                     tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                >  1:589503/1‑149 (MQ=255)
                                                     tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                <  1:719531/149‑1 (MQ=255)
                                                     tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                <  2:301920/149‑1 (MQ=255)
                                                      cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTg                                                                                               <  1:509427/149‑1 (MQ=255)
                                                       ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                                              <  2:529200/149‑1 (MQ=255)
                                                             cACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGGATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGc                                                                                        <  1:831854/149‑1 (MQ=255)
                                                                           aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCttt                                                                          <  2:589503/149‑1 (MQ=255)
                                                                              agaTTGCCGAACGTTCCCCAGTACAAAAGGAGAGATTAATGACTATTGCTGCACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGaa                                                                       <  1:563704/149‑1 (MQ=255)
                                                                                aTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAAGAACAAGACATTAAAATTGCTATGAGAATTCTTTAGCAGAGCTTTGATCg                                                                     >  2:15979/1‑149 (MQ=21)
                                                                                 ttGCCGAACGTTCACCCCTACAAAAGGAGAGATTACTGACTATTGCTGCACGGCTACGGCAGGCAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGt                                                                    <  1:139351/149‑1 (MQ=255)
                                                                                  tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGATTGAACAATGCATTTAAAATGCTTTTACCATTCTGTAGCACGGGTTTGAACGTg                                                                   >  1:850513/1‑149 (MQ=18)
                                                                                   gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACTAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                  <  1:574673/149‑1 (MQ=255)
                                                                                     cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaat                                                                                                                         >  1:756258/1‑89 (MQ=255)
                                                                                       aaCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAtt                                                              <  1:40448/149‑1 (MQ=255)
                                                                                          gcgcaccacTACACAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATtgtg                                                           <  1:655417/146‑1 (MQ=255)
                                                                                           ttCACCACTTCAAAAGGAGATATTAATGACTATTGCTGAAAGGCTACTGGAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGc                                                          <  2:83839/149‑1 (MQ=25)
                                                                                                  cTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAg                                                                                                                               <  1:591485/73‑1 (MQ=9)
                                                                                                  cTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAg                                                                                                                               >  2:591485/1‑73 (MQ=255)
                                                                                                    aCAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAAGAAGCCATTAAAATTGCTTTGCTCATTCCTGAGCTGAGCTTTGATCTTGAATTTCTGCTggacgaag                                                 >  2:122755/1‑143 (MQ=18)
                                                                                                      aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAcc                                               >  2:1001126/1‑149 (MQ=255)
                                                                                                             gagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGTCTGGCAGGAAGTTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTTTATCAGGGCTTTGAACGTGAGATTGTGCTGGCTACAACCCAACTCa                                        <  2:449848/149‑1 (MQ=37)
                                                                                                                aTTAATGACTCTTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTCTGCATGAACAAGCCAATAAAACTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTg                                     <  1:584733/149‑1 (MQ=255)
                                                                                                                               tGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGa                      >  2:566384/1‑149 (MQ=255)
                                                                                                                                gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa                     <  1:452507/149‑1 (MQ=255)
                                                                                                                                    ggCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAATGTATGCATGAACAAGCCATTAAATTTGCTTTGCGCATGCTGGAGCAGGGCTGTGAACGTGAGATTGTGCGGGCGACAACCCAACTAACTGATGCTGATATTCCGAACTGt                 <  1:188756/149‑1 (MQ=38)
                                                                                                                                        aCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATTCTGGAGCAGGGATTTGAACTTGAGATTGTGATGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAtt             <  2:53673/149‑1 (MQ=255)
                                                                                                                                            cAGGAAGGGCATCAAATTGTCTGGCAGGAAGCAATGCATGAACAAGCCATTAAAATTGCTTTGCGCATCCTTGAGCAGGGCTTTTAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaaa         <  2:797977/149‑1 (MQ=38)
                                                                                                                                            cAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTaaaa         <  1:20443/149‑1 (MQ=255)
                                                                                                                                                   ggCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgagaga  <  1:924520/149‑1 (MQ=255)
                                                                                                                                                   ggCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgagaga  >  1:243342/1‑149 (MQ=255)
                                                                                                                                                   ggCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgagaga  <  1:106401/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGA  >  NC_000913/2357517‑2357812

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