Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,331 C→T 9.4% intergenic (+42/‑80) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3310CT 9.4% 159.9 / 6.7 54intergenic (+42/‑80)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (30/18);  new base T (2/3);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGG  >  NC_000913/4612194‑4612478
                                                                                                                                         |                                                                                                                                                    
cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCtt                                                                                                                                           <  2:814027/148‑1 (MQ=255)
cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCtt                                                                                                                                           <  2:853067/148‑1 (MQ=255)
   gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTg                                                                                                                                        <  1:436426/148‑1 (MQ=255)
      ggTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGt                                                                                                                                    <  3:562372/149‑1 (MQ=255)
             cTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCt                                                                                                                                            >  1:867384/1‑134 (MQ=255)
             cTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCt                                                                                                                                            <  2:867384/134‑1 (MQ=255)
                          tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                                                                                 >  1:613929/1‑148 (MQ=255)
                             aaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGt                                                                                                              >  1:335373/1‑148 (MQ=255)
                               ttGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTAc                                                                                                           >  4:510209/1‑149 (MQ=255)
                                     ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                      <  1:26866/148‑1 (MQ=255)
                                     ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                      >  2:40732/1‑148 (MQ=255)
                                           tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCt                                                                                                >  2:252309/1‑148 (MQ=255)
                                              ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                             >  1:699688/1‑148 (MQ=255)
                                                attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                          >  2:217569/1‑101 (MQ=255)
                                                attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                                                                          >  3:696910/1‑149 (MQ=255)
                                                attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCttt                                                                                                                                          <  1:217569/101‑1 (MQ=255)
                                                      cTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                         >  1:303497/1‑128 (MQ=255)
                                                      cTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                         <  2:303497/128‑1 (MQ=255)
                                                      cTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaa                                                                                     >  2:857633/1‑148 (MQ=255)
                                                         ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc                                                                                  <  1:641893/148‑1 (MQ=255)
                                                               gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCTAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                           >  3:263749/1‑149 (MQ=255)
                                                               gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                            <  2:144729/148‑1 (MQ=255)
                                                                tGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                   >  2:614478/1‑140 (MQ=255)
                                                                tGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                   <  1:614478/140‑1 (MQ=255)
                                                                    tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                      >  4:706814/1‑149 (MQ=255)
                                                                    tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                       <  2:335373/148‑1 (MQ=255)
                                                                     ttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                      >  2:788052/1‑148 (MQ=255)
                                                                     ttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGA‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGAAAAGCagag                                                                      >  2:50001/1‑148 (MQ=255)
                                                                           aTGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGc                                                               >  4:647233/1‑149 (MQ=255)
                                                                               gCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTTCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                            >  2:809778/1‑148 (MQ=255)
                                                                               gCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                            <  1:809778/148‑1 (MQ=255)
                                                                                 cGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGa                                                          <  1:857633/148‑1 (MQ=255)
                                                                                            cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAAc                                               >  2:358300/1‑148 (MQ=255)
                                                                                            cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                              >  4:453824/1‑149 (MQ=255)
                                                                                             tAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGc                                             >  4:389977/1‑149 (MQ=255)
                                                                                                aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                          <  3:510209/149‑1 (MQ=255)
                                                                                                aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                          <  3:453824/149‑1 (MQ=255)
                                                                                                 tttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTAGGGCAGCGAATACCTGTAACGCTTg                                          <  2:613929/148‑1 (MQ=255)
                                                                                                  ttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTgg                                         >  1:707115/1‑148 (MQ=255)
                                                                                                         cATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                                 >  4:176931/1‑149 (MQ=255)
                                                                                                                aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  <  2:872079/109‑1 (MQ=255)
                                                                                                                aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  >  1:872079/1‑109 (MQ=255)
                                                                                                                  aaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGt                        >  4:740061/1‑149 (MQ=255)
                                                                                                                  aaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGt                        <  3:176931/149‑1 (MQ=255)
                                                                                                                           cAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  >  4:196919/1‑98 (MQ=255)
                                                                                                                           cAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                                  <  3:196919/98‑1 (MQ=255)
                                                                                                                            aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  3:304427/1‑149 (MQ=255)
                                                                                                                            aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGc               >  1:337637/1‑148 (MQ=255)
                                                                                                                             aaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  1:588379/1‑148 (MQ=255)
                                                                                                                             aaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  1:33544/1‑148 (MQ=255)
                                                                                                                             aaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  2:623032/1‑148 (MQ=255)
                                                                                                                                      cAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTTCAGcc     <  1:40732/148‑1 (MQ=255)
                                                                                                                                       aGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg   >  3:7999/1‑149 (MQ=255)
                                                                                                                                        gC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTgg  >  3:479073/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
CTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGG  >  NC_000913/4612194‑4612478

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