Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,600,471 G→A 12.9% A221T (GCC→ACC)  zntA → zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,600,4710GA12.9% 69.3 / 3.5 33A221T (GCC→ACC) zntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base A (2/2);  total (17/16)
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.85e-01

TCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGAAGTCGCCGCAGGTGGGCG  >  CP009273/3600334‑3600613
                                                                                                                                         |                                                                                                                                              
tCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCcgcgcg                                                                                                                                     >  1:175014/1‑149 (MQ=255)
  aaaTCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCcgcgcgcg                                                                                                                                   >  2:106220/1‑149 (MQ=255)
   aaTCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGt                                                                                                                                  >  1:49740/1‑149 (MQ=255)
   aaTCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCGGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCcgcgcgcgg                                                                                                                                  >  1:376665/1‑148 (MQ=255)
      ccGGCAGCTACTTCGCCATTGAAACCGTAATGACCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCcgcgcgctgcag                                                                                                                               <  2:56917/149‑6 (MQ=255)
                 ttCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAgcgc                                                                                                                    <  1:306188/149‑1 (MQ=255)
                  tCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGcgcg                                                                                                                   >  1:373263/1‑149 (MQ=255)
                  tCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATAGGCGCAACTGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGcgcg                                                                                                                   <  1:130779/149‑1 (MQ=255)
                       aTTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAAt                                                                                                              <  1:75085/149‑1 (MQ=255)
                           aaaCCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGgcg                                                                                                          >  1:199769/1‑149 (MQ=255)
                              ccTTAATGAGCGTAGCCGCTATTGGGGCACGGTTTATTGGCGCAACGGCTGAAGCTGCGATGGGGGTGCTGCTGGTTTTGGTTGGTGAACCGCTGGAAGGCGCGGCCGCCAGCCGCGCGCGGCAGGGCGGTAGCGCGTTAATGGCGCTg                                                                                                       >  2:242152/1‑149 (MQ=255)
                                  aaTGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAAc                                                                                                   <  2:270457/149‑1 (MQ=255)
                                     gggCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAg                                                                                                <  2:246860/147‑1 (MQ=255)
                                       gCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGaa                                                                                              <  1:106220/149‑1 (MQ=255)
                                                         cACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGt                                                                            <  2:293173/149‑1 (MQ=255)
                                                           cTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTaa                                                                          <  1:321948/149‑1 (MQ=255)
                                                              tttATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCACCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAggg                                                                       >  2:18016/1‑149 (MQ=255)
                                                              tttATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCACCAGCCGCGCGCGTCAGGGCGTTAGCGCGATAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAggg                                                                       >  2:240440/1‑149 (MQ=255)
                                                                        cAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGaa                                                             >  1:127351/1‑149 (MQ=255)
                                                                                        gATGGTGTTGCTGCTGTTTTTGGGTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGATACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACa                                             <  2:300726/149‑1 (MQ=255)
                                                                                            gtgtTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCt                                         >  1:238773/1‑149 (MQ=255)
                                                                                                         ttttGAGTTGTCAACGCCCGGAAGGTCGGTCCCCCAGCCGCCCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCCGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGAt                            <  2:15978/149‑1 (MQ=255)
                                                                                                           ttGATTGGTGAACGACTGGAAGGCTGGGCCACCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGAtgt                          <  1:18016/149‑1 (MQ=255)
                                                                                                           ttGATTGGTGAACGACTGGAAGGCTGGGCCACCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGCGGCGCTTAACAGCCTGCGCCCTGGCGAtgt                          <  1:240440/149‑1 (MQ=255)
                                                                                                                tGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTg                     <  1:21693/149‑1 (MQ=255)
                                                                                                                  gTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGaa                   >  1:288010/1‑149 (MQ=255)
                                                                                                                  gTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGaa                   <  2:397114/149‑1 (MQ=255)
                                                                                                                       cGACTGGAAGGCAGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCCCGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGAAGTcgc              <  1:80223/149‑1 (MQ=255)
                                                                                                                         aCTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGCTGTGATTGAAGTcgccg            >  2:102934/1‑149 (MQ=255)
                                                                                                                           tGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGAAGTCGCCGCa          >  1:125246/1‑149 (MQ=255)
                                                                                                                                 gCTGGGCCGCCAGCCGCCCGCGTCAGGGCGTTAGCGCGTTGGCGCCGCTGAAACCCGAAAGCTACACAAACTTGCGTAAGGGTGACCGGGAAGAGGGGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGAAGTCGCCGCAGGTggg    >  2:450514/1‑149 (MQ=255)
                                                                                                                                  cTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTCATTGAAGTCGCCGCAGGTGGGc   >  1:72224/1‑149 (MQ=255)
                                                                                                                                   tGGGCCGCGAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGGGAGCGGGAAGAGGTGGCACTTAACGGCCTGCGACCTGTCGATGTGATTGAACTCCCCGCACGTGGGCg  >  2:10993/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                              
TCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTATTGGTGCACTGTTTATTGGCGCAACGGCTGAAGCTGCGATGGTGTTGCTGCTGTTTTTGATTGGTGAACGACTGGAAGGCTGGGCCGCCAGCCGCGCGCGTCAGGGCGTTAGCGCGTTAATGGCGCTGAAACCAGAAACCGCCACGCGCCTGCGTAAGGGTGAGCGGGAAGAGGTGGCGATTAACAGCCTGCGCCCTGGCGATGTGATTGAAGTCGCCGCAGGTGGGCG  >  CP009273/3600334‑3600613

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