Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,613,222 C→G R238P (CGC→CCC)  lplA ← lipoate‑‑protein ligase LplA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,613,2220CG95.5% 137.7 / ‑0.3 44R238P (CGC→CCC) lplAlipoate‑‑protein ligase LplA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (23/19);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

ACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGC  >  NZ_CP009273/4613141‑4613307
                                                                                 |                                                                                     
aCGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAg                                                                               >  1:936218/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTc                                                                          >  1:1235609/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTc                                                                          >  2:139937/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTc                                                                          >  2:1127300/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTc                                                                          >  1:245846/1‑90 (MQ=255)
               aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt                                                                >  1:1362770/1‑90 (MQ=255)
               aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt                                                                >  2:1351272/1‑90 (MQ=255)
                 gCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTgg                                                              >  1:1266758/1‑90 (MQ=255)
                   gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                            <  2:936218/90‑1 (MQ=255)
                   gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                            <  2:911077/90‑1 (MQ=255)
                   gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                            >  2:761390/1‑90 (MQ=255)
                    ttCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCaa                                                           >  2:1009173/1‑90 (MQ=255)
                       aTCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTc                                                        <  2:426984/90‑1 (MQ=255)
                          cagcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTgg                                                     <  2:127850/90‑1 (MQ=255)
                           agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGc                                                    >  2:938642/1‑90 (MQ=255)
                           agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGc                                                    >  2:825799/1‑90 (MQ=255)
                            gcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCg                                                   >  1:841906/1‑90 (MQ=255)
                               gATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGttt                                                <  2:1362770/90‑1 (MQ=255)
                                   cGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGtt                                            <  1:139937/90‑1 (MQ=255)
                                   cGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGtt                                            >  1:191255/1‑90 (MQ=255)
                                      gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg                                         >  2:686704/1‑90 (MQ=255)
                                      gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg                                         >  1:279789/1‑90 (MQ=255)
                                         tGCCGGAGCCTGACCGAAGTTCCATTCTCAGCTACTCTGGGGGGCAAAGGGTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGga                                      >  1:1256032/1‑90 (MQ=255)
                                         tGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGga                                      >  2:1207021/1‑90 (MQ=255)
                                          gCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGgag                                     >  2:1567594/1‑90 (MQ=255)
                                              gAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAgatga                                 <  1:1046612/90‑1 (MQ=255)
                                                     aCCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGc                          <  2:852449/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAGGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:807206/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:955444/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  1:1422750/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:465275/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:180084/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:1109161/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:1067578/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  1:825799/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:245846/90‑1 (MQ=255)
                                                                aTTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCt               <  2:1266758/90‑1 (MQ=255)
                                                                   cccAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGc            >  2:1574248/1‑90 (MQ=255)
                                                                     cAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGtttt                                               >  2:896217/1‑53 (MQ=255)
                                                                     cAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGtttt                                               <  1:896217/53‑1 (MQ=255)
                                                                       gCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATa        >  2:601585/1‑90 (MQ=255)
                                                                             ctGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGc  >  1:1547167/1‑90 (MQ=255)
                                                                             ctGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGc  >  1:742495/1‑90 (MQ=255)
                                                                             ctGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGAAGATTTCCGCTTCCACGCGCTCGCCATAATGGGc  <  2:1534533/90‑1 (MQ=255)
                                                                                 |                                                                                     
ACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGC  >  NZ_CP009273/4613141‑4613307

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