Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,450,499 C→G 9.3% G284A (GGC→GCC)  tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,450,4990CG 9.3% 132.6 / 6.5 43G284A (GGC→GCC) tynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base G (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

ACGCGGTCACGGCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGG  >  NC_000913/1450364‑1450643
                                                                                                                                       |                                                                                                                                                
aCGCGGTCACGGCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATc                                                                                                                                     <  5:285224/149‑1 (MQ=255)
      tCACGGCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTATTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACAGACATCAAGATa                                                                                                                               <  8:348147/149‑1 (MQ=255)
         cGGCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCt                                                                                                                            <  4:223129/149‑1 (MQ=255)
         cGGCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCt                                                                                                                            <  2:29580/149‑1 (MQ=255)
          ggCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCt                                                                                                                            >  7:208380/1‑148 (MQ=255)
            ccATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCtgat                                                                                                                         <  8:40222/149‑1 (MQ=255)
            ccATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCtgat                                                                                                                         >  2:279403/1‑149 (MQ=255)
             cATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGcc                                                                                                                                                 >  2:8324/1‑124 (MQ=255)
             cATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGcc                                                                                                                                                 <  1:8324/124‑1 (MQ=255)
                 aaaTGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACtt                                                                                                                    <  4:52381/149‑1 (MQ=255)
                  aaTGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACttt                                                                                                                   >  2:349122/1‑149 (MQ=255)
                      ggCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGc                                                                                                               <  4:84164/149‑1 (MQ=255)
                            cGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAAAGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCgg                                                                                                         <  8:67745/149‑1 (MQ=255)
                               tCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGAACGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCa                                                                                                      >  8:198684/1‑149 (MQ=255)
                                     gCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGGCATCACCGACATCAAGATAg                                                                                                                              <  2:371232/119‑1 (MQ=255)
                                     gCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGGCATCACCGACATCAAGATAg                                                                                                                              <  2:376566/119‑1 (MQ=255)
                                     gCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGGCATCACCGACATCAAGATAg                                                                                                                              >  1:376566/1‑119 (MQ=255)
                                     gCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGGCATCACCGACATCAAGATAg                                                                                                                              >  1:371232/1‑119 (MQ=255)
                                     gCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTg                                                                                                 <  4:107537/148‑1 (MQ=255)
                                        ccGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTc                                                                                             >  6:4844/1‑149 (MQ=255)
                                         cGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGttt                                                                                              >  2:110346/1‑147 (MQ=255)
                                         cGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCa                                                                                            <  1:349122/149‑1 (MQ=255)
                                         cGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCa                                                                                            <  1:110346/149‑1 (MQ=255)
                                          ggAACTACCGGACCTTCTTCAATCTAAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGATTTCGATGCGATGTGCCCAGTAGTTGCCATCACCGCCATCAAGATAGCTGATGACTTTGAGCAACCGGCCATCTAGTTTCAg                                                                                           >  7:154733/1‑149 (MQ=255)
                                                 ccGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCAt                                                                                      <  1:296322/147‑1 (MQ=255)
                                                              aaTCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGaaa                                                                        >  8:343167/1‑148 (MQ=255)
                                                                  ttAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAAcc                                                                   >  1:87104/1‑149 (MQ=255)
                                                                  ttAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAAcc                                                                   >  8:142163/1‑149 (MQ=255)
                                                                    aaCGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTa                                                                 >  4:325119/1‑149 (MQ=255)
                                                                    aaCGATTGTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAAc                                                                    <  2:209933/146‑1 (MQ=255)
                                                                      cGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACg                                                               <  3:246526/149‑1 (MQ=255)
                                                                       gATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACgg                                                              >  5:247625/1‑149 (MQ=255)
                                                                                cTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAgcggcg                                                     >  1:51238/1‑149 (MQ=255)
                                                                                     cTAAATCAACGACCGCCACCAGGTTATCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGt                                                 >  2:203059/1‑148 (MQ=255)
                                                                                      tAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTaa                                               >  8:184641/1‑149 (MQ=255)
                                                                                          tCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCAc                                           <  7:178259/149‑1 (MQ=255)
                                                                                                 ccGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGCCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTc                                    >  2:80860/1‑149 (MQ=255)
                                                                                                  cGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCg                                   <  7:142163/149‑1 (MQ=255)
                                                                                                    ccaccaGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGAAACCGGGCATCTTGTTTCAGGCAATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCCCTTTTTTCGCa                                 >  3:290750/1‑149 (MQ=255)
                                                                                                                 gATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACcgc                    >  7:94842/1‑149 (MQ=255)
                                                                                                                      gATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTttc               >  2:299538/1‑149 (MQ=255)
                                                                                                                             ccAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACgg        <  8:208380/149‑1 (MQ=255)
                                                                                                                              cAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACgg        >  8:174712/1‑148 (MQ=255)
                                                                                                                                   gtTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCgg  >  2:289868/1‑149 (MQ=255)
                                                                                                                                    ttGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCgg  >  1:221885/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                                
ACGCGGTCACGGCCATCAAATGGGCGTGCGGTCATTGGCACCGGAACTACCGGACCTTCTTCAATCTTAACGATTTTTTTCTGTTCTAAATCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGG  >  NC_000913/1450364‑1450643

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