Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,510,121 +G intergenic (‑358/+151) yhfZ ← / ← trpS hypothetical protein/tryptophan tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,510,1201.G96.4% 211.6 / 2.7 84intergenic (‑357/+152)yhfZ/trpShypothetical protein/tryptophan tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (51/29);  total (51/33)
Fisher's exact test for biased strand distribution p-value = 5.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAAT  >  NC_000913/3510046‑3510190
                                                                            |                                                                      
gagaTGTGGAAACGGCGAGGCACTTCACAGTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGCTTCCCCg                                                                       >  1:1665620/1‑76 (MQ=255)
 agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCgg                                                                      >  2:2570013/1‑75 (MQ=255)
 agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCgg                                                                      >  1:3053577/1‑75 (MQ=255)
 agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCgg                                                                      <  2:4034266/76‑2 (MQ=255)
   aTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGGGGATAAAAAAACCGGGTTTCCCCggtt                                                                    <  2:3623376/76‑4 (MQ=255)
   aTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCggtt                                                                    >  1:5200779/1‑73 (MQ=255)
   aTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCggtt                                                                    <  1:2747591/76‑4 (MQ=255)
   aTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCggtt                                                                    >  1:5174696/1‑73 (MQ=255)
       ggAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCag                                                                 >  1:4600316/1‑76 (MQ=255)
        gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCaga                                                                >  2:1540152/1‑76 (MQ=255)
          aaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGt                                                              <  1:3045396/76‑1 (MQ=255)
          aaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGt                                                              <  1:495169/76‑1 (MQ=255)
           aCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtg                                                             <  1:1051224/76‑1 (MQ=255)
            cGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtga                                                            >  1:1232292/1‑76 (MQ=255)
             ggCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgat                                                           <  2:1094379/76‑1 (MQ=255)
              gCGAGGCAGTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatc                                                          >  2:2049893/1‑75 (MQ=255)
              gCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatg                                                          >  2:3948940/1‑76 (MQ=255)
                 aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATa                                                       >  2:5264717/1‑76 (MQ=255)
                   gCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaa                                                     <  2:460859/76‑1 (MQ=255)
                    cACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaaa                                                    >  1:4921836/1‑76 (MQ=255)
                    cACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaaa                                                    >  2:1930310/1‑76 (MQ=255)
                     aCTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAg                                                   >  1:2224093/1‑76 (MQ=255)
                      cTTCACACTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGc                                                  <  1:3956157/76‑1 (MQ=255)
                       ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCa                                                 >  1:3636900/1‑76 (MQ=255)
                       ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCa                                                 <  1:5771408/76‑1 (MQ=255)
                        tCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaa                                                <  2:2265150/76‑1 (MQ=255)
                        tCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaa                                                >  1:1099175/1‑76 (MQ=255)
                          acaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaaa                                              >  2:2633182/1‑76 (MQ=255)
                          acaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaaa                                              <  1:267050/76‑1 (MQ=255)
                           caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAt                                             >  1:3473021/1‑76 (MQ=255)
                            aTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAtt                                            <  1:3837559/76‑1 (MQ=255)
                             ttttttCGGATTATTGCAAATGAGGATAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGc                                          >  2:1643095/1‑76 (MQ=255)
                              tttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGc                                          >  1:4790522/1‑76 (MQ=255)
                               ttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGcc                                         <  1:1189200/76‑1 (MQ=255)
                               ttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGGTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGcc                                         >  1:5777848/1‑76 (MQ=255)
                                 ttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTg                                       <  1:343873/76‑1 (MQ=255)
                                  tCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                      >  2:2265941/1‑76 (MQ=255)
                                  tCGGACTATTGAAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCACAGTGATGATAAAAGCAAAATTGCCTGa                                      >  2:4011378/1‑76 (MQ=255)
                                   cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                     >  1:3419179/1‑76 (MQ=255)
                                   cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                     <  2:4243257/76‑1 (MQ=255)
                                   cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                     >  2:4736048/1‑76 (MQ=255)
                                     gATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgc                                   >  1:3906639/1‑76 (MQ=255)
                                      attattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcg                                  >  2:286072/1‑76 (MQ=255)
                                       ttattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                 >  1:4614971/1‑76 (MQ=255)
                                         attGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTa                               >  2:3258614/1‑76 (MQ=255)
                                          ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                              >  2:2672518/1‑76 (MQ=255)
                                          ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                              >  2:639628/1‑76 (MQ=255)
                                          ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                              >  2:5175001/1‑76 (MQ=255)
                                          ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                              <  2:1218413/76‑1 (MQ=255)
                                           tGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACg                             >  1:5713366/1‑76 (MQ=255)
                                             cAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCt                           <  1:3920675/76‑1 (MQ=255)
                                             cAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCt                           <  1:5417816/76‑1 (MQ=255)
                                              aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCtt                          >  2:2931291/1‑76 (MQ=255)
                                               aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa                         >  2:1360427/1‑76 (MQ=255)
                                               aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa                         <  1:2530116/76‑1 (MQ=255)
                                                aTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTAt                        >  1:556783/1‑76 (MQ=255)
                                                 tGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATc                       <  2:1738320/76‑1 (MQ=255)
                                                 tGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATc                       >  2:5190899/1‑76 (MQ=255)
                                                 tGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATc                       >  2:801782/1‑76 (MQ=255)
                                                  gAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCa                      <  1:3290226/76‑1 (MQ=255)
                                                     gATAAAAAAACCGGGTCTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAGGCGCTACGCTTATCAGGc                   >  2:2833366/1‑76 (MQ=25)
                                                      aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGcc                  <  2:1665620/76‑1 (MQ=255)
                                                      aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGcc                  >  1:4958650/1‑76 (MQ=255)
                                                         aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc               >  2:1779282/1‑76 (MQ=255)
                                                         aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc               >  1:4542799/1‑76 (MQ=255)
                                                          aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa              >  2:941255/1‑76 (MQ=255)
                                                           aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt             <  2:5098440/76‑1 (MQ=255)
                                                           aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt             >  1:2792490/1‑76 (MQ=255)
                                                            aagccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            <  1:1841111/73‑1 (MQ=21)
                                                            aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            >  1:1264330/1‑76 (MQ=255)
                                                            aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            <  1:1813668/76‑1 (MQ=255)
                                                            aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            <  2:1374431/76‑1 (MQ=255)
                                                            aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            >  2:3164417/1‑76 (MQ=255)
                                                            aaaCCGGGTGTCACCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            <  1:1841244/76‑1 (MQ=18)
                                                            aaaCCGGGTGTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            <  1:1841214/76‑1 (MQ=21)
                                                            aaaCCGGGTGTACCAGGTGTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt            <  1:1841221/76‑1 (MQ=17)
                                                             aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAttt           >  2:5292414/1‑76 (MQ=255)
                                                              aCCTGGTTTCCTCTGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTCTCAGGCCTCCATTTc          <  1:599592/76‑1 (MQ=18)
                                                              aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTc          >  2:163676/1‑76 (MQ=255)
                                                              aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTc          >  1:2698028/1‑76 (MQ=255)
                                                              aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTc          <  2:3359827/76‑1 (MQ=255)
                                                               ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTcc         >  1:3759670/1‑76 (MQ=255)
                                                                 gggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCtt       >  1:3427141/1‑76 (MQ=255)
                                                                 gggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCtt       <  1:5190900/76‑1 (MQ=255)
                                                                  ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTg      >  2:1716162/1‑76 (MQ=255)
                                                                  ggTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGTCCTTGc     <  2:5200778/76‑1 (MQ=18)
                                                                   gTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGc     <  2:5174695/76‑1 (MQ=255)
                                                                   gTTCCCCGGTTTCAGAGTGATTATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTTCCTTGCa    <  2:3532045/76‑1 (MQ=18)
                                                                    tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCa    >  1:3058859/1‑76 (MQ=255)
                                                                    tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCa    >  2:1448702/1‑76 (MQ=255)
                                                                     ttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCaa   >  1:676559/1‑76 (MQ=255)
                                                                      tCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAat  >  2:3720303/1‑76 (MQ=255)
                                                                            |                                                                      
GAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAAT  >  NC_000913/3510046‑3510190

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