Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,225 T→A 9.3% intergenic (+20/+32) yfeD → / ← gltX DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2250TA 9.3% 179.0 / 8.9 54intergenic (+20/+32)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base A (3/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGC  >  NC_000913/2519079‑2519370
                                                                                                                                                  |                                                                                                                                                 
gTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                                                                 <  1:649331/149‑1 (MQ=255)
   tCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtctc                                                                                                                                              <  2:717285/149‑1 (MQ=255)
             atcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                                                                          >  2:39675/1‑143 (MQ=255)
             atcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                                                                          <  1:39675/143‑1 (MQ=255)
              tcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCt                                                                                                                                            >  3:299685/1‑140 (MQ=255)
              tcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCt                                                                                                                                            <  4:299685/140‑1 (MQ=255)
                         gataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                                                                 <  3:579502/124‑1 (MQ=255)
                         gataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                                                                 >  4:579502/1‑124 (MQ=255)
                             atatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                                                                    >  3:953464/1‑149 (MQ=255)
                                   cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                                                                       <  1:1433695/124‑1 (MQ=255)
                                   cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                                                                       >  2:1433695/1‑124 (MQ=255)
                                       caTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGa                                                                                                          <  3:1319782/149‑1 (MQ=255)
                                        aTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGcc                                                                                                                                         >  3:709021/1‑117 (MQ=255)
                                        aTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGcc                                                                                                                                         <  4:709021/117‑1 (MQ=255)
                                                 gAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                                                                       >  4:1115173/1‑110 (MQ=255)
                                                 gAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                                                                       <  3:1115173/110‑1 (MQ=255)
                                                      aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                                                                           >  4:1095795/1‑149 (MQ=255)
                                                          gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                                                                       >  4:284018/1‑149 (MQ=255)
                                                          gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                                                                       >  2:771279/1‑149 (MQ=255)
                                                              aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                                                                   <  1:771279/149‑1 (MQ=255)
                                                              aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                                                                   <  3:284018/149‑1 (MQ=255)
                                                                tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                                                                 <  1:1301801/85‑1 (MQ=255)
                                                                tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                                                                 >  2:1301801/1‑85 (MQ=255)
                                                                         cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                                                                        <  4:953464/149‑1 (MQ=255)
                                                                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                                             >  3:446154/1‑110 (MQ=255)
                                                                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                                             <  4:446154/110‑1 (MQ=255)
                                                                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                                                                      >  1:323553/1‑149 (MQ=255)
                                                                             aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                    >  1:186093/1‑149 (MQ=255)
                                                                                gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                                                                              <  4:1328524/88‑1 (MQ=255)
                                                                                gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                                                                              >  3:1328524/1‑88 (MQ=255)
                                                                                gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                                                                          >  2:688158/1‑92 (MQ=255)
                                                                                gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                                                                          <  1:688158/92‑1 (MQ=255)
                                                                                          tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                                                              >  2:1396621/1‑110 (MQ=255)
                                                                                          tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                                                              <  1:1396621/110‑1 (MQ=255)
                                                                                             ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                    <  3:1095795/149‑1 (MQ=255)
                                                                                                   ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTTATTACTGCTGATTTTCGCGTTCAGCAAAAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt                                                <  4:339906/147‑1 (MQ=255)
                                                                                                   ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt                                                >  3:339906/1‑147 (MQ=255)
                                                                                                     gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                                            >  4:1276804/1‑149 (MQ=255)
                                                                                                        tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa                                              <  3:999050/144‑1 (MQ=255)
                                                                                                        tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa                                              >  4:999050/1‑144 (MQ=255)
                                                                                                          gATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                                                                 >  2:566764/1‑107 (MQ=255)
                                                                                                          gATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                                                                 <  1:566764/107‑1 (MQ=255)
                                                                                                                gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGaa                                 <  2:323553/149‑1 (MQ=255)
                                                                                                                    aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGt                             <  3:1276804/149‑1 (MQ=255)
                                                                                                                         gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                    >  3:93104/1‑121 (MQ=255)
                                                                                                                         gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                    <  4:93104/121‑1 (MQ=255)
                                                                                                                             aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                                                                       <  1:589904/66‑1 (MQ=255)
                                                                                                                             aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                                                                       >  2:589904/1‑66 (MQ=255)
                                                                                                                                taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGc                 >  3:477082/1‑149 (MQ=255)
                                                                                                                                        ggCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                    >  1:202053/1‑106 (MQ=255)
                                                                                                                                        ggCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                    <  2:202053/106‑1 (MQ=255)
                                                                                                                                        ggCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACt         >  3:111809/1‑149 (MQ=255)
                                                                                                                                         gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGc                                      <  3:61568/119‑1 (MQ=255)
                                                                                                                                         gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGc                                      >  4:61568/1‑119 (MQ=255)
                                                                                                                                             gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                                                                            <  4:577179/77‑1 (MQ=255)
                                                                                                                                             gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                                                                            >  3:577179/1‑77 (MQ=255)
                                                                                                                                             gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                                       <  2:1052110/82‑1 (MQ=255)
                                                                                                                                             gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                                       >  1:1052110/1‑82 (MQ=255)
                                                                                                                                             gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGcccc    <  2:186093/149‑1 (MQ=255)
                                                                                                                                                attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGc  <  2:111610/148‑1 (MQ=255)
                                                                                                                                                attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGc  >  1:111610/1‑148 (MQ=255)
                                                                                                                                                 ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAc                                >  4:1395486/1‑117 (MQ=255)
                                                                                                                                                 ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAc                                <  3:1395486/117‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGC  >  NC_000913/2519079‑2519370

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