Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,138,500 T→A 16.9% F300Y (TTC→TAC) ‡ czcC → CzcC family metal RND transporter outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,138,5000TA16.9% 52.9 / 3.9 24F300Y (TTC→TAC) ‡czcCCzcC family metal RND transporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (11/8);  new base A (2/2);  total (14/10)
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.11e-01

CCATGACCCGTGAAGCACAAGCCCAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAGCTGACCTGGCGGAATACCAACGCCTACAACGGGCAG  >  NC_002947/6138356‑6138643
                                                                                                                                                |                                                                                                                                               
ccATGACCCGTGAAGCACAAGCCCAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAAGTCAGCAACCTGGTAATTGGCGAGGTCAGTTTTCAATTGCCGCTGTTCAGCg                                                                                                                                            >  2:40714/1‑150 (MQ=255)
 cATGACCCGTGAAGCACAAGCCCAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCgg                                                                                                                                           >  2:329884/1‑150 (MQ=255)
      cccGTGAAGCACAAGCCCAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGAACAGCGGcaatt                                                                                                                                      <  2:488335/150‑5 (MQ=255)
                  aaGCCCAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCga                                                                                                                          <  1:40714/150‑1 (MQ=255)
                      ccATGTGCATCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGAACAGCGGcaa                                                                                                                                        <  2:489842/132‑3 (MQ=255)
                      ccAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGAACAGCGGcaa                                                                                                                                        >  1:489842/1‑130 (MQ=255)
                         ggTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCgc                                                                                                                   <  2:155371/150‑1 (MQ=255)
                                 agGCCATCGCCGAGAAAAAGTCGGATTGGGGCACGACGGTCTACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGcc                                                                                                           <  2:501416/150‑1 (MQ=255)
                                      aTCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCgat                                                                                                                         >  1:319703/1‑131 (MQ=255)
                                      aTCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCgat                                                                                                                         <  2:319703/131‑1 (MQ=255)
                                           cGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGt                                                                                                 >  1:244355/1‑150 (MQ=255)
                                                aaaaGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTc                                                                                            >  2:238352/1‑150 (MQ=255)
                                                        gATTGGGGCTGGCAGGTCGACTAGCAACGCCGGGGGCCCGCATTAAGCAACATGGTCAGTGTCGAGGTCAGATTTCAATTGCCGCTGTTCAGCGGCCCGAGGTAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGTGTCGGCTggag                                                                                    >  2:79544/1‑150 (MQ=255)
                                                                        tCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCgg                                                                    <  1:75691/150‑1 (MQ=255)
                                                                                  aCGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGa                                                          <  1:238352/150‑1 (MQ=255)
                                                                                                     aGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGaa                                       >  2:425194/1‑150 (MQ=255)
                                                                                                     aGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGaa                                       >  2:461102/1‑150 (MQ=255)
                                                                                                      gCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAg                                      <  1:329884/150‑1 (MQ=255)
                                                                                                                     tCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAGCTGACCTGGCGGAAt                       >  2:247202/1‑150 (MQ=255)
                                                                                                                     tCCAGGGCAGTTTGCAATTGCCGCTGTGCAGCGGCTCGCGGCAAGACCCGACCATCGCTGCTCCCCCGGCCCACGTGCGACCGCTGGACGACCAACACGAAGCGGCCCTGCGTGAGCAGACCGCGCAAAGCGAAGCTGACCTGGCGTaaa                       >  2:485725/1‑149 (MQ=255)
                                                                                                                          gTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGAGGCACTGCGTGAGCACACCGCGCAACTCGTAGCTGACCTGGCGGAATACCaa                  >  2:55744/1‑150 (MQ=255)
                                                                                                                             aGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAGCTGACCTGGCGGAATACCAACGc               <  1:385594/150‑1 (MQ=255)
                                                                                                                                     cgagCCGCTGAACAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAGCTGACCTGGCGGAATACCAACGCCTACAACg       >  1:304090/4‑150 (MQ=255)
                                                                                                                                          cGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTGCGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAGCTGACCTGGCGGAATACCAACGCCTACAACCGGCAg  >  1:236726/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CCATGACCCGTGAAGCACAAGCCCAGGTGCAGCAGGCCATCGCCGAGAAAAAGTCGGATTGGGGCTGGCAGGTCGACTACCAACGCCGCGGGCCCGAATTCAGCAACATGGTCAGTGTCCAGGTCAGTTTTCAATTGCCGCTGTTCAGCGGCTCGCGGCAAGACCCGATGATCGCTGCTCGCCGGGCCCAGGTGCGTCAGCTGGAGGACGAACAGGAAGCGGCACTGCGTGAGCACACCGCGCAACTCGAAGCTGACCTGGCGGAATACCAACGCCTACAACGGGCAG  >  NC_002947/6138356‑6138643

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