Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,443,919 T→A D563V (GAT→GTT)  tynA ← primary‑amine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,443,9190TA96.9% 115.0 / ‑3.4 32D563V (GAT→GTT) tynAprimary‑amine oxidase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (16/15);  total (17/15)
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.82e-01

GGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGC  >  NZ_CP009273/1443837‑1444004
                                                                                  |                                                                                     
ggTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCg                                                                                >  2:608119/1‑90 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGATCGAGGCGGAAAt                                                                           >  2:924681/1‑90 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGa                                                                              <  2:512852/87‑1 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGa                                                                              >  1:512852/1‑87 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAAt                                                                           >  1:272651/1‑90 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAAt                                                                           >  1:367897/1‑90 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAAt                                                                           >  1:391657/1‑90 (MQ=255)
     gTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAAt                                                                           >  2:326823/1‑90 (MQ=255)
                gCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGt                                                                <  1:326823/90‑1 (MQ=255)
                gCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGt                                                                <  1:618859/90‑1 (MQ=255)
                 cAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGtt                                                               <  1:255899/90‑1 (MQ=255)
                          gTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTgtag                                                      <  2:1000315/90‑1 (MQ=255)
                           ttttACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTgtagt                                                     <  1:641929/90‑1 (MQ=255)
                             ttACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTgtag                                                      <  1:752190/87‑1 (MQ=255)
                             ttACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTgtag                                                      >  2:752190/1‑87 (MQ=255)
                             ttACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTAc                                                   >  2:623949/1‑90 (MQ=255)
                              tACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTAcc                                                  <  1:623949/90‑1 (MQ=255)
                               aCCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATaaa                                                                     <  1:682547/70‑1 (MQ=255)
                               aCCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATaaa                                                                     >  2:682547/1‑70 (MQ=255)
                                    tGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTAcc                                                  <  2:277504/84‑1 (MQ=255)
                                    tGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTAcc                                                  >  1:277504/1‑84 (MQ=255)
                                    tGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGat                                            >  2:765145/1‑90 (MQ=255)
                                           aTCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGa                                     <  2:399035/90‑1 (MQ=255)
                                                 accaGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATa                               >  1:963630/1‑90 (MQ=255)
                                                 accaGGCTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATa                               <  2:939692/90‑1 (MQ=255)
                                                   caGGCTGTTATCCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAg                             <  1:596364/90‑1 (MQ=255)
                                                       cTGTTATTCTCGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTg                         >  1:375772/1‑90 (MQ=255)
                                                               ctcGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAgcg                 >  2:961801/1‑90 (MQ=255)
                                                                 cGCCATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGcgcg               >  2:799528/1‑90 (MQ=255)
                                                                    cATCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCgtgt            <  1:799528/90‑1 (MQ=255)
                                                                      tCTACATCCAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCa          >  1:904466/1‑90 (MQ=255)
                                                                              cAGAACGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGc  <  1:912056/90‑1 (MQ=255)
                                                                                  |                                                                                     
GGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACTGGGTCCATCGCCACCAGGCTGTTATTCTCGCCATCTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGC  >  NZ_CP009273/1443837‑1444004

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