Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,671,907 G→A H46Y (CAC→TAC)  USA300HOU_2527 ← possible thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,671,9070GA97.6% 129.9 / ‑5.0 41H46Y (CAC→TAC) USA300HOU_2527possible thioesterase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/19);  total (22/19)
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.99e-01

TCTTGCTTGCTTTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTGCTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCTG  >  CP000730/2671775‑2672035
                                                                                                                                    |                                                                                                                                
tCTTGCTTGCTTTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATtaa                                                                                                                          >  1:558774/1‑141 (MQ=255)
             aaCTAAATAAATACGACTTACTGGTTTGTCATGTATTGTGACAGGTACTTTATTCTCTGCTTCCGTTCTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTtct                                                                                                             <  1:280797/141‑1 (MQ=255)
                 aaGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTTAGTAATTCATGTACTGATTAATAAACGCTGTTCTCaaa                                                                                                         <  1:873953/141‑1 (MQ=255)
                 aaGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCaaa                                                                                                         <  1:652785/141‑1 (MQ=255)
                 aaGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCaaa                                                                                                         <  1:710152/141‑1 (MQ=255)
                  aGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAAt                                                                                                        <  1:417461/141‑1 (MQ=255)
                    tAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCt                                                                                                      <  1:820394/141‑1 (MQ=255)
                    tAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCt                                                                                                      <  1:13787/141‑1 (MQ=255)
                     aaaTACGACCTACTGGTTCGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCtt                                                                                                     <  1:319183/141‑1 (MQ=255)
                      aaTACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTc                                                                                                    >  1:370415/1‑141 (MQ=255)
                            aCCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCTAACCa                                                                                              >  1:308500/1‑141 (MQ=255)
                               tACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCGAATCTTCGAACCATCt                                                                                           <  1:475434/141‑1 (MQ=255)
                                 cTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTaa                                                                                         >  1:231000/1‑141 (MQ=255)
                                  tGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAat                                                                                        >  1:925305/1‑141 (MQ=255)
                                                    gACAGGTACTTTATACTCTACTTCAACTGGCAGAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAAc                                                                      >  1:682447/1‑141 (MQ=255)
                                                            cTTTATACTCTGCTTCCGTATTCATAAGTATAGGCGAAATGCGTTGATAGATCATTGTTGAGTAATTGATGTGCAGATTAATAGACGCTGGTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTGACGACACCCa                                                              >  1:693282/1‑141 (MQ=255)
                                                                ataCTCTGCCTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGc                                                          <  1:504475/141‑1 (MQ=255)
                                                                      tGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAAt                                                    <  1:358541/141‑1 (MQ=255)
                                                                       gCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATg                                                   >  1:362113/1‑141 (MQ=255)
                                                                            cGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTa                                              >  1:517769/1‑141 (MQ=255)
                                                                               tttCATAAGTATATGTGAAATGCCTTGTTTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAAcc                                           <  1:394868/141‑1 (MQ=255)
                                                                                 tCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGGCGTAACCAt                                         >  1:92748/1‑141 (MQ=255)
                                                                                 tCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAAGAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAAccgt                                         >  1:157449/1‑139 (MQ=255)
                                                                                             gTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAttt                             <  1:34170/141‑1 (MQ=255)
                                                                                                  aTGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAtttctttc                        <  1:904787/141‑1 (MQ=255)
                                                                                                  aTGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAtttctttc                        >  1:758312/1‑141 (MQ=255)
                                                                                                    gCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATGTTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTtctc                      >  1:26004/1‑141 (MQ=255)
                                                                                                    gCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACAAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTtctc                      >  1:719214/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCGTTGTTGAGCAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctct                     >  1:479666/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGCTTTCTTTctct                     >  1:873496/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctct                     >  1:847970/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctct                     >  1:364440/1‑141 (MQ=255)
                                                                                                      cTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTg                    >  1:187064/1‑141 (MQ=255)
                                                                                                      cTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTg                    >  1:139711/1‑141 (MQ=255)
                                                                                                       ttgattgaTCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGc                   >  1:113968/1‑141 (MQ=255)
                                                                                                       ttgattgaTCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGc                   >  1:189861/1‑141 (MQ=255)
                                                                                                                   tgttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGAcc       <  1:702265/141‑1 (MQ=255)
                                                                                                                    gttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCt      <  1:331832/141‑1 (MQ=255)
                                                                                                                     ttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCtt     <  1:187109/141‑1 (MQ=255)
                                                                                                                     ttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCtt     <  1:142304/141‑1 (MQ=255)
                                                                                                                        aGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCtg  <  1:252844/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                
TCTTGCTTGCTTTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTGCTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCTG  >  CP000730/2671775‑2672035

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