Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,094,293 T→C 7.5% E283E (GAA→GAG)  insF1 ← IS3 transposase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,094,2930TC 7.5% 165.9 / 3.6 53E283E (GAA→GAG) insF1IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base T (29/20);  new base C (2/2);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATT  >  NC_000913/1094148‑1094386
                                                                                                                                                 |                                                                                             
gCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCa                                                                                            >  4:290197/1‑149 (MQ=255)
  tGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCaaa                                                                                          <  2:342025/149‑1 (MQ=255)
  tGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGCCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCaaa                                                                                          >  2:106816/1‑149 (MQ=255)
   gATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAt                                                                                         >  2:42760/1‑149 (MQ=255)
   gATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAt                                                                                         >  4:487841/1‑149 (MQ=255)
     tAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTg                                                                                       <  3:311187/149‑1 (MQ=255)
         gCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGCCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTcc                                                                                   <  1:106816/149‑1 (MQ=255)
           aCTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCgg                                                                                 >  3:723994/1‑149 (MQ=255)
             tttACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTCTCAAATTGTTCCGGAc                                                                               >  4:399323/1‑149 (MQ=255)
             tttACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTCTCAAATTGTTCCGGAc                                                                               >  2:258037/1‑149 (MQ=255)
                 cGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                           <  1:136221/149‑1 (MQ=255)
                 cGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                           <  3:290197/149‑1 (MQ=255)
                   aGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgc                                                                         <  1:682015/149‑1 (MQ=255)
                    ggCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTCTCAAATTGTTCCGGACTGAGgcc                                                                        <  3:399323/149‑1 (MQ=255)
                    ggCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTCTCAAATTGTTCCGGACTGAGgcc                                                                        <  1:258037/149‑1 (MQ=255)
                      cGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                                      >  2:659998/1‑149 (MQ=255)
                       gAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgcc                                                                     >  2:242954/1‑149 (MQ=255)
                       gAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgcc                                                                     >  3:705682/1‑149 (MQ=255)
                       gAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgcc                                                                     >  4:51726/1‑149 (MQ=255)
                            agaCCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACAcc                                                                >  4:196706/1‑149 (MQ=255)
                               ccGGGCGTAATAAGGTGGTTGTGAGTGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAGTGTTCCGGACTGAGGCCGCCACACCAAc                                                             <  4:139170/149‑1 (MQ=255)
                               ccGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAAc                                                             <  3:196706/149‑1 (MQ=255)
                                 gggCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTATGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtg                                                           >  4:464164/1‑149 (MQ=255)
                                  ggCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtgt                                                          <  2:652051/149‑1 (MQ=255)
                                  ggCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtgt                                                          <  4:364662/149‑1 (MQ=255)
                                         taaGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgccg                                                   >  3:252941/1‑149 (MQ=255)
                                         taaGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgccg                                                   >  2:754874/1‑149 (MQ=255)
                                         taaGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgccg                                                   >  2:612709/1‑149 (MQ=255)
                                             gtggtTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                               >  2:605427/1‑149 (MQ=255)
                                                  tgtgAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGAtt                                          <  1:242954/149‑1 (MQ=255)
                                                          aTGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATcac                                  >  2:157968/1‑149 (MQ=255)
                                                           tgtgGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                     <  1:260882/145‑1 (MQ=255)
                                                           tgtgGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                     >  2:260882/1‑145 (MQ=255)
                                                               gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTcc                                                                                   <  4:595879/95‑1 (MQ=255)
                                                               gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTcc                                                                                   >  3:595879/1‑95 (MQ=255)
                                                               gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCg                             >  4:269342/1‑149 (MQ=255)
                                                               gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCg                             >  2:126052/1‑149 (MQ=255)
                                                                tgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa                            >  1:330664/1‑149 (MQ=255)
                                                                    tttttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGAtt                                          <  4:369413/131‑1 (MQ=255)
                                                                    tttttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGAtt                                          >  3:369413/1‑131 (MQ=255)
                                                                          tgAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTaaa                  <  2:379997/149‑1 (MQ=255)
                                                                              tCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACAcc              <  2:171656/149‑1 (MQ=255)
                                                                                aCCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGt            >  3:780454/1‑149 (MQ=255)
                                                                                  cATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTg          >  3:165294/1‑149 (MQ=40)
                                                                                   aTAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGc         >  4:720176/1‑149 (MQ=40)
                                                                                   aTAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGc         >  4:33131/1‑149 (MQ=40)
                                                                                    tAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGcc        >  2:743935/1‑149 (MQ=39)
                                                                                    tAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGc         >  1:675835/1‑148 (MQ=39)
                                                                                    tAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGc         <  2:675835/148‑1 (MQ=39)
                                                                                        gCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCa    <  1:612709/149‑1 (MQ=35)
                                                                                        gCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCa    <  1:754874/149‑1 (MQ=35)
                                                                                        gCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCa    <  3:51726/149‑1 (MQ=35)
                                                                                          ggCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCatt  >  4:206764/1‑149 (MQ=34)
                                                                                                                                                 |                                                                                             
GCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATT  >  NC_000913/1094148‑1094386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑