Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,077,206 G→A R659R (CGC→CGT)  fdnG ← formate dehydrogenase‑N subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,077,2060GA95.3% 150.2 / ‑1.5 43R659R (CGC→CGT) fdnGformate dehydrogenase‑N subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (23/18);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

GCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTT  >  NZ_CP009273/4077075‑4077321
                                                                                                                                   |                                                                                                                   
gCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAgg                                                                                                               >  2:588180/1‑138 (MQ=255)
gCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGa                                                                                                              <  1:589391/139‑1 (MQ=255)
 ccAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAgg                                                                                                               >  2:393670/1‑137 (MQ=255)
               ttGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGa                                                                                               >  1:35348/1‑139 (MQ=255)
                  ctctCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGAtct                                                                                            >  2:20144/1‑139 (MQ=255)
                               cttcttCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGtcagtcag                                                                               >  2:11380/1‑139 (MQ=255)
                                     ctGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCa                                                                           <  2:707221/137‑1 (MQ=255)
                                       gAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATcc                                                                       >  2:615670/1‑139 (MQ=255)
                                             ggCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCgggg                                                                 >  1:292515/1‑139 (MQ=255)
                                               cTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCg                                                               >  1:637324/1‑139 (MQ=255)
                                                   tGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGa                                                                                                              <  2:690174/88‑1 (MQ=255)
                                                   tGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGa                                                                                                              >  1:690174/1‑88 (MQ=255)
                                                      ggAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCAc                                                        <  2:35348/139‑1 (MQ=255)
                                                          tGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCttt                                                    >  1:665740/1‑139 (MQ=255)
                                                             cGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTc                                                   >  2:323401/1‑137 (MQ=255)
                                                             cGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGGTCGCGCCACCATGTTCGGCATACATCTTGCGCAAACGCAGGAAGACACCGGAGAGGATCTAGCCGTCCGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTc                                                   >  2:658086/1‑137 (MQ=255)
                                                               tAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCg                                                                                          <  1:714656/96‑1 (MQ=255)
                                                               tAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCg                                                                                          >  2:714656/1‑96 (MQ=255)
                                                                 gTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTg                                              <  2:637324/138‑1 (MQ=255)
                                                                 gTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTg                                              <  2:558987/138‑1 (MQ=255)
                                                                       tgtcatgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTc                                                   <  1:60903/127‑1 (MQ=255)
                                                                       tgtcatgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTc                                                   >  2:60903/1‑127 (MQ=255)
                                                                            tgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAacg                                                                                                                   <  1:4807/58‑1 (MQ=255)
                                                                            tgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAacg                                                                                                                   >  2:4807/1‑58 (MQ=255)
                                                                            tgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCggg                                                                  >  1:346753/1‑107 (MQ=255)
                                                                            tgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCggg                                                                  <  2:346753/107‑1 (MQ=255)
                                                                            tgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAAc                                    >  1:531850/1‑137 (MQ=255)
                                                                            tgtTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAAc                                    >  2:149193/1‑137 (MQ=255)
                                                                              ttCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAg                                 <  1:20144/138‑1 (MQ=255)
                                                                                   cACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGga                                                                                                    >  2:151681/1‑66 (MQ=255)
                                                                                   cACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGga                                                                                                    <  1:151681/66‑1 (MQ=255)
                                                                                   cACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCt                           >  1:633541/1‑139 (MQ=255)
                                                                                        ggTCCGGTTTCGCGCCCCCCTGCTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGCGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGtt                      <  2:327820/139‑1 (MQ=39)
                                                                                        ggTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGt                       >  2:292133/1‑138 (MQ=255)
                                                                                        ggTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGt                       >  1:474545/1‑138 (MQ=255)
                                                                                         gTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGtt                      <  2:39033/138‑1 (MQ=255)
                                                                                         gTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGtt                      <  2:531850/138‑1 (MQ=255)
                                                                                             gggTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGa                 <  1:11380/139‑1 (MQ=255)
                                                                                                   gcgcCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGaa            <  1:323401/138‑1 (MQ=255)
                                                                                                    cgcCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGACACTGTAACCAGCGACCTGAGTTAACGATAGAAcc          >  1:327820/1‑139 (MQ=255)
                                                                                                     gcCACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCg         <  1:766675/139‑1 (MQ=255)
                                                                                                       cACCCTGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGtt       <  2:633541/139‑1 (MQ=255)
                                                                                                            tGTTCGGCATACATCTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCtt  >  2:685525/1‑139 (MQ=255)
                                                                                                                          cTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGtt                      <  2:253364/105‑1 (MQ=255)
                                                                                                                          cTTGCGCAAACGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGtt                      >  1:253364/1‑105 (MQ=255)
                                                                                                                                   |                                                                                                                   
GCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTT  >  NZ_CP009273/4077075‑4077321

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