Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,179,484 A→G E30G (GAA→GGA)  rcnA → nickel/cobalt efflux protein RcnA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,179,4840AG96.8% 103.4 / NA 31E30G (GAA→GGA) rcnAnickel/cobalt efflux protein RcnA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (14/16);  minor base T (0/1);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGAACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGTTAATTGCCTTTGGCGGGA  >  NZ_CP009273/2179354‑2179618
                                                                                                                                  |                                                                                                                                      
gAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAggg                                                                                                                                 >  2:100957/1‑138 (MQ=255)
   tAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACt                                                                                                                             <  2:292119/139‑1 (MQ=255)
          aGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAAc                                                                                                                      <  1:100957/139‑1 (MQ=255)
           gATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAAc                                                                                                                      >  1:572913/1‑138 (MQ=255)
            aTTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACga                                                                                                                    >  1:268145/1‑139 (MQ=255)
             ttGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAggg                                                                                                                                 >  1:654309/1‑125 (MQ=255)
             ttGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAggg                                                                                                                                 <  2:654309/125‑1 (MQ=255)
             ttGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACgat                                                                                                                   <  2:268145/139‑1 (MQ=255)
                              aGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTatca                                                                                                  >  1:715671/1‑139 (MQ=255)
                                          aTGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAg                                                                                      >  2:729965/1‑139 (MQ=255)
                                          aTGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAg                                                                                      >  1:850172/1‑139 (MQ=255)
                                                            cttcttCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATggcg                                                                                                            <  1:58528/99‑1 (MQ=255)
                                                            cttcttCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATggcg                                                                                                            >  2:58528/1‑99 (MQ=255)
                                                               cttcAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGt                                                                   <  2:41969/137‑1 (MQ=255)
                                                                  cAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGc                                                              >  1:605908/1‑139 (MQ=255)
                                                                                   gTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACt                                             <  2:605908/139‑1 (MQ=255)
                                                                                              ccAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCaa                                               >  2:184142/1‑126 (MQ=255)
                                                                                              ccAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCaa                                               <  1:184142/126‑1 (MQ=255)
                                                                                                aGCGCCACCTTACTTGGTGCGCTTCCTGGCCTGGGACCAGGGCACTCAAAAACGATGATGCCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCAGTGATGATCGGACTATCCGCAACTCTTTCGCCTCCCg                                >  2:812560/1‑139 (MQ=255)
                                                                                                  cgcCATCTTACTTGGTGCGCTTCATGGCCTGGTACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCa                              <  2:708380/139‑1 (MQ=255)
                                                                                                   gcCATCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAg                             <  1:106620/139‑1 (MQ=255)
                                                                                                       tCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGt                         >  1:61006/1‑139 (MQ=255)
                                                                                                       tCTTACTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGt                         >  1:472281/1‑139 (MQ=255)
                                                                                                           aCTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGt                         <  2:329192/135‑1 (MQ=255)
                                                                                                            cTTGGTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGt                    <  2:715671/139‑1 (MQ=255)
                                                                                                                gTGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGTTaa                 <  2:259622/138‑1 (MQ=255)
                                                                                                                 tGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGt                         >  2:192015/1‑129 (MQ=255)
                                                                                                                 tGCGCTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGt                         <  1:192015/129‑1 (MQ=255)
                                                                                                                   cgcTTCATGGCCTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGTTAATTg              <  2:650577/138‑1 (MQ=255)
                                                                                                                              cTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGTTAATTGCCTTTGGCGGGa  <  1:69664/139‑1 (MQ=255)
                                                                                                                              cTGGGACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGTTAATTGCCTTTGGCGGGa  <  2:83612/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
GAGTAGAATCAGATTGCCGAATTAATACTAAGAATTATTATCATGACCGAATTTACAACTCTTCTTCAGCAAGGAAACGCCTGGTTCTTCATCCCCAGCGCCATCTTACTTGGTGCGCTTCATGGCCTGGAACCAGGGCACTCAAAAACGATGATGGCGGCGTTTATCATCGCCATCAAAGGCACCATTAAACAAGCGGTGATGCTCGGACTGGCAGCAACTATTTCGCATACCGCAGTGGTCTGGTTAATTGCCTTTGGCGGGA  >  NZ_CP009273/2179354‑2179618

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