Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,093,562 A→G I396I (ATT→ATC gndA ← NADP‑dependent phosphogluconate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,093,5620AG90.0% 109.0 / 6.2 40I396I (ATT→ATCgndANADP‑dependent phosphogluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (25/11);  total (27/13)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

AGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTT  >  NZ_CP009273/2093431‑2093696
                                                                                                                                   |                                                                                                                                      
aGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCt                                                                                                                                 >  1:286465/1‑139 (MQ=255)
aGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCt                                                                                                                                 >  1:480406/1‑139 (MQ=255)
                        cGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGt                                                                                                         >  2:199234/1‑139 (MQ=255)
                           tAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAg                                                                                                      >  1:744694/1‑139 (MQ=255)
                           tAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAg                                                                                                      >  2:440857/1‑139 (MQ=255)
                             gCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCg                                                                                                    <  1:440857/139‑1 (MQ=255)
                                gTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTg                                                                                                                               <  2:610494/109‑1 (MQ=255)
                                gTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTg                                                                                                                               >  1:610494/1‑109 (MQ=255)
                                     aaTAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCtgtg                                                                                             >  1:457423/1‑138 (MQ=255)
                                          gCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGAtttt                                                                                       <  1:46955/139‑1 (MQ=255)
                                           cAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTc                                                                                      >  2:106144/1‑139 (MQ=255)
                                           cAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTc                                                                                      >  1:533221/1‑139 (MQ=255)
                                            aaCCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCg                                                                                     >  1:812203/1‑139 (MQ=255)
                                            aaCCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCg                                                                                     >  2:330279/1‑139 (MQ=255)
                                              ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAg                                                                                                                 >  2:277173/1‑109 (MQ=255)
                                              ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAg                                                                                                                 <  1:277173/109‑1 (MQ=255)
                                              ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCgg                                                                                    >  1:488143/1‑138 (MQ=255)
                                                    cGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAg                                                                              >  1:10631/1‑138 (MQ=255)
                                                            gTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCaa                                                                                                              <  2:111618/98‑1 (MQ=255)
                                                            gTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCaa                                                                                                              >  1:111618/1‑98 (MQ=255)
                                                               ggAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGAttt                                                                  <  1:625839/139‑1 (MQ=255)
                                                                            cGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTg                                                     <  2:480406/139‑1 (MQ=255)
                                                                                cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATa                                                                                >  2:532446/1‑108 (MQ=255)
                                                                                cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATa                                                                                <  1:532446/108‑1 (MQ=255)
                                                                                cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTgcgc                                                  <  2:533221/138‑1 (MQ=255)
                                                                                 tGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTCCTGCAGGAACTGCGCAc                                                >  1:579399/1‑139 (MQ=255)
                                                                                   tACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACgg                                              <  2:286465/139‑1 (MQ=255)
                                                                                        cATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                         >  2:716370/1‑139 (MQ=255)
                                                                                        cATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                         >  2:611386/1‑139 (MQ=255)
                                                                                           aaGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGatgatg                                        >  2:560556/1‑137 (MQ=255)
                                                                                             gCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGaa                                                        >  1:672854/1‑119 (MQ=255)
                                                                                             gCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGaa                                                        <  2:672854/119‑1 (MQ=255)
                                                                                             gCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagc                                    >  1:704322/1‑139 (MQ=255)
                                                                                             gCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagc                                    >  2:862472/1‑139 (MQ=255)
                                                                                                  gACATCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCa                               >  1:600961/1‑139 (MQ=255)
                                                                                                      tCACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGa                           <  2:600961/139‑1 (MQ=255)
                                                                                                       cACGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGa                           >  2:62843/1‑138 (MQ=255)
                                                                                                         cGCAGCGCCTGCTGGTAGTCATCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGaaaa                        >  2:567283/1‑139 (MQ=255)
                                                                                                                             aTCGGCGATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAg    >  1:755863/1‑139 (MQ=255)
                                                                                                                                ggCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGtt  <  1:611386/138‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
AGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTT  >  NZ_CP009273/2093431‑2093696

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