Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 584,251 T→C Y191H (TAT→CAT)  appY → DNA‑binding transcriptional activator AppY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913584,2510TC96.5% 183.1 / ‑2.2 57Y191H (TAT→CAT) appYDNA‑binding transcriptional activator AppY
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (28/27);  total (29/28)
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.31e-01

TATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGA  >  NC_000913/584124‑584386
                                                                                                                               |                                                                                                                                       
tatTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAAc                                                                                                                             >  2:733708/1‑140 (MQ=255)
    gAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGATATATACTAGGCTGAGGCATGCAAAAAAACTCAt                                                                                                                         >  2:367869/1‑140 (MQ=255)
    gAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCAt                                                                                                                         >  1:140681/1‑140 (MQ=255)
    gAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCAt                                                                                                                         >  2:814425/1‑140 (MQ=255)
    gAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCAt                                                                                                                         >  1:100363/1‑140 (MQ=255)
    gAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCAt                                                                                                                         >  1:141239/1‑140 (MQ=255)
    gAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAAACTCGGGTGAGGCGTGCAAAAAAACTCAt                                                                                                                         >  1:403292/1‑140 (MQ=255)
       cGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGGGAGATACGAGGGTGAGGCATGGAAAAAAACTCATAAc                                                                                                                      >  1:236086/1‑140 (MQ=255)
         gCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACtt                                                                                                                    >  2:408151/1‑140 (MQ=255)
           aaTGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCa                                                                                                                  <  2:718935/140‑1 (MQ=255)
           aaTGGCATTTAAAAGATATTGCAGAATTGATATATACGAGTGAAAGTTAAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCa                                                                                                                  <  2:299339/140‑1 (MQ=255)
            aTGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCaa                                                                                                                 >  2:607229/1‑140 (MQ=255)
                cATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAActc                                                                                                             >  2:150777/1‑140 (MQ=255)
                   ttAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTACGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTa                                                                                                          <  2:499265/140‑1 (MQ=255)
                   ttAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTa                                                                                                          <  2:636805/140‑1 (MQ=255)
                   ttAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTa                                                                                                          <  1:476479/140‑1 (MQ=255)
                    tAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGACAAGCAAAAAAACTCATAACTTCAAACTCTTAt                                                                                                         >  2:398993/1‑140 (MQ=255)
                       aaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCt                                                                                                      <  1:814425/140‑1 (MQ=255)
                       aaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCt                                                                                                      >  1:702477/1‑140 (MQ=255)
                                 ggAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                                                            <  2:702477/140‑1 (MQ=255)
                                 ggAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                                                            <  1:733708/140‑1 (MQ=255)
                                 ggAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                                                            <  1:589680/140‑1 (MQ=255)
                                 ggAATTGATTTATACCAGTGAATGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                                                            <  1:367869/140‑1 (MQ=255)
                                   aaTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTa                                                                                          <  2:141239/140‑1 (MQ=255)
                                     ttGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                                                            <  2:463295/136‑1 (MQ=255)
                                     ttGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                                                            >  1:463295/1‑136 (MQ=255)
                                        aTTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTACCAATGTCGTAGCCCa                                                                                     >  1:309831/1‑140 (MQ=255)
                                             taCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAGCTTGAAAATCTTATTCGATCAGTGTCGTAGCCCAGAAAt                                                                                >  2:317891/1‑140 (MQ=255)
                                             taCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaac                                                                                >  1:86649/1‑139 (MQ=255)
                                                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTTATAACTTCAAACTCTTTTTCTATCACTGTCCGAGCCCAGAAATGTGGc                                                                           >  2:369487/1‑140 (MQ=255)
                                                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                                                           >  2:427349/1‑140 (MQ=255)
                                                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                                                           >  1:200607/1‑140 (MQ=255)
                                                   tGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCt                                                                          >  1:647118/1‑140 (MQ=255)
                                                     aaaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                                                                        <  2:140681/140‑1 (MQ=255)
                                                       aGTTTAATTAAAAAAAGATTAACGGATGAATGATCTTCATTTACTGAATTATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                                                                      <  2:403292/140‑1 (MQ=255)
                                                       aGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                                                                      <  1:468783/140‑1 (MQ=255)
                                                         tttAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAcg                                                                    <  1:408151/140‑2 (MQ=255)
                                                          ttAATAAAAAAAAGATTAAGGGTTGAAGGAACGTCATTTCCTGAAATAATGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAcgg                                                                   <  2:309831/140‑3 (MQ=255)
                                                            aataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTa                                                                 <  1:427349/140‑1 (MQ=255)
                                                                                    aGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTa                                         <  2:647118/140‑1 (MQ=255)
                                                                                    aGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTa                                         <  1:742363/140‑1 (MQ=255)
                                                                                    aGGAACGTCATTTACTGAAATAGTGAGAGATACTAGGGTGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTa                                         <  1:336334/140‑1 (MQ=255)
                                                                                      gAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTaaa                                       <  1:607229/140‑1 (MQ=255)
                                                                                      gAACGTCATTTACTGAAATATTAAGAGATCCTAGGATTAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTaaa                                       <  1:317891/140‑1 (MQ=255)
                                                                                         cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGat                                    >  2:546573/1‑140 (MQ=255)
                                                                                                   tGAAATATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGgtg                          <  1:150777/140‑1 (MQ=255)
                                                                                                     aaaTATTGAGAGATACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTc                        >  1:791629/1‑140 (MQ=255)
                                                                                                        tatTGAGAGATACTAGGATGAGGCAAGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACg                     <  1:372982/140‑1 (MQ=255)
                                                                                                             agagaTACTAGGATGAGACATGTAAAAAAACGCATTACTTCAAACTCTAATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCAtc                <  2:236086/140‑1 (MQ=255)
                                                                                                               agaTACTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCAtctc              >  2:112879/1‑140 (MQ=255)
                                                                                                               agaTACTAGGATGAGACATGCATAAAAACGCATAACTGCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCAtctc              <  1:398993/140‑1 (MQ=255)
                                                                                                                    cTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACg                     >  1:536183/1‑128 (MQ=255)
                                                                                                                    cTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACg                     <  2:536183/128‑1 (MQ=255)
                                                                                                                    cTAGGATGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCattat         >  2:203286/1‑140 (MQ=255)
                                                                                                                          tGAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTgag   >  2:495909/1‑140 (MQ=255)
                                                                                                                           gAGGCATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACGGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTgaga  >  2:411642/1‑140 (MQ=255)
                                                                                                                            aGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCa                  >  2:503809/1‑123 (MQ=255)
                                                                                                                            aGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCa                  <  1:503809/123‑1 (MQ=255)
                                                                                                                               |                                                                                                                                       
TATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTTTGAGA  >  NC_000913/584124‑584386

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