Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,986 T→A 7.3% D108V (GAT→GTT)  ydgJ ← putative oxidoreductase YdgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,9860TA 7.3% 141.8 / 2.1 53D108V (GAT→GTT) ydgJputative oxidoreductase YdgJ
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base A (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATC  >  NC_000913/1703847‑1704125
                                                                                                                                           |                                                                                                                                           
gCGGTCAAAATGAGACTCAAAGTAGGCAAAAAATTTCAACGCGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCa                                                                                                                                          <  2:445519/143‑1 (MQ=255)
 cGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAg                                                                                                                                         >  1:248610/1‑143 (MQ=255)
       aaaTGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCGTCCAGCTCTcg                                                                                                                                   >  1:819018/1‑143 (MQ=255)
        aaTGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCCGCCCTcgc                                                                                                                                  >  1:21800/1‑143 (MQ=255)
           gagaCTCGGTGTTAGCAACGTCACCCTGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCTCTATCCCAGTGACGGTTAGACCATCCAGACATCACACGCCCCGCGCTTTTTGCCCGTGCATCCAGCTCTCGCGCt                                                                                                                               <  2:578970/143‑1 (MQ=255)
             gaCTCAAAGTACGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTGTAGCGTCAAGAAATCGCTAGCCCAGCGACGGTTATGGAACACAGACAGCACACGCCCCAGGCCTTTTGCCAGCGCATCCAGCTCCCGCGCTtg                                                                                                                             <  1:169075/143‑1 (MQ=255)
              aCTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCCGGCTTTTTGCCAGCGCATCCAGCACTCGCGCTtgt                                                                                                                            >  2:516928/1‑143 (MQ=255)
                tCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAGGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGa                                                                                                                          <  2:739039/143‑1 (MQ=255)
                tCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCGGCTCTCGCGCTTGTGa                                                                                                                          >  1:365014/1‑143 (MQ=255)
                  aaaGTAAGCAACTTCACCCAGCACGCCTTCCGCGCGCAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACa                                                                                                                        <  1:689056/143‑1 (MQ=255)
                  aaaGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAAAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACa                                                                                                                        >  1:529588/1‑143 (MQ=255)
                                  cccAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAggg                                                                                                        >  1:397229/1‑143 (MQ=255)
                                   ccAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGt                                                                                                       <  2:397229/143‑1 (MQ=255)
                                       cACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGTACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAt                                                                                                   <  1:320673/143‑1 (MQ=255)
                                       cACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTAt                                                                                                   >  1:492248/1‑143 (MQ=255)
                                           ccTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAAc                                                                                               >  1:955340/1‑143 (MQ=255)
                                           ccTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTAGCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAAc                                                                                               <  1:92604/143‑1 (MQ=255)
                                                 gcgAGTAAACCTTGGAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCAAGACTTCCGTCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACcac                                                                                         <  2:938193/143‑1 (MQ=255)
                                                 gcgAGTAAACCTCTGAGCGTCAAGAAACCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACcac                                                                                         <  1:482625/143‑1 (MQ=255)
                                                        aaCCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTa                                                                                  <  2:298800/143‑1 (MQ=255)
                                                          ccTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGGTCTCGCGCTTGTGTCAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTAcc                                                                                >  1:634960/1‑143 (MQ=255)
                                                            ttttAGCCTCAGGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACGCGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCg                                                                              <  2:346872/143‑1 (MQ=255)
                                                             tttAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGc                                                                             <  2:293697/143‑1 (MQ=255)
                                                                 gCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCACCGACCACATGTTTACACGCCTCa                                                                         >  1:356862/1‑143 (MQ=255)
                                                                 gCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCa                                                                         >  2:609221/1‑143 (MQ=255)
                                                                  cGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTcca                                                                        >  1:492116/1‑141 (MQ=255)
                                                                  cGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCaa                                                                        <  1:270451/143‑1 (MQ=255)
                                                                  cGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTCACAGGGGAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCaa                                                                        <  1:193708/143‑1 (MQ=255)
                                                                   gTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCAACCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCacg                                                                       >  2:922151/1‑141 (MQ=255)
                                                                   gTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTGTTGCCAGCGCAACCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAg                                                                       >  2:920469/1‑143 (MQ=255)
                                                                         aactCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCCTCCGCCGTTCTCTCTTGGGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTAACCGCCTCAAGcgccgg                                                                 >  2:884618/4‑142 (MQ=255)
                                                                                   cccGGCGGCGGTGAAGGAAGACAGAAAGCACCCGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTGGTGCCAGTGTCACGGTAAAGGGTTTATCAACGACCCCATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAAc                                                       <  1:506737/143‑1 (MQ=255)
                                                                                    ccAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGCCAGTGTCACGGTAAAGGGTTTCTCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACg                                                      >  1:292933/1‑143 (MQ=255)
                                                                                         gccGGTTATGGAATACATACCGCACACGCCCCAGGCTGCTTGGCCGCGAATCCTGCCCCAGCGATTGTGACAGTCTCACGGTAAAGGGTTTATCAACGACCACACGATTTCCCGCCTCAAGCGCCGCTTTGGCTAACGGGaaa                                                 <  2:497187/141‑1 (MQ=255)
                                                                                                        cagacagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTACCGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGgt                                  <  2:492248/143‑1 (MQ=255)
                                                                                                         agacagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCgccgggcg                                 >  2:675268/1‑136 (MQ=255)
                                                                                                            cagCACCCGCCCCAGGGATGTGGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAg                              <  1:470340/143‑1 (MQ=255)
                                                                                                                  acGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGGATCGTTGGGTGTAGGAATGa                        >  2:697178/1‑143 (MQ=255)
                                                                                                                     ccccAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGATATGGGTATCGTTGGGTGTGGGGATGaata                     >  2:893423/1‑140 (MQ=255)
                                                                                                                           gCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGt               <  1:609221/143‑1 (MQ=255)
                                                                                                                             tttttCCCAGCTCTTCCAGTTCTCGCGCTTGTGGCTTTGTCACGGTAAAGGGTGTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCt             <  1:928310/143‑1 (MQ=255)
                                                                                                                             tttttATCAGCGTATCTAGTTCTCGCTCTTGTGACATTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTAGCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCt             <  1:989257/143‑1 (MQ=255)
                                                                                                                                ttGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATg          >  1:592838/1‑143 (MQ=255)
                                                                                                                                  gCCAGCGCATCCAGCTCTCGCGCTAGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGtt        >  1:798222/1‑143 (MQ=255)
                                                                                                                                   ccAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTg       <  2:955340/143‑1 (MQ=255)
                                                                                                                                    taGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTgg      >  1:60044/2‑143 (MQ=255)
                                                                                                                                    cAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTgg      >  1:967943/1‑143 (MQ=255)
                                                                                                                                    cAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGGGTAGGAATGACAATCAGGGCTATGTTgg      >  2:996587/1‑143 (MQ=255)
                                                                                                                                     ggcgcATCCCGCTCTCGCGCTTGGGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTGCCCGCCTCAAGCTACGCGTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTggg     <  2:938464/142‑1 (MQ=255)
                                                                                                                                      tcgcATCCAGCTCTCTCGTTTGTTTCTGTGTTCTTGTAACTGGTTTATGAACGACCACATGTTTACTCGTCTCAAGCGTCGCTTTGGCTAACGGGAAATGGGTATCGTTGTGTGTAGGAATGACAATCAGGTCTATGTTGGGa    <  2:787615/142‑1 (MQ=255)
                                                                                                                                        gcATGCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATc  <  1:575935/143‑1 (MQ=255)
                                                                                                                                        gcAACCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATc  <  1:920469/143‑1 (MQ=255)
                                                                                                                                        gcAACCAGCTCTCGCGCGTGTGACAGGGGCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATc  <  1:922151/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATC  >  NC_000913/1703847‑1704125

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