Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 980,640 G→A L176L (CTG→TTG)  aspC ← aspartate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273980,6400GA95.0% 135.9 / ‑0.8 40L176L (CTG→TTG) aspCaspartate transaminase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (16/22);  total (17/23)
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.69e-01

CGGGCAAAACCCTGGTAAGCGAAGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAGCACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACGCTCTTATGGTTCGG  >  NZ_CP009273/980510‑980775
                                                                                                                                  |                                                                                                                                       
cGGGCAAAACCCTGGTAAGCGAAGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAGCACTACGt                                                                                                                                 >  1:656315/1‑139 (MQ=255)
  ggCAAAACCCTGGTAAGCGAAGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCg                                                                                                                               <  2:370890/139‑1 (MQ=255)
         cccTGGTAAGCGAAGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCt                                                                                                                        <  2:729101/139‑1 (MQ=255)
                     aaGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCttcattca                                                                                                            <  2:683344/139‑1 (MQ=255)
                     aaGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCttcattca                                                                                                            <  1:676114/139‑1 (MQ=255)
                       gTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAgg                                                                                                          <  2:366845/139‑1 (MQ=255)
                       gTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAgg                                                                                                          <  2:907640/139‑1 (MQ=255)
                                  gTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCttcattca                                                                                                            >  2:497297/1‑126 (MQ=255)
                                  gTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCttcattca                                                                                                            <  1:497297/126‑1 (MQ=255)
                                  gTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAg                                                                                               <  2:1076220/139‑1 (MQ=255)
                                        aGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATc                                                                                         <  2:662278/139‑1 (MQ=255)
                                                    cGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAgt                                                                             >  2:956762/1‑139 (MQ=255)
                                                              tgCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTcc                                                                   <  1:907361/139‑1 (MQ=255)
                                                                ccAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGc                                                                 <  1:425111/139‑1 (MQ=255)
                                                                 cAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGca                                                                <  2:294982/139‑1 (MQ=255)
                                                                  aGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAGCACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGcat                                                               <  2:656315/139‑1 (MQ=255)
                                                                        tGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTTATCTGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCataa                                                          >  2:492596/1‑138 (MQ=255)
                                                                                          gggTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACtt                                       >  2:540404/1‑139 (MQ=255)
                                                                                             tCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTc                                      >  2:51257/1‑137 (MQ=255)
                                                                                             tCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCa                                    <  2:370621/139‑1 (MQ=255)
                                                                                                aTACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTcc                                                                   <  1:557732/105‑1 (MQ=255)
                                                                                                aTACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTcc                                                                   >  2:557732/1‑105 (MQ=255)
                                                                                                aTACCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGAc                                 >  1:814237/1‑139 (MQ=255)
                                                                                                  aCCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGc                                                                                                                              <  2:349441/44‑1 (MQ=255)
                                                                                                  aCCGGTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGc                                                                                                                              >  1:349441/1‑44 (MQ=255)
                                                                                                     ggTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTg                              >  2:976329/1‑137 (MQ=255)
                                                                                                     ggTTGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCa                            <  2:437471/139‑1 (MQ=255)
                                                                                                       ttGGGTTATGGCAGCAGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCa                            >  1:238571/1‑137 (MQ=255)
                                                                                                                 gcagcaGGCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTAAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACg                <  1:492596/139‑1 (MQ=255)
                                                                                                                 gcagcaGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACg                <  2:502727/139‑1 (MQ=255)
                                                                                                                 gcagcaGCCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACg                <  1:956762/139‑1 (MQ=255)
                                                                                                                    gcagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGa                                                                         >  2:113382/1‑79 (MQ=255)
                                                                                                                    gcagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGa                                                                         <  1:113382/79‑1 (MQ=255)
                                                                                                                    gcagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGa                  <  1:979184/134‑1 (MQ=255)
                                                                                                                    gcagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGa                  >  2:979184/1‑134 (MQ=255)
                                                                                                                    gcagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACGctc             >  2:286088/1‑139 (MQ=255)
                                                                                                                    gcagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACGctc             >  2:1012584/1‑139 (MQ=255)
                                                                                                                     cagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCaa                                                                                 <  2:656076/70‑1 (MQ=255)
                                                                                                                     cagcCATGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCaa                                                                                 >  1:656076/1‑70 (MQ=255)
                                                                                                                          aTGGAACAACACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACGCTCTTATgg       >  2:1186713/1‑139 (MQ=255)
                                                                                                                               acaacaCTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACGCTCTTATGGTTCgg  <  1:976329/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CGGGCAAAACCCTGGTAAGCGAAGTCAAACAGCGGTAACCAGCCTTTCTCAACGGAGAGTTGTGCCAGTGTTTGCCATTGTTCCAGCGTAGGGTCGATACCGGTTGGGTTATGGCAGCAGCCATGGAACAGCACTACGTCGCCAGCCTGAGCTTCATTCAGGCTGTTAATCAGTGCATCGAAGTCAAGAGTGTGATTTTCCGCATCATAATAAGCGTATTCACGAACTTCCAGACCTGCAGAGTTAAAGACGCTCTTATGGTTCGG  >  NZ_CP009273/980510‑980775

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