Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,705,165 G→A A463T (GCC→ACC)  nadB → L‑aspartate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,705,1650GA97.6% 152.1 / ‑4.2 42A463T (GCC→ACC) nadBL‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (18/23);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAAT  >  NZ_CP009273/2705032‑2705292
                                                                                                                                     |                                                                                                                               
gTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTg                                                                                                                            <  1:748887/139‑1 (MQ=255)
       ccgTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTgcg                                                                                                                          <  1:805914/134‑1 (MQ=255)
       ccgTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTgcg                                                                                                                          >  2:805914/1‑134 (MQ=255)
       ccgTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGa                                                                                                                     <  2:811478/139‑1 (MQ=255)
       ccgTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGa                                                                                                                     <  2:1159710/139‑1 (MQ=255)
       ccgTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGa                                                                                                                     <  1:286488/139‑1 (MQ=255)
                      cgcgTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCcaac                                                                                                      <  2:492742/139‑1 (MQ=255)
                             agaACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGaaa                                                                                                >  1:783165/1‑138 (MQ=255)
                                    tGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGaa                                                                                        >  2:952256/1‑139 (MQ=255)
                                           cGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACg                                                                                 <  2:372656/139‑1 (MQ=255)
                                            gggTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGc                                                                                <  1:759494/139‑1 (MQ=255)
                                                    aTTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGaaa                                                                                                <  1:418132/115‑1 (MQ=255)
                                                    aTTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGaaa                                                                                                >  2:418132/1‑115 (MQ=255)
                                                        aGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGt                                                                    >  1:1170739/1‑139 (MQ=255)
                                                          cATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGtct                                                                  >  2:368706/1‑139 (MQ=255)
                                                                       aGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGAGCACAACGAAGCGCCTGGGACGAACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGAGTCTCAAATAATTTGCt                                                     >  1:593737/1‑139 (MQ=255)
                                                                            cGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAg                                                 <  2:783165/138‑1 (MQ=255)
                                                                            cGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGc                                                >  1:940353/1‑139 (MQ=255)
                                                                                  tttATGTGGGATTACGTTGGCATTGTGCGCACCACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGt                                           >  1:1115945/1‑138 (MQ=255)
                                                                                  tttATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCg                                            <  1:565048/137‑1 (MQ=255)
                                                                                  tttATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCg                                            <  2:1286346/137‑1 (MQ=255)
                                                                                  tttATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGt                                           >  1:524767/1‑138 (MQ=255)
                                                                                  tttATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGt                                           >  1:1097756/1‑138 (MQ=255)
                                                                                    tATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTg                                                                                                                            <  1:681766/55‑1 (MQ=255)
                                                                                    tATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTg                                                                                                                            >  2:681766/1‑55 (MQ=255)
                                                                                     aTGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAAt                                        <  2:719736/138‑1 (MQ=255)
                                                                                     aTGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAAt                                        <  2:524767/138‑1 (MQ=255)
                                                                                      tgtgGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATc                                       <  1:368706/138‑1 (MQ=255)
                                                                                          ggATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGt                                   <  1:952256/138‑1 (MQ=255)
                                                                                          ggATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTa                                  >  2:1137988/1‑139 (MQ=255)
                                                                                            aTTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTa                                  >  2:927725/1‑137 (MQ=255)
                                                                                               aCGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAgg                              <  1:1137988/138‑1 (MQ=255)
                                                                                                     gCATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGcc                         >  1:740526/1‑137 (MQ=255)
                                                                                                      cATTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAg                      >  2:979918/1‑139 (MQ=255)
                                                                                                       aTTGTGCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGt                     <  1:969324/139‑1 (MQ=255)
                                                                                                           tgCGCACAACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGAt                 >  2:1029900/1‑139 (MQ=255)
                                                                                                               cacaACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGt                     <  2:262543/131‑1 (MQ=255)
                                                                                                               cacaACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGt                     >  1:262543/1‑131 (MQ=255)
                                                                                                                 caACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCg           <  2:452549/139‑1 (MQ=255)
                                                                                                                 caACGAAGCGCCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCg           <  2:1170739/139‑1 (MQ=255)
                                                                                                                    cGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCtg        >  1:699462/1‑139 (MQ=255)
                                                                                                                          gcCTGGAACGCACCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAat  <  2:940353/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
GTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCCATTTCCGCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAAT  >  NZ_CP009273/2705032‑2705292

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