Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,451,536 A→G 9.1% intergenic (‑187/+61) tynA ← / ← paaZ tyramine oxidase, copper‑requiring/oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,451,5360AG 9.1% 210.9 / 9.3 66intergenic (‑187/+61)tynA/paaZtyramine oxidase, copper‑requiring/oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (3/3);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGA  >  NC_000913/1451414‑1451676
                                                                                                                          |                                                                                                                                            
tCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaaa                                                                                                                     >  3:608325/1‑148 (MQ=255)
      tATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacat                                                                                                                                     >  3:1043198/1‑125 (MQ=255)
      tATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacat                                                                                                                                     <  4:1043198/126‑2 (MQ=255)
      tATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAg                                                                                                                              <  2:690878/133‑1 (MQ=255)
      tATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAg                                                                                                                              >  1:690878/1‑133 (MQ=255)
           gTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacat                                                                                                                                     <  1:765671/121‑2 (MQ=255)
             aCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCGTAATGTGACAACGTCACAAAACtt                                                                                                                                >  3:943161/1‑124 (MQ=255)
             aCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                               <  2:989410/125‑1 (MQ=255)
             aCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                               >  1:989410/1‑125 (MQ=255)
             aCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATCTACACATCATGCATAATGTGACAACGTCACAAAACtt                                                                                                                                <  4:943161/124‑1 (MQ=255)
                 aTGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaa                                                                                                                                    >  1:416680/1‑116 (MQ=255)
                 aTGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaa                                                                                                                                    <  2:974267/116‑1 (MQ=255)
                 aTGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaa                                                                                                                                    <  2:416680/116‑1 (MQ=255)
                 aTGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaa                                                                                                                                    >  1:974267/1‑116 (MQ=255)
                    cAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTg                                                                                                 <  4:608325/148‑1 (MQ=255)
                      aaGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGaa                                                                                                                          >  1:534626/1‑121 (MQ=255)
                      aaGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGaa                                                                                                                          <  2:534626/121‑1 (MQ=255)
                         gggATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcaca                                                                                                                                      >  4:878760/1‑106 (MQ=255)
                         gggATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcaca                                                                                                                                      <  3:878760/106‑1 (MQ=255)
                                  tGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCa                                                                                                               <  3:631248/120‑1 (MQ=255)
                                  tGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCa                                                                                                               >  4:631248/1‑120 (MQ=255)
                                   gCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacatta                                                                                                                                   >  3:399552/1‑96 (MQ=255)
                                        tGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAAt                                                                             <  1:320171/148‑1 (MQ=255)
                                        tGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAAt                                                                             <  3:248217/148‑1 (MQ=255)
                                             ggTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAgg                                                                                                                  <  3:1033567/106‑1 (MQ=255)
                                             ggTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAgg                                                                                                                  >  4:1033567/1‑106 (MQ=255)
                                                  tttGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTAtt                                                                                                         <  1:1054356/110‑1 (MQ=255)
                                                  tttGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTAtt                                                                                                         >  2:1054356/1‑110 (MQ=255)
                                                        aGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGaa                                                                                                                          <  4:796433/87‑1 (MQ=255)
                                                        aGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGaa                                                                                                                          >  3:796433/1‑87 (MQ=255)
                                                                cAGCTTTTCATACGGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAg                                                                                                                   >  2:450089/1‑86 (MQ=255)
                                                                cAGCTTTTCATACGGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAg                                                                                                                   <  1:450089/86‑1 (MQ=255)
                                                                cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCa                                                                  >  2:672005/1‑135 (MQ=255)
                                                                cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCa                                                                  <  1:672005/135‑1 (MQ=255)
                                                                 aGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                               <  3:939825/73‑1 (MQ=255)
                                                                 aGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                               >  4:939825/1‑73 (MQ=255)
                                                                 aGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCAcc                                                                >  1:324319/1‑136 (MQ=255)
                                                                 aGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCAcc                                                                <  2:324319/136‑1 (MQ=255)
                                                                                ggATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAAt                                     >  1:231244/1‑148 (MQ=255)
                                                                                   tGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTAtt                                                                                                         <  2:7503/77‑1 (MQ=255)
                                                                                   tGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTAtt                                                                                                         >  1:7503/1‑77 (MQ=255)
                                                                                    gCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCttc                                                                                           >  3:847038/1‑90 (MQ=255)
                                                                                    gCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCttc                                                                                           <  4:847038/90‑1 (MQ=255)
                                                                                         gAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAg                            >  2:873072/1‑148 (MQ=255)
                                                                                              tAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAg                                                                                                                              <  3:472889/45‑1 (MQ=255)
                                                                                              tAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAg                                                                                                                              >  4:472889/1‑45 (MQ=255)
                                                                                              tAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCa                                                                  <  3:1098701/105‑1 (MQ=255)
                                                                                              tAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCa                                                                  >  4:1098701/1‑105 (MQ=255)
                                                                                                gATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTaa                                                                              >  1:165176/1‑91 (MQ=255)
                                                                                                gATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTaa                                                                              <  3:863625/91‑1 (MQ=255)
                                                                                                gATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTaa                                                                              >  4:863625/1‑91 (MQ=255)
                                                                                                gATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTaa                                                                              <  2:165176/91‑1 (MQ=255)
                                                                                                   gTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCggg                  <  2:231244/148‑1 (MQ=255)
                                                                                                      cacaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACaaa                                                                      >  1:9445/1‑93 (MQ=255)
                                                                                                      cacaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACaaa                                                                      <  2:9445/93‑1 (MQ=255)
                                                                                                      cacaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGcc                                                        >  2:77301/1‑107 (MQ=255)
                                                                                                      cacaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGcc                                                        <  1:77301/107‑1 (MQ=255)
                                                                                                       acaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATc                                                                   <  1:105855/95‑1 (MQ=255)
                                                                                                       acaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATc                                                                   >  2:105855/1‑95 (MQ=255)
                                                                                                           catGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATg          <  1:873072/148‑1 (MQ=255)
                                                                                                             tGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTaa                                                                              >  4:567450/1‑78 (MQ=255)
                                                                                                             tGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTaa                                                                              <  3:567450/78‑1 (MQ=255)
                                                                                                               cATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGa      >  4:928461/1‑148 (MQ=255)
                                                                                                                   aTGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGa  <  3:928461/148‑1 (MQ=255)
                                                                                                                      tgACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGa            >  3:961479/1‑135 (MQ=255)
                                                                                                                      tgACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGa            <  4:961479/135‑1 (MQ=255)
                                                                                                                          |                                                                                                                                            
TCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGA  >  NC_000913/1451414‑1451676

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