Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,754,310 C→A R382M (AGG→ATG)  yiaY ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,754,3100CA99.1% 320.5 / ‑5.5 107R382M (AGG→ATG) yiaYputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (54/52);  total (55/52)
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.66e-01

Only 100 of 107 total aligned reads displayed.

ACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTG‑‑CATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTT  >  NC_000913/3754238‑3754380
                                                                        |                                                                         
aCGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCg                                                                        >  2:3204650/1‑76 (MQ=255)
aCGGGCAAGGGCTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCg                                                                        >  2:1971792/1‑76 (MQ=255)
  gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCt                                                                      <  2:4696980/76‑1 (MQ=255)
  gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCt                                                                      >  2:3842367/1‑76 (MQ=255)
  gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCt                                                                      <  1:2778447/76‑1 (MQ=255)
  gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCt                                                                      >  1:539117/1‑76 (MQ=255)
   ggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                                                                     >  2:4076227/1‑76 (MQ=255)
    gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑c                                                                  >  1:17198/1‑76 (MQ=255)
     cAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAGCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cg                                                                 <  1:4318877‑M1/76‑2 (MQ=255)
      aaGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑ccc                                                                >  2:133942/1‑74 (MQ=255)
      aaGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgc                                                                >  1:1477848‑M1/1‑74 (MQ=255)
      aaGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgc                                                                <  2:4556567‑M1/76‑3 (MQ=255)
       aGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcg                                                               >  1:1522653‑M1/1‑73 (MQ=255)
         ggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgat                                                             <  1:1438291‑M1/76‑6 (MQ=255)
            tACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaa                                                          >  1:766626‑M1/1‑68 (MQ=255)
            tACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaa                                                          <  2:4648615‑M1/76‑9 (MQ=255)
            tACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaa                                                          <  1:1665417‑M1/76‑9 (MQ=255)
             aCAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaat                                                         >  1:3321013‑M1/1‑67 (MQ=255)
             aCAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaat                                                         <  1:3214684‑M1/76‑10 (MQ=255)
              cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatc                                                        <  1:536374‑M1/76‑11 (MQ=255)
              cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatc                                                        <  1:3845614‑M1/76‑11 (MQ=255)
              cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatc                                                        >  1:3725158‑M1/1‑66 (MQ=255)
              cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatc                                                        >  1:3440743‑M1/1‑66 (MQ=255)
               agagATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcg                                                       >  2:2943249‑M1/1‑65 (MQ=255)
                tagatAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgc                                                      >  2:353335‑M1/2‑64 (MQ=255)
                  gaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgcca                                                    >  2:1121179‑M1/1‑62 (MQ=255)
                  gaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgcca                                                    >  1:1129562‑M1/1‑62 (MQ=255)
                  gaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgcca                                                    <  2:1671099‑M1/76‑15 (MQ=255)
                   aTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccac                                                   <  1:2426054‑M1/76‑16 (MQ=255)
                     aGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaa                                                 >  2:4384743‑M1/1‑59 (MQ=255)
                     aGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaa                                                 >  1:4452739‑M1/1‑59 (MQ=255)
                      gTTGAAAAAGCGGCTAACAATTGGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaat                                                <  2:3713683‑M1/76‑19 (MQ=255)
                       ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatt                                               >  2:3874819‑M1/1‑57 (MQ=255)
                       ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatt                                               >  1:2883835‑M1/1‑57 (MQ=255)
                        tGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaattt                                              >  2:4402639‑M1/1‑56 (MQ=255)
                          aaaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttct                                            <  2:178706‑M1/76‑23 (MQ=255)
                          aaaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttct                                            <  2:1641586‑M1/76‑23 (MQ=255)
                           aaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCt                                                                      >  1:1244300/1‑51 (MQ=255)
                           aaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCt                                                                      <  2:1244299/51‑1 (MQ=255)
                           aaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttctt                                           >  2:4157667‑M1/1‑53 (MQ=255)
                             aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttctttg                                         >  1:3982915‑M1/1‑51 (MQ=255)
                             aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcg                                         <  2:462822‑M1/76‑26 (MQ=255)
                             aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcg                                         >  1:3948557‑M1/1‑51 (MQ=255)
                              aGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgatgactcgccccaatttcttcgt                                        >  2:4033164‑M1/1‑50 (MQ=255)
                              aGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgt                                        >  1:1900584‑M1/1‑50 (MQ=255)
                              aGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgt                                        >  1:2975176‑M1/1‑50 (MQ=255)
                              aGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTA‑‑Cgcgataaatcgccacaatttcctcgt                                        >  2:4028030‑M1/1‑50 (MQ=255)
                                cGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtga                                      <  1:550160‑M1/76‑29 (MQ=255)
                                 ggCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgag                                     <  1:1589435‑M1/76‑30 (MQ=255)
                                  gCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagt                                    >  1:4631057‑M1/1‑46 (MQ=255)
                                   cTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagtt                                   >  1:332536‑M1/1‑45 (MQ=255)
                                   cTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagtt                                   <  1:3615502‑M1/76‑32 (MQ=255)
                                     aaaaatttGCCAGCAGTTGTGGAAAGGATGATTACATCGCGG‑‑Cgcgataaatcgccacaatttcttcgtgagttgc                                 <  1:2561106‑M1/73‑34 (MQ=255)
                                     aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgcgctgc                                 >  2:258672‑M1/1‑43 (MQ=255)
                                     aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgc                                 <  1:3007587‑M1/76‑34 (MQ=255)
                                     aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgc                                 <  1:2515177‑M1/76‑34 (MQ=255)
                                     aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgc                                 <  1:1556862‑M1/76‑34 (MQ=255)
                                      aCAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgaggtgcc                                >  2:2222942‑M1/1‑42 (MQ=255)
                                      aCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcc                                >  2:2205392‑M1/1‑42 (MQ=255)
                                        aaTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctg                              <  2:2883835‑M1/76‑37 (MQ=255)
                                        aaTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctg                              >  1:841216‑M1/1‑40 (MQ=255)
                                        aaTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctg                              >  1:3695225‑M1/1‑40 (MQ=255)
                                         aTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctgg                             >  1:2240089‑M1/1‑39 (MQ=255)
                                         aTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctgg                             <  2:4185035‑M1/76‑38 (MQ=255)
                                         aTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctgg                             >  2:2573645‑M1/1‑39 (MQ=255)
                                          tttGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctgga                            <  1:1145620‑M1/76‑39 (MQ=255)
                                          gttGCCAGCCGTTGGGGTAACGAAGATTACATCGCTGAACGATAACACGACACAACTTCGTCGTGAGTTGCCTGGa                              >  2:3996056/2‑76 (MQ=255)
                                           ttGCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggat                           <  1:2709292‑M1/76‑40 (MQ=255)
                                             gCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcg                         >  1:2082038‑M1/1‑35 (MQ=255)
                                             gCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcg                         >  1:4535508‑M1/1‑35 (MQ=255)
                                             gCCAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcg                         <  2:940260‑M1/76‑42 (MQ=255)
                                              ccAGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgg                        <  2:4074760‑M1/76‑43 (MQ=255)
                                              ccAGCCGTTGTGGAAATAATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgg                        >  1:3578029‑M1/1‑34 (MQ=255)
                                                ggcCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggt                      <  2:950474‑M1/75‑45 (MQ=255)
                                                agtcgttGTGGATATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggt                      <  1:379769‑M1/73‑45 (MQ=255)
                                                aggcgttGTGGAAATGTTGCTTACATCGCGG‑‑Cgggataaatcgccacaatttcttcgtgagttgcctggatcgggt                      <  2:1551954‑M1/73‑45 (MQ=255)
                                                aGCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggacagggt                      >  2:4739270‑M1/1‑32 (MQ=255)
                                                aGCCGGTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggt                      <  2:1829569‑M1/76‑45 (MQ=255)
                                                 gCCGTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggtt                     <  1:4683465‑M1/76‑46 (MQ=255)
                                                  ccGTTGTGGAAATGATGATTACATCGCTG‑‑Cacgataaatcgccacaatttcttcgtgagttgcctggatcgggtta                    <  2:2815067‑M1/76‑47 (MQ=255)
                                                    gTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagt                  <  2:2376510‑M1/76‑49 (MQ=255)
                                                    gTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagt                  >  1:3009108‑M1/1‑28 (MQ=255)
                                                    gTTGTGGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggaccgggttagt                  >  2:2134664‑M1/1‑28 (MQ=255)
                                                        tgGAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaag              <  1:4466076‑M1/76‑53 (MQ=255)
                                                          gAAATGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcc            >  2:2266733‑M1/1‑22 (MQ=255)
                                                           aaaTGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcca           >  2:4558577‑M1/1‑21 (MQ=255)
                                                           aaaTGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcca           <  2:3520047‑M1/76‑56 (MQ=255)
                                                            aaTGATGATTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccac          >  1:11585‑M1/1‑20 (MQ=255)
                                                            aaTGATGATTACATCGATG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgtgttagtcaagcaca          >  1:45833‑M1/1‑20 (MQ=255)
                                                              tgatgaTTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacag        <  1:4076227‑M1/76‑59 (MQ=255)
                                                                atgatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggc      >  1:2526291‑M1/1‑16 (MQ=255)
                                                                 tgatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca     <  2:4758701‑M1/76‑62 (MQ=255)
                                                                 tgatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca     <  2:3253929‑M1/76‑62 (MQ=255)
                                                                 tgatTACATCGCGG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca     <  2:551533‑M1/76‑62 (MQ=255)
                                                                  gatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcat    <  1:4384743‑M1/76‑63 (MQ=255)
                                                                  gatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcat    <  2:4332128‑M1/76‑63 (MQ=255)
                                                                  gatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcat    >  2:4643771‑M1/1‑14 (MQ=255)
                                                                  gatTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcat    <  2:2880016‑M1/76‑63 (MQ=255)
                                                                   atTACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc   >  1:883872‑M1/1‑13 (MQ=255)
                                                                    ttACATCGCTG‑‑Cgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatct  >  2:1368029‑M1/1‑12 (MQ=255)
                                                                        |                                                                         
ACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTG‑‑CATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTT  >  NC_000913/3754238‑3754380

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