Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,609,216 A→G intergenic (‑194/+13) yneE ← / ← uxaB bestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,609,2160AG98.2% 194.9 / ‑4.2 54intergenic (‑194/+13)yneE/uxaBbestrophin family putative inner membrane protein/altronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (37/16);  total (37/17)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGCCAAAAC  >  NC_000913/1609071‑1609357
                                                                                                                                                 |                                                                                                                                             
cacaTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTaa                                                                                                                                           >  1:196894/1‑150 (MQ=255)
 acaTTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACg                                                                                                                                         <  2:754021/151‑1 (MQ=255)
       tcAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGtt                                                                                                                                    >  2:446680/1‑150 (MQ=255)
        cAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGtt                                                                                                                                    <  2:749528/149‑1 (MQ=255)
         aaGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGttatt                                                                                                                                 >  2:135994/1‑151 (MQ=255)
                 aCTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACaa                                                                                                                          <  1:454702/150‑1 (MQ=255)
                     aCAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGAc                                                                                                                     >  2:207335/1‑151 (MQ=255)
                      cAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGAc                                                                                                                     <  1:723957/150‑1 (MQ=255)
                       aaaTTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGt                                                                                                                   >  2:93323/1‑151 (MQ=255)
                                  aaGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTcgcg                                                                                                        >  1:299140/1‑151 (MQ=255)
                                     cTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCAt                                                                                                     >  1:243311/1‑151 (MQ=255)
                                          aacaTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCtt                                                                                                >  2:364310/1‑151 (MQ=255)
                                            caTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCtttt                                                                                              >  2:737457/1‑151 (MQ=255)
                                             aTTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATGAGCACAACGGACGTACAGCTTCGCGCATCCCttttt                                                                                             <  1:29852/151‑1 (MQ=255)
                                                            aTGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCa                                                                              <  1:5869/151‑1 (MQ=255)
                                                                   aCCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCAtt                                                                       >  1:470152/1‑151 (MQ=255)
                                                                        ttAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCa                                                                   >  2:7152/1‑150 (MQ=255)
                                                                           aTGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCt                                                               >  1:287291/1‑151 (MQ=255)
                                                                               tAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTc                                                           >  2:618925/1‑151 (MQ=255)
                                                                                aaaGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCg                                                          >  1:367176/1‑151 (MQ=255)
                                                                                 aaGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCga                                                         <  2:243335/151‑1 (MQ=255)
                                                                                 aaGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCg                                                          >  2:453300/1‑150 (MQ=255)
                                                                                  aGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgac                                                        <  2:549301/151‑1 (MQ=255)
                                                                                      atatAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccaga                                                    <  1:742499/151‑1 (MQ=255)
                                                                                        ataGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacc                                                  <  1:738386/151‑1 (MQ=255)
                                                                                         taGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacc                                                  >  2:451116/1‑150 (MQ=255)
                                                                                          aGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCgaccagacc                                                  >  2:120167/1‑149 (MQ=255)
                                                                                            gATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGc                                              >  1:350562/1‑151 (MQ=255)
                                                                                            gATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGc                                              >  2:631213/1‑151 (MQ=255)
                                                                                             aTGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGc                                              >  2:516180/1‑150 (MQ=255)
                                                                                              tgttgtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCAc                                            >  2:809120/1‑151 (MQ=255)
                                                                                               gttgtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCAc                                            <  2:645133/150‑1 (MQ=255)
                                                                                                  gtAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGc                                        >  1:81375/1‑151 (MQ=255)
                                                                                                    aaCTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGt                                      >  2:191366/1‑151 (MQ=255)
                                                                                                    aaCTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGt                                      <  2:47865/151‑1 (MQ=255)
                                                                                                     aCTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc                                     >  2:528324/1‑151 (MQ=255)
                                                                                                     aCTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc                                     >  2:240079/1‑151 (MQ=255)
                                                                                                      cTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTc                                     >  1:469805/1‑150 (MQ=255)
                                                                                                          tatGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTc                                >  2:257501/1‑151 (MQ=255)
                                                                                                                gAAATAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTc                          >  2:510188/1‑151 (MQ=255)
                                                                                                                    tAAGAAACCCGGTCGAAACCGTGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTccc                        >  1:292588/1‑149 (MQ=255)
                                                                                                                    tAAGAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTccc                        >  1:364478/1‑149 (MQ=255)
                                                                                                                       gAAACCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTgg                   >  2:416366/1‑151 (MQ=255)
                                                                                                                         aaCCCGATCGACACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGGACAGCTTAGCGCATCCCTTTATCGAGAATCGCATCCAGGCCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGCGGt                  >  1:612459/1‑150 (MQ=255)
                                                                                                                          aCCCGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTca                <  1:510692/151‑2 (MQ=255)
                                                                                                                           cccGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCtt               >  2:738972/1‑151 (MQ=255)
                                                                                                                           cccGGTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCtt               >  1:89176/1‑151 (MQ=255)
                                                                                                                              ggTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCttttt            <  2:595225/151‑1 (MQ=255)
                                                                                                                               gTCGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTc           <  1:135971/151‑1 (MQ=255)
                                                                                                                                 cGAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCt          <  1:564420/150‑1 (MQ=255)
                                                                                                                                  gAAACCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGc        >  1:125212/1‑151 (MQ=255)
                                                                                                                                    aaCCGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGcc       >  2:854451/1‑150 (MQ=255)
                                                                                                                                       cGGGTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGCCaaa    >  2:145889/1‑150 (MQ=255)
                                                                                                                                        gggTTCAGAGGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGGCAGGTCCTGCTCCCAGTGGTCTTTTTCTGCCAAAAc  <  2:542775/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCATATAGGATGTTGTAACTAATATGGTGAAATAAGAAACCCGGTCGAAACCGGGTTCAGAAGTAACGGTGTTATTAGCACAACGGACGTACAGCTTCGCGCATCCCTTTTTCGAGAATCGCATCCAGGTCATTAGCAACCTGCTCGACCAGACCAGGCACTTGCGTCAGGTCCTGCTCCCAGTGGTCTTTTTCTGCCAAAAC  >  NC_000913/1609071‑1609357

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