Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,303,545 T→C D328G (GAT→GGT)  pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,303,5450TC83.3% 84.4 / 13.3 36D328G (GAT→GGT) pnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (18/12);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGTTCGTTTTCGTCCAGGGTTTCGT  >  NZ_CP009273/3303414‑3303674
                                                                                                                                   |                                                                                                                                 
gTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGAc                                                                                                                             >  1:183094/1‑138 (MQ=255)
 ttCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCa                                                                                                                           <  2:449435/139‑1 (MQ=255)
            gTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATc                                                                                                                <  2:46583/139‑1 (MQ=255)
             ttCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCt                                                                                                               >  1:780969/1‑139 (MQ=255)
                   aaGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCa                                                                                                         <  2:627523/139‑1 (MQ=255)
                       aCCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGAc                                                                                                     <  2:183094/139‑1 (MQ=255)
                           gcgcGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTc                                                                                                 >  1:848972/1‑139 (MQ=255)
                            cgcgTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTc                                                                                                 >  1:813624/1‑138 (MQ=255)
                                        aGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCGCGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCg                                                                                    <  1:40672/139‑1 (MQ=255)
                                               aGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGc                                                                               >  2:312465/1‑137 (MQ=255)
                                                gCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCa                                                                                                                           <  2:85547/92‑1 (MQ=255)
                                                gCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCa                                                                                                                           >  1:85547/1‑92 (MQ=255)
                                                  gTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTa                                                                          <  1:540003/139‑1 (MQ=255)
                                                    tGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCAGGGATCATATCTTTTTAACGACCGTCGTTACGCGGTCCGCCTGccctgccc                                                                        >  1:326846/1‑133 (MQ=255)
                                                     gCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTAcg                                                                        <  2:93485/138‑1 (MQ=255)
                                                      cGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCa                                                                                                                           >  2:272551/1‑86 (MQ=255)
                                                      cGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCa                                                                                                                           <  1:272551/86‑1 (MQ=255)
                                                      cGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTAcgcg                                                                      >  1:842778/1‑139 (MQ=255)
                                                      cGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTAcgcg                                                                      >  1:208485/1‑139 (MQ=255)
                                                      cGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTAcgcg                                                                      >  2:855400/1‑139 (MQ=255)
                                                                   cTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACg                                                         <  1:720215/139‑1 (MQ=255)
                                                                      cGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGttt                                                      >  1:451418/1‑139 (MQ=255)
                                                                      cGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGtt                                                       <  2:498307/138‑1 (MQ=255)
                                                                                   gTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGt                                         >  1:842043/1‑139 (MQ=255)
                                                                                           cgcAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAAtt                                 >  2:657518/1‑139 (MQ=255)
                                                                                           cgcAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAAtt                                 >  1:140285/1‑139 (MQ=255)
                                                                                           cgcAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAAtt                                 >  1:483163/1‑139 (MQ=255)
                                                                                           cgcAGAACCGTGAGTACGAGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAAtt                                 >  1:702100/1‑139 (MQ=255)
                                                                                                     tGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGtt                       >  1:410125/1‑139 (MQ=255)
                                                                                                     tGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGtt                       >  2:503292/1‑139 (MQ=255)
                                                                                                     tGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGtt                       >  2:549834/1‑139 (MQ=255)
                                                                                                     tGAGTACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGtt                       >  2:379806/1‑139 (MQ=255)
                                                                                                         tACGCGGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGTTCGtt                   <  1:312465/139‑1 (MQ=255)
                                                                                                           cgcgGCAGCACGCCAGTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGTTCGttt                  <  2:601623/138‑1 (MQ=255)
                                                                                                                cagcaCGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGTTCGTTTTCGTcc            <  1:657518/139‑1 (MQ=255)
                                                                                                                          gTACGCACACCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGTTCGTTTTCGTCCAGGGTTTCGt  <  1:751007/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
GTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTACGAACAACGTTTTTCTCGATCGCGTGCAGAATTTCACCCAGTTCGTTTTCGTCCAGGGTTTCGT  >  NZ_CP009273/3303414‑3303674

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