Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,066,248 C→T 8.0% intergenic (‑94/+57) yeeH ← / ← insH1 pseudogene/IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,066,2480CT 8.0% 137.2 / 3.0 50intergenic (‑94/+57)yeeH/insH1pseudogene/IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (26/20);  new base T (2/2);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

AAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCAC  >  NC_000913/2066100‑2066323
                                                                                                                                                    |                                                                           
aaaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                                             >  5:324837/1‑149 (MQ=255)
 aaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                                            >  1:166059/1‑149 (MQ=255)
 aaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                                            >  5:417218/1‑149 (MQ=255)
       cTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                                                      >  6:270048/1‑149 (MQ=255)
        ttGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGAGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                                                      >  3:214499/1‑148 (MQ=255)
          gTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                                                    <  7:33979/148‑1 (MQ=255)
          gTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                                                    <  3:102304/148‑1 (MQ=255)
           tAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                                                  <  5:394231/149‑1 (MQ=255)
             aaaCAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGcc                                                                 >  4:76144/1‑148 (MQ=255)
                caacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                                             <  1:190832/149‑1 (MQ=255)
                caacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGAGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                                             <  4:214499/149‑1 (MQ=255)
                  acGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                                           >  6:205599/1‑149 (MQ=255)
                     tAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                                                        <  5:270048/149‑1 (MQ=255)
                       aTAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                                                       <  2:166059/148‑1 (MQ=255)
                         aGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGtt                                                    >  3:43685/1‑149 (MQ=255)
                           aaTTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                                                  >  5:172510/1‑149 (MQ=255)
                                tCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGc                                             >  6:5932/1‑149 (MQ=255)
                                tCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGc                                             >  6:12502/1‑149 (MQ=255)
                                 ccGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCga                                             >  6:286301/1‑147 (MQ=255)
                                  cGGAATGGTGGCGACTATGCACTAGGGAAGGGGCGAATAAGTGGGGAAAGTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCa                                           >  5:346110/1‑149 (MQ=255)
                                      aTGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATATCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                                        >  2:381148/1‑148 (MQ=255)
                                       tggtggCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                                        >  3:422598/1‑147 (MQ=255)
                                             cGACTATGCAATAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTTAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt                                <  6:346110/149‑1 (MQ=255)
                                                 tATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc                            >  4:158140/1‑149 (MQ=38)
                                                 tATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc                            >  1:310147/1‑149 (MQ=38)
                                                   tGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                                             <  1:288649/114‑1 (MQ=38)
                                                   tGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                                             >  2:288649/1‑114 (MQ=38)
                                                     cACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc                                   <  3:275597/138‑1 (MQ=34)
                                                     cACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc                                   >  4:275597/1‑138 (MQ=34)
                                                     cACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa                         <  3:424125/148‑1 (MQ=34)
                                                     cACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAg                        >  3:136620/1‑149 (MQ=34)
                                                      aCTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGTCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGt                       >  8:442393/1‑149 (MQ=18)
                                                      aCTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGTCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCATTATCCCCAGt                       >  8:448666/1‑149 (MQ=18)
                                                      aCTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGt                       >  4:201756/1‑149 (MQ=33)
                                                        tatggAAGGTGCGAATAAGTGGTGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGCGCCGATTTTTTCTGCCGTAAATGCCTTGAATCAGCCTATTTAGCCCGTGTCTTCGCCATTTAAGGCGTTATCCCCAGttt                     >  6:314183/4‑149 (MQ=1)
                                                        tAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGCCTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGttt                     >  1:88418/1‑149 (MQ=17)
                                                          gggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGtttt                    <  3:384168/148‑1 (MQ=32)
                                                          gggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGtttt                    >  3:241336/1‑148 (MQ=32)
                                                           ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGTCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCTTTATCCCCAGTTTTTa                  <  7:448666/149‑1 (MQ=17)
                                                           ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGTCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTa                  <  7:442393/149‑1 (MQ=17)
                                                           ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTa                  <  5:205599/149‑1 (MQ=32)
                                                           ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTa                  <  1:146528/149‑1 (MQ=32)
                                                           ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTa                  <  5:5932/149‑1 (MQ=32)
                                                           ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTa                  <  5:228014/149‑1 (MQ=32)
                                                             aaGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGt                <  2:310147/149‑1 (MQ=32)
                                                               ggTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTga              <  4:241336/149‑1 (MQ=32)
                                                                 tGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGt                >  8:144790/1‑145 (MQ=32)
                                                                  gCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGt                >  5:396079/1‑144 (MQ=32)
                                                                  gCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTgaga            >  1:66725/1‑148 (MQ=32)
                                                                     aaTAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtctc        >  1:69696/1‑149 (MQ=32)
                                                                       tAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTTAATCAGCCTATTTAGACCGTTGCTTCGCCATTTCAGGCGTTATCCCCAGTTTTTAGTGAGAtctctc      <  1:397280/149‑1 (MQ=12)
                                                                           tGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCAc  <  3:201756/149‑1 (MQ=11)
                                                                           tGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCAc  >  6:119416/1‑149 (MQ=11)
                                                                                                                                                    |                                                                           
AAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCAC  >  NC_000913/2066100‑2066323

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