Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,078,566:1 +T 100% coding (253/681 nt) USA300HOU_1935 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,078,5661.T92.0% 71.9 / 0.7 25coding (253/681 nt)USA300HOU_1935hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (16/7);  total (18/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GATATGACAAATGCATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTG‑‑TTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAAATAATCAATATTGTTGGACTTATAAAAGA  >  CP000730/2078427‑2078692
                                                                                                                                              |                                                                                                                               
gATATGACAAATGCATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTG‑‑t                                                                                                                               <  1:57126/141‑1 (MQ=255)
  tatGACAAATGCATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTG‑‑ttt                                                                                                                             >  2:92524/1‑141 (MQ=255)
  tatGACAAATGCATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAAAATAGTTGACATTAAGCAATTG‑‑ttt                                                                                                                             >  2:257927/1‑141 (MQ=255)
  tatGACAAATGCATATGCAGTTGAAAAATAGATTGAATTTGGTTAACTCACATCTGCTTAGTAAGAATATAAAAAAATCACACGCGTTCCCATTAATGTTAAGACAATGAATGTAATATAGTTGGCATTAAGCAATTG‑‑ttt                                                                                                                             >  2:77129/1‑141 (MQ=255)
            gCATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTG‑‑ttttttttgaacc                                                                                                                   >  2:31719/1‑135 (MQ=255)
             cATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACt                                                                                                                   >  1:249602/1‑141 (MQ=255)
               tatGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTgg                                                                                                                 <  2:250935/141‑1 (MQ=255)
                 tGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTa                                                                                                               >  1:83674/1‑141 (MQ=255)
                    aGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAAc                                                                                                            >  2:181770/1‑141 (MQ=255)
                         aaaaaTAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGAATTAAGCGATGTT‑‑TTTTTTTGACACTGGGaaagtaatt                                                                                                       >  1:335009/1‑135 (MQ=255)
                          aaaaTAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGTTAAAATCataa                                                                                                      >  2:340619/1‑140 (MQ=255)
                          aaaaTAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCatat                                                                                                      >  1:273559/1‑141 (MQ=255)
                          aaaaTAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCatat                                                                                                      >  1:357153/1‑141 (MQ=255)
                          aaaaTAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGAAACTGGTAAACTCatat                                                                                                      >  1:148242/1‑141 (MQ=255)
                                      aTTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATg                                                                                          >  1:156569/1‑141 (MQ=255)
                                                         gCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCt                                                                       >  1:395702/1‑141 (MQ=255)
                                                              aTAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGAt                                                                  >  1:253119/1‑141 (MQ=255)
                                                                    atatAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTc                                                            <  2:389767/141‑1 (MQ=255)
                                                                    atatAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTc                                                            <  1:18264/141‑1 (MQ=255)
                                                                    aaaaaatttaataaCAACCGTTCCAATTAAGGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTc                                                            <  1:244649/133‑1 (MQ=255)
                                                                    aaaaaaaaTAATAAAAAGCATTCCAATTAAACTTAAAACAATGCATGTAATATAGTTAGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTc                                                            <  1:18563/137‑1 (MQ=255)
                                                                                  cAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCa                                              <  1:232519/141‑1 (MQ=255)
                                                                                     gCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATagagtacgagca                                           >  2:275750/1‑130 (MQ=255)
                                                                                       gTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATc                                         >  1:148010/1‑141 (MQ=255)
                                                                                                  tGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCaaa                              >  1:351179/1‑141 (MQ=255)
                                                                                                    ttAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAaata                            <  1:181770/141‑1 (MQ=255)
                                                                                                               cATGTAATATAGTTGGCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAAATAATCAATAttgt                 >  1:8840/1‑141 (MQ=255)
                                                                                                                         aGTTGGCATTAAGCCATTGG‑‑TTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAAATAATCAATATTGTTGGACTtata       >  1:217773/1‑141 (MQ=255)
                                                                                                                             ggCATTAAGCAATTGTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAAATAATCAATATTGTTGGACTTATAAAAg   >  1:9146/1‑141 (MQ=255)
                                                                                                                               cATTAAGCAATTGTTTTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAAATAATCAATATTGTTGGACTTATAAAAGa  >  1:69530/1‑141 (MQ=255)
                                                                                                                                              |                                                                                                                               
GATATGACAAATGCATATGCAGTTGAAAAATAGATTGAATTTGGTTCAATCACATCTGCTTCATAAGCATATAAACTAATAACAAGCGTTCCAATTAATGTTAAAACAATGCATGTAATATAGTTGGCATTAAGCAATTG‑‑TTTTTTTGACACTGGTAAACTCATATAAAAATAATATGAATCATTACCACCTAATCTACGATTGAGTCTAAATAAGTGACCAGCATCTAAAATCCAAATAATCAATATTGTTGGACTTATAAAAGA  >  CP000730/2078427‑2078692

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