Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,367,563 T→C intergenic (‑24/+65) blc ← / ← ampC lipocalin Blc/BlaEC family class C beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,367,5630TC84.0% 51.6 / 4.8 25intergenic (‑24/+65)blc/ampClipocalin Blc/BlaEC family class C beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (9/12);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA  >  NZ_CP009273/4367477‑4367647
                                                                                      |                                                                                    
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTa                                                                                   <  2:280909/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:505515/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:491389/90‑1 (MQ=255)
           cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                         >  2:562465/1‑89 (MQ=255)
           cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                         <  1:562465/89‑1 (MQ=255)
              cGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAg                                                                     <  2:179832/90‑1 (MQ=255)
              cGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAg                                                                     <  2:231461/90‑1 (MQ=255)
                 ccAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCaa                                                                  <  1:650495/90‑1 (MQ=255)
                          cAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCg                                                         >  2:562817/1‑90 (MQ=255)
                              tGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGgaga                                                     >  2:18261/1‑90 (MQ=255)
                                                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                 <  1:38538/90‑1 (MQ=255)
                                                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                 >  2:104153/1‑90 (MQ=255)
                                                      cAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa                             >  1:397582/1‑90 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                            <  1:780602/90‑1 (MQ=255)
                                                            cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGCAAATTCGGACCCGATGGAAt                       <  1:526418/90‑1 (MQ=255)
                                                            cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                       <  1:104153/90‑1 (MQ=255)
                                                             aTAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt                      >  2:248303/1‑90 (MQ=255)
                                                             aTAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt                      >  2:580817/1‑90 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                     <  1:102909/90‑1 (MQ=255)
                                                                     ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa              <  1:248303/90‑1 (MQ=255)
                                                                      tCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag             >  2:625568/1‑90 (MQ=255)
                                                                        cTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag           <  1:832123/90‑1 (MQ=255)
                                                                        cTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag           >  2:470564/1‑90 (MQ=255)
                                                                        cTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag           <  2:761959/90‑1 (MQ=255)
                                                                          tACTGGTTTTTTCCTAAGTGTAGATGAAAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCg         >  2:702809/1‑90 (MQ=255)
                                                                                 ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGaa  <  1:331407/90‑1 (MQ=255)
                                                                                 ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGaa  <  2:710857/90‑1 (MQ=255)
                                                                                      |                                                                                    
AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA  >  NZ_CP009273/4367477‑4367647

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