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,9300AT98.2% 172.7 / ‑3.9 56coding (1239/1341 nt)USA300HOU_1379hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (21/34);  minor base . (1/0);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAAC  >  CP000730/1513794‑1514051
                                                                                                                                         |                                                                                                                         
atttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTg                                                                                                                       <  2:359901/141‑1 (MQ=255)
    atttGTATTGATCGTTGTATTAGCAGCCATTTGTATGGGTTTTGTCGCTTGGATAACGCGTACAAAAAATACTGCAATAAAAGTTCCACATTGCGTATAAATTCCCGCATCTATATCCCTCATAATATGTTTTTTTGTAtt                                                                                                                   <  2:548373/141‑1 (MQ=255)
     tttggttgtatCGCTTTATTAGGAGGCATTGGCATGGGGTTTGTCCCTTGGATGAAGCGTACAAATATTACTGCAATTATAGTAAAAAAAAAAAAAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTa                                                                                                                  <  2:475530/135‑1 (MQ=255)
        ttttgtatCGCGTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGAGAAAAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGa                                                                                                               <  2:522876/138‑1 (MQ=255)
            ggatCGATGTATTATTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGAATTCATCGCACGGCAATAAAAGTACCCCATCGCGTAACAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGattat                                                                                                           <  1:409946/139‑1 (MQ=255)
              atCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCTCGTAAAACACCCATCATCTATATCACCCACAATCTGTTTTATTGTATTAGGattatga                                                                                                         >  2:41043/1‑139 (MQ=255)
              atCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGAttatta                                                                                                         >  1:794389/1‑141 (MQ=255)
                 gCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATg                                                                                                      <  2:459455/141‑1 (MQ=255)
                     tattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCCCATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGGATTAGGATTATTAATGCCat                                                                                                  >  1:657870/1‑141 (MQ=255)
                      attaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCAatt                                                                                                 >  1:631070/1‑141 (MQ=255)
                        tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATCTCACTCATAATATGTTTTATTGTATTAGGATTATTAATGACattat                                                                                               >  1:422325/1‑141 (MQ=255)
                        tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                               <  2:682454/141‑1 (MQ=255)
                        tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                               >  1:438396/1‑141 (MQ=255)
                        tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                               >  1:506954/1‑141 (MQ=255)
                            agCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                                                           <  1:58960/141‑1 (MQ=255)
                                cttGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                                                       <  1:27092/140‑1 (MQ=255)
                                cttGGCATGGGGTTTGTCTCTTGGATAAAGCGAACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                                                       <  1:28423/140‑1 (MQ=255)
                                aTTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACCGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                                                       <  1:201575/141‑1 (MQ=255)
                                    gCATGGGGTTTGTCTCTTCGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCAc                                                                                   <  1:313826/141‑1 (MQ=255)
                                     cATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                                                  <  1:34856/141‑1 (MQ=255)
                                      aTGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCACGTAAAAAAACCAGCATCTATATCCACTTATTGACTGTGTATTGTATTAGGATTATTAATGCCATAATTTGGATTATCACtt                                                                                 >  1:848159/1‑141 (MQ=255)
                                         gggTGTGTCTCTTGAATAAAGCTTACAAAAAATACTGCAATAAAACTACCAAATCTTATAAAAAAACCACCATCTATATCACTCATAATATGTTTTATTGTATTAAAATTATTAATGCCATTATTTTTATTATCACTTATc                                                                              <  2:170087/141‑1 (MQ=255)
                                                  tctTGGATAAAGCGTACAAAAAAAACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAttt                                                                     <  1:773666/141‑1 (MQ=255)
                                                   ctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATGATTAGTGCCATTATTTGGATTATCACTTATCCTTGTATTta                                                                    <  2:204018/141‑1 (MQ=255)
                                                    ttGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTtat                                                                   <  1:169338/141‑1 (MQ=255)
                                                     tGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATCATATGTTTTATTGTATTAGGATTATTCATGCCATTATTTGGAATATCACTTATCCTTGTATTtata                                                                  >  2:364615/1‑141 (MQ=255)
                                                           aaGCCTACAAAAAATACTACAACAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGaa                                                            <  2:731076/141‑1 (MQ=255)
                                                            aGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                                                           <  2:762941/141‑1 (MQ=255)
                                                             gcttACAAAAAATACTGCAACAAAAGTACCACATCGCGTAAAAAAACCAGCATCGATATCACTCATAATATTTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                                                          <  2:116527/138‑1 (MQ=255)
                                                               gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                                                        >  2:522975/1‑141 (MQ=255)
                                                                   aaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatt                                                    <  1:35418/141‑1 (MQ=255)
                                                                    aaaaaTACGGAAATAAAAGTACCACATCGCATAAAAAAACCAGAATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatta                                                   <  2:839312/141‑1 (MQ=255)
                                                                     aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat                                                  <  1:673211/141‑1 (MQ=255)
                                                                     aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat                                                  >  1:775633/1‑141 (MQ=255)
                                                                     aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTAACTTAAttttat                                                  >  2:141678/1‑141 (MQ=255)
                                                                      aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtata                                                 <  2:733478/141‑1 (MQ=255)
                                                                             gCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaaa                                          >  2:422203/1‑141 (MQ=255)
                                                                               caTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTATTGCCATTATTTGGATTTTCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGa                                        <  1:382628/140‑1 (MQ=255)
                                                                               aaTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGa                                        <  1:382634/141‑1 (MQ=255)
                                                                                aTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt                                       >  1:552562/1‑141 (MQ=255)
                                                                                aTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGGTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAATTATAATTGAATTAATATTATATATTAAAgct                                       >  1:579961/1‑139 (MQ=255)
                                                                                  aaaaGTACCACTTCGCGTAACAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGGATTTATAATTGAATTAATATTATATATTAAAGAtcg                                     <  2:547628/141‑1 (MQ=255)
                                                                                     aGTACCACGTCCGCTAAAAAAACCATCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgtcg                                  <  2:785988/141‑1 (MQ=255)
                                                                                     aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTTAATTAATATTATATATTAAAGAtcgtct                                  >  1:256057/1‑140 (MQ=255)
                                                                                          cacaTCGCGTAAAAATACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTa                             >  1:645019/1‑141 (MQ=255)
                                                                                              tCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTAGTGCTAAACa                         >  1:690260/1‑141 (MQ=255)
                                                                                                 gttAAAAAGACCAGCGTTTTTATCGATCACATTATGTTTTATCGTATTAGGATTATTAATGCCATTATCTGGATTATCACTTATCCTTGGATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaata                      <  2:300861/139‑1 (MQ=255)
                                                                                                   tAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaataat                    <  2:422325/141‑1 (MQ=255)
                                                                                                   tAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaataat                    <  2:438396/141‑1 (MQ=255)
                                                                                                      aaaaaCCAGCATCTATATCACTCCTAATTTTTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCa                 <  2:389612/141‑1 (MQ=255)
                                                                                                      aaaaaCAAGCATCTATATCACTCATAATATGGTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTCTCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCa                 <  2:389568/141‑1 (MQ=255)
                                                                                                        aaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt               <  2:84287/141‑1 (MQ=255)
                                                                                                          aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTa             >  2:801527/1‑141 (MQ=255)
                                                                                                                 tCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGtttt      >  2:409747/1‑141 (MQ=255)
                                                                                                                    atatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGaa   >  2:324171/1‑141 (MQ=255)
                                                                                                                     tatCACTCATAATATGTTTTATTGGATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAAc  >  1:550064/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                         
ATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAAC  >  CP000730/1513794‑1514051

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: