Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,653 A→G 9.5% L250L (CTA→CTG)  yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6530AG 9.5% 124.3 / 3.3 52L250L (CTA→CTG) yfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (19/27);  new base G (2/3);  total (21/31)
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

GAGGCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATT  >  NC_000913/2357505‑2357801
                                                                                                                                                    |                                                                                                                                                    
gagGCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTa                                                                                                                                                      <  1:694828/149‑1 (MQ=255)
     ttCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCa                                                                                                                                                 >  2:20443/1‑149 (MQ=255)
         cTAATGCCAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGaa                                                                                                                                             >  2:719531/1‑149 (MQ=255)
            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                                                                                     <  2:301920/149‑1 (MQ=255)
                                                                 tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                     <  1:719531/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)
                                                                                   aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCgg                                                                                                                                                   <  2:812463/69‑1 (MQ=255)
                                                                                   aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCgg                                                                                                                                                   >  1:812463/1‑69 (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)
                                                                                                                                                    |                                                                                                                                                    
GAGGCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATT  >  NC_000913/2357505‑2357801

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