Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,513,930 A→T 100% L413F (TTA→TTT USA300HOU_1379 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,513,9300AT100.0% 145.9 / NA 47L413F (TTA→TTTUSA300HOU_1379hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (26/21);  total (26/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAC  >  CP000730/1513801‑1514064
                                                                                                                                 |                                                                                                                                      
tgtatgtaTCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAgg                                                                                                                             <  2:657878/141‑1 (MQ=255)
          gCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATg                                                                                                                   >  1:538833/1‑141 (MQ=255)
            tttATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCTCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGca                                                                                                                 >  2:217257/1‑140 (MQ=255)
            tttATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGcc                                                                                                                 <  1:736028/141‑1 (MQ=255)
            tttATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTAGTGTAGTAGGCTTATTCATccc                                                                                                                 >  2:324541/1‑141 (MQ=255)
             ttattaGTAGGCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                                                                <  2:589101/141‑1 (MQ=255)
             ttattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                                                                >  2:93901/1‑141 (MQ=255)
              tattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                                                                               >  1:581080/1‑141 (MQ=255)
               attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatt                                                                                                              >  2:619055/1‑141 (MQ=255)
               attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatt                                                                                                              >  1:55322/1‑141 (MQ=255)
               attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatt                                                                                                              >  2:50820/1‑141 (MQ=255)
               attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatt                                                                                                              >  1:385958/1‑141 (MQ=255)
                ttaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAAACAGCATCTATATCACTGATACTATGTTTTATTGTATTCGGATTATTAATGccgtga                                                                                                             >  2:254743/1‑137 (MQ=255)
                  agtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                                                                           >  1:273185/1‑141 (MQ=255)
                  agtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                                                                           <  1:648035/141‑1 (MQ=255)
                  aggagCAATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                                                                           <  2:376659/138‑1 (MQ=255)
                   gtagCCATTGGCATGGGGTTTGTCTCTTGGGTAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                                                                          >  2:658277/1‑141 (MQ=255)
                   gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                                                                          >  2:677417/1‑141 (MQ=255)
                   gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                                                                          >  2:677507/1‑141 (MQ=255)
                             gCATGGGGGTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTAACACATCGCGTAAAAAAACCAGCATCTATATCACTCATCATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCAc                                                                                                <  1:5513/141‑1 (MQ=255)
                              cATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAATGTACCACATCGCGTAAAAAAACCAGTATCTATATCACTCATAATATGTTTTGTTGTATTACGATTCCTGACCCCATTATTTGGATCCTCACt                                                                                               <  1:571216/141‑1 (MQ=255)
                                tGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                                                                             <  1:714711/141‑1 (MQ=255)
                                 ggtgTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTACTCCCATTATTTGGATTATCACTTAt                                                                                            <  1:771529/138‑1 (MQ=255)
                                 ggggTTTGTTTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTAt                                                                                            <  2:573326/141‑1 (MQ=255)
                                    gTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTCTCCCTTATCCt                                                                                         >  1:91698/1‑141 (MQ=255)
                                      ttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAACCCAGAATCTATATCACTCATGATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGTTTATCACTTTTCCttt                                                                                       >  2:2651/1‑140 (MQ=255)
                                       tGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACAGCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGt                                                                                      <  1:279825/141‑1 (MQ=255)
                                           tctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAttt                                                                                  >  1:84445/1‑141 (MQ=255)
                                                gATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATCTCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                                                                             <  1:522361/141‑1 (MQ=255)
                                                        gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                                                                     >  1:776072/1‑141 (MQ=255)
                                                                 aTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTAAATTAATATtatata                                                            >  1:609346/1‑141 (MQ=255)
                                                                     tGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaa                                                        <  1:539677/141‑1 (MQ=255)
                                                                            aaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                                 <  1:69654/141‑1 (MQ=255)
                                                                             acgTACCACATCGCCTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTCTTAATGCCATTATTTGGATTACCACTTATCCTTGTATTTATAAGTGAATTAATATCATATATTAAAGAtcgtc                                                <  1:267304/139‑1 (MQ=255)
                                                                                     caTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTaaa                                        >  1:671827/1‑141 (MQ=255)
                                                                                      aTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAAc                                       <  1:93901/141‑1 (MQ=255)
                                                                                         gcgTAAAAACACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaat                                    <  1:677886/141‑1 (MQ=255)
                                                                                         gcatAAAAAATCCTGCATCTTTATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaat                                    <  2:252345/138‑1 (MQ=255)
                                                                                            tAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaataat                                 <  1:773179/141‑1 (MQ=255)
                                                                                            tAAAAAAACCAGCATCTATATCACGCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaataat                                 <  1:234899/141‑1 (MQ=255)
                                                                                              aaaaaaaaaaaaTCAATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGc                               <  2:400095/130‑1 (MQ=255)
                                                                                                aaaaCGAGGATCAATAGCTCGCATAATATGTTGTATTGTATTGCGATTATTAATTTCTTTTTTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGGTCGTCGTGCTAAACAATAATGCaa                             >  2:625316/1‑140 (MQ=255)
                                                                                                                 cACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGa            >  2:520320/1‑141 (MQ=255)
                                                                                                                         atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATCTTATATGTTAAAGAGCGTCGTACTAAACAATAATGCACTTAAAGTTTTGAAGTGCCGAAAatt    >  2:21597/1‑141 (MQ=255)
                                                                                                                         atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAttt    >  2:774987/1‑141 (MQ=255)
                                                                                                                          tatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTa   >  2:132406/1‑141 (MQ=255)
                                                                                                                           atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTaa  >  2:383835/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAC  >  CP000730/1513801‑1514064

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