Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,822,446 T→A 8.8% E65V (GAA→GTA) ‡ recX ← regulatory protein for RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,822,4460TA 8.8% 187.6 / 8.4 57E65V (GAA→GTA) ‡recXregulatory protein for RecA
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base A (2/3);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TCACATTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCA  >  NC_000913/2822299‑2822589
                                                                                                                                                   |                                                                                                                                               
tCACATTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTc                                                                                                                                                >  4:26424/1‑149 (MQ=255)
   caTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATg                                                                                                                                             <  3:39376/149‑1 (MQ=255)
    aTTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATgg                                                                                                                                            <  1:149650/149‑1 (MQ=255)
        aCGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCAcc                                                                                                                                        <  3:54717/149‑1 (MQ=255)
         cGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCa                                                                                                                                       <  1:400671/149‑1 (MQ=255)
            aTCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGc                                                                                                                                    >  2:80923/1‑149 (MQ=255)
             tCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCAt                                                                                                                                              <  2:394/138‑1 (MQ=255)
             tCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTACGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCAt                                                                                                                                              >  1:394/1‑138 (MQ=255)
                cTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATa                                                                                                                                >  2:438050/1‑149 (MQ=255)
                                      aaTACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCg                                                                                                          >  1:39755/1‑149 (MQ=255)
                                            tttCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATCAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGcc                                                                                                           <  4:342639/142‑1 (MQ=255)
                                            tttCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGcc                                                                                                           >  3:342639/1‑142 (MQ=255)
                                            tttCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCAt                                                                                                    >  2:302823/1‑149 (MQ=255)
                                                 gATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAAtc                                                                                               >  2:467246/1‑149 (MQ=255)
                                                  aTTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAAtct                                                                                              >  4:300937/1‑149 (MQ=255)
                                                   ttCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAg                                                                                                                                      <  1:198948/108‑1 (MQ=255)
                                                   ttCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAg                                                                                                                                      >  2:198948/1‑108 (MQ=255)
                                                        ttCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTg                                                                                        <  2:77783/149‑1 (MQ=255)
                                                         tCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTgg                                                                                       <  1:438050/149‑1 (MQ=255)
                                                            tGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGcc                                                                                    >  3:416991/1‑149 (MQ=255)
                                                                gAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAttt                                                                                >  3:410171/1‑149 (MQ=255)
                                                                  aTACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAttttt                                                                              >  3:193456/1‑149 (MQ=255)
                                                                    aCGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGc                                                                            <  2:39755/149‑1 (MQ=255)
                                                                    aCGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGc                                                                            <  1:302823/149‑1 (MQ=255)
                                                                          aGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATaa                                                                      >  3:413571/1‑149 (MQ=255)
                                                                                 tAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAttttt                                                                              >  3:355132/1‑134 (MQ=255)
                                                                                 tAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCAttttt                                                                              <  4:355132/134‑1 (MQ=255)
                                                                                 tAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGc                                                               <  3:300937/149‑1 (MQ=255)
                                                                                              tACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAAt                                                                     >  4:13241/1‑130 (MQ=255)
                                                                                              tACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAAt                                                                     <  3:13241/130‑1 (MQ=255)
                                                                                              tACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTa                                                     >  2:302575/1‑146 (MQ=255)
                                                                                              tACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTa                                                     <  1:302575/146‑1 (MQ=255)
                                                                                                         aTAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCGGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTg                                       <  4:410171/149‑1 (MQ=255)
                                                                                                         aTAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTg                                       <  2:49367/149‑1 (MQ=255)
                                                                                                         aTAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTg                                       >  1:49367/1‑149 (MQ=255)
                                                                                                         aTAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTg                                       <  3:60866/149‑1 (MQ=255)
                                                                                                          tAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATGACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGTGAGTTTACGTCGCAGTTCTTGc                                      >  1:408532/1‑149 (MQ=255)
                                                                                                                       aaaTCTGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATcgc                         <  4:416991/149‑1 (MQ=255)
                                                                                                                            ggCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAcc                    >  1:405204/1‑149 (MQ=255)
                                                                                                                              cTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGc                  <  4:413571/149‑1 (MQ=255)
                                                                                                                                 tCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAg               <  1:53187/149‑1 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATcgc                         >  1:300528/1‑137 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATcgc                         <  2:300528/137‑1 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAc                     <  1:13208/141‑1 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAc                     >  4:109801/1‑141 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAc                     >  2:13208/1‑141 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAc                     <  3:109801/141‑1 (MQ=255)
                                                                                                                                   aTCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGaa             >  3:447313/1‑149 (MQ=255)
                                                                                                                                    tCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAAt            >  4:135090/1‑149 (MQ=255)
                                                                                                                                      gagaTAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGc          <  1:71547/149‑1 (MQ=255)
                                                                                                                                       ctggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTAc                                                    >  3:404239/6‑106 (MQ=255)
                                                                                                                                       ctggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTAc                                                    <  4:404239/101‑1 (MQ=255)
                                                                                                                                        tggtgCCATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGt        >  4:469813/5‑149 (MQ=255)
                                                                                                                                         aTAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTa       <  3:135090/149‑1 (MQ=255)
                                                                                                                                             ccATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGc   <  1:163469/149‑1 (MQ=255)
                                                                                                                                             ccATGAACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGc   <  3:469813/149‑1 (MQ=255)
                                                                                                                                              cATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCa  >  1:317504/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCACATTCACGCATCGCTTTTTCTGTCGCTTCGCGGGAAATACCTTTCTGATTCAGTTCCTGGCGAATACGCGCAGGTCCATAACCTTTGCGGCTACGGCTGGCGATAAAGCGCGCAACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCA  >  NC_000913/2822299‑2822589

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