Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,709,600 (T)5→4 100% coding (1327/1812 nt) USA300HOU_2562 ← possible acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,709,5960T.100.0% 149.6 / NA 33coding (1331/1812 nt)USA300HOU_2562possible acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/13);  total (20/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATTGCGATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCATGTTGTTTGCCAAGTGCA  >  CP000730/2709458‑2709727
                                                                                                                                           |                                                                                                                                   
tatTGCGATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACttt                                                                                                                                   <  2:195041/141‑1 (MQ=255)
  ttGCGATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATAACtttta                                                                                                                                 >  1:228835/1‑140 (MQ=255)
   tGCGATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGTTGACttttca                                                                                                                                >  2:29693/1‑139 (MQ=255)
   tGCGATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGGCACCAATTAATAGTGTTGATGACttttaa                                                                                                                                >  2:266174/1‑139 (MQ=255)
      gATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATTACCGAGTCACCAATTAATAGTGGTGATGACTTTTAatat                                                                                                                            >  1:7363/1‑141 (MQ=255)
         ttCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGGGGTGATGACTTTTAATATCct                                                                                                                         >  1:242452/1‑141 (MQ=255)
          tCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCctc                                                                                                                        >  2:166599/1‑141 (MQ=255)
                     gtAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCg                                                                                                             >  1:29383/1‑141 (MQ=255)
                          aTCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTc                                                                                                        <  1:236353/141‑1 (MQ=255)
                           tCAACCAGTTGCCGTCCAACTTAACCCCAATATTGTGCATATGGTATTTTCTTAGTAAAGACATTAGCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCt                                                                                                       >  1:210821/1‑141 (MQ=255)
                                      ccGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTGTGCTATTCGCTGTCTGCTTATCctct                                                                                            <  2:171450/141‑1 (MQ=255)
                                       cGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCtt                                                                                           >  1:152378/1‑141 (MQ=255)
                                          ccAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCt                                                                                        <  1:131271/141‑1 (MQ=255)
                                             aCTTTACCATCATTTTGTTCATTTGGTATTTTCATAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTAt                                                                                     <  2:19414/141‑1 (MQ=255)
                                               tttACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCt                                                                                   >  1:272908/1‑141 (MQ=255)
                                                                gCATTTGGTATTTTCTTAGTAAAGACATTACCAATATTCACAATGACCGAGTGACCAAATAATAGTGGGGATTACTTTTAATATCCTCTTTGCTATTCGCTGCCTGCTTTACGTCTTTCTTATCTTTTTTAACGCCttttt                                                                   >  1:60542/1‑141 (MQ=255)
                                                                  aTTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTttc                                                                 >  2:198141/1‑141 (MQ=255)
                                                                    ttGGGACTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTttctt                                                               <  1:26142/141‑1 (MQ=255)
                                                                            ttCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTTCTTATCTTTTTTAACGACTTTTttcttgctttttt                                                       >  2:63883/1‑133 (MQ=255)
                                                                              cTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTa                                                     <  2:43380/141‑1 (MQ=255)
                                                                               ttAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGGTGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATACTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTTTTTTaa                                                    >  2:83380/1‑141 (MQ=255)
                                                                                  gTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTTTTTTAAAtt                                                 >  2:125341/1‑141 (MQ=255)
                                                                                       gACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGtc                                            >  1:130154/1‑141 (MQ=255)
                                                                                        aCATTACCAATATCCACCATGACCGAGTCACAAAATAAAAGTGGTGATGACTTTTAATATCCTATTTGCTATTAGCTGTCTGCGTATCCTCTTGCTTATCTTTTTTAACGTCTTTTTTCTTCGTTGTTTTAAATTCCGtct                                           >  2:210034/1‑141 (MQ=255)
                                                                                             aCCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGtcttcttc                                      <  1:287510/141‑1 (MQ=255)
                                                                                              ccAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTGTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGtctttttct                                     >  2:182267/1‑141 (MQ=255)
                                                                                                 atatCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCtctt                                  <  1:287148/141‑1 (MQ=255)
                                                                                                      caccaTGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTttcttc                             <  1:263680/141‑1 (MQ=255)
                                                                                                       cccaTGACCGAGCCCCCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTtcttct                            <  1:178330/140‑1 (MQ=255)
                                                                                                       accaTGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTATTTTTCTCTTtcttct                            <  2:138662/141‑1 (MQ=255)
                                                                                                           tGACCGAGTCACCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGc                        <  1:256551/141‑1 (MQ=255)
                                                                                                                   tCACCAATTAATAGTTGTGATGACTTTTAATATCCTCTTTTCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCAtgtt                <  1:166599/141‑1 (MQ=255)
                                                                                                                     aCCAATTAATAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCGtgttgt              >  2:171880/1‑141 (MQ=255)
                                                                                                                            aaTAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAGATTCCGTCTTCTTCTCTTTCTTCTTGGCTTCATGTTGTTTGCCag       >  1:155583/1‑140 (MQ=255)
                                                                                                                              tAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCATGTTGTTTGCCAAGt     >  1:271062/1‑141 (MQ=255)
                                                                                                                              tAGTGGTGATGACTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCATGTTGTTTGCCAAGt     >  1:64778/1‑141 (MQ=255)
                                                                                                                                 tggtgATGACTTTTGATATACTCTTAGCTATTCGCTGTCTGCTTGTCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCATGTTGTTTGCCAAGttct  >  1:156647/1‑138 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TATTGCGATTTCACAATTGGTGTAGCATCAACGAGTTGCCGTCCAACTTTACCATCAATTTGTGCATTTGGTATTTTCTTAGTAAAGACATTACCAATATCCACCATGACCGAGTCACCAATTAATAGTGGTGATGACTTTTTAATATCCTCTTTGCTATTCGCTGTCTGCTTATCCTCTTGCTTATCTTTTTTAACGACTTTTTTCTTCGTTGTTTTAAATTCCGTCTTCTTCTCTTTCTTCTCGGCTTCATGTTGTTTGCCAAGTGCA  >  CP000730/2709458‑2709727

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