Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000731 25,669 (A)7→6 coding (322/333 nt) USA300HOU_pUSA300HOUMR0031 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073125,6630A.97.9% 213.7 / ‑1.9 49coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (30/17);  minor base T (1/0);  total (31/18)
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

CACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTA  >  CP000731/25538‑25782
                                                                                                                             |                                                                                                                       
cacTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTTAAGGATCTTGGTTTTCTAAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                                                      >  1:65602/1‑129 (MQ=255)
 acTAGTGCTTAAATTCACAGGTTTGTCACCAAGTATTGCAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                                                     <  2:12720/129‑1 (MQ=255)
  cTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATacaaaa                                                                                                                    >  2:116828/1‑129 (MQ=255)
     gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTATCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                                                 >  1:103976/1‑126 (MQ=255)
     gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTAATTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                                                 >  1:186114/1‑126 (MQ=255)
         ttAAATTTACAGGTTTGTCACCAAGTATTGAAGGCTCTTGCTTATCTCAATATTTAGATATCTTCTATATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                                                            >  1:172965/1‑129 (MQ=255)
         ttAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATACAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACAt                                                                                                            >  1:125609/1‑129 (MQ=255)
                   aGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGATTa                                                                                                  >  2:108077/1‑129 (MQ=255)
                       ttGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTATCACTATTTCTTCTGATGTTTAATTACTATACTACAATTCAAAAATATACATCTAAACGTTCAAtt                                                                                              >  2:153996/1‑129 (MQ=255)
                       ttGTCACCAAGTATTCAAGGATCTTGGTTTTATCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAAtt                                                                                              <  1:228445/129‑1 (MQ=255)
                       ttGTCAACAAGTATTTAAGGATCTTTGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAAtt                                                                                              <  2:25056/129‑1 (MQ=255)
                        tGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTa                                                                                             <  1:3246/129‑1 (MQ=255)
                        tGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTa                                                                                             <  1:131004/129‑1 (MQ=255)
                            actaaGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                                                         <  2:147764/126‑1 (MQ=255)
                            aCCAAGTATTGAAAGATCTTCGTTTTCTTAATATCCAGAAATTTTTTTTTTAATTGTATTACTATTTCTTATGCTGTGTAAGTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGa                                                                                         <  1:68673/129‑1 (MQ=255)
                                   aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                                                  >  1:10125/1‑129 (MQ=255)
                                    ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTAATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                 >  1:112731/1‑129 (MQ=255)
                                        aGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATATAAAGGTTAAATTAATGAAAAATAGAATTa                                                                             >  2:6896/1‑129 (MQ=255)
                                           aTCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATAagag                                                                          >  2:152357/1‑129 (MQ=255)
                                            tCTTGGTTTTCTCAATATCCAGAAATTTTCTTGATAATTGTAGCACTATTTCTTATGCTGTTTAATTATTATACTACAAATAAAAAATCTCCTTCTTTAGGTTAATTTAATGAAAAATAGAAttctagt                                                                         >  2:82411/1‑124 (MQ=255)
                                            tCTTGGTTTTCTCAATATCCAGAAATCTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTACAGCTTAAACTAATTAAAACCAGCATTaccga                                                                         >  2:137274/1‑125 (MQ=255)
                                                 gTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATTCTGTTTAATTACTATATTAAAAATAAAAAATAAACAAATAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAtata                                                                    <  1:106088/129‑1 (MQ=255)
                                                  ttttCTCAATATCCAGAAATTTTCTTTATAATTATAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                                                   <  1:39438/129‑1 (MQ=255)
                                                    ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                                                 >  1:12426/1‑129 (MQ=255)
                                                     tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTTATGACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                                                >  2:115363/1‑129 (MQ=255)
                                                     tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                                                >  2:76455/1‑129 (MQ=255)
                                                     tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                                                >  1:155470/1‑129 (MQ=255)
                                                     tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                                                >  1:188271/1‑129 (MQ=255)
                                                     tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATACTGATTAATTACAATACAACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                                                >  1:38998/1‑129 (MQ=255)
                                                       tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAc                                                              >  1:4549/1‑129 (MQ=255)
                                                            atCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                                                         >  2:34309/1‑129 (MQ=255)
                                                              ccAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                                                       <  2:183379/129‑1 (MQ=255)
                                                               cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACACAAAGGa                                                      >  2:118460/1‑129 (MQ=255)
                                                               cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                                                      <  2:67195/129‑1 (MQ=255)
                                                                    aTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                                                 >  2:33878/1‑129 (MQ=255)
                                                                    aTTTTCTGTATAAGTGTAGTACTATTACTTATTCTTTATAAATAATATATTACAAATAAAAAAGATACATCTAAAGATTAAATTAATTAAAAATAGAATTAGAGAATATATAAAACAAAAAGGAATTtc                                                 <  2:203184/129‑1 (MQ=255)
                                                                        tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTAAAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAGTATAGAAAACAAAAAGGAATTTCTCaa                                             <  2:209868/129‑1 (MQ=255)
                                                                         cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTATAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg                                            <  2:228836/129‑1 (MQ=255)
                                                                          tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGTTTAAATTAATGAAAAATAGAATTACAGAATATAGAAAACAAAAAGAAATTTATAAAGa                                           >  2:185191/1‑129 (MQ=255)
                                                                            tataATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaaa                                         >  1:137793/1‑129 (MQ=255)
                                                                              taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                       >  2:97179/1‑129 (MQ=255)
                                                                              taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt                                       >  1:80159/1‑129 (MQ=255)
                                                                               aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                      >  2:134713/1‑129 (MQ=255)
                                                                               aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                      >  2:7183/1‑129 (MQ=255)
                                                                                   gTAGGACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAg                                  <  1:213544/129‑1 (MQ=255)
                                                                                     aGCACTATTTCTTATGCTGTTTAATAACTATAATACAAATAAAAAATATACATGTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                <  2:43675/129‑1 (MQ=255)
                                                                                        aCTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAg                             >  1:49502/1‑129 (MQ=255)
                                                                                              tCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAAAATTAAAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaa                       >  1:172536/1‑129 (MQ=255)
                                                                                               ctaaTGCTGGTTAATTACGATACTACAAATAAAAAATATACATCTAAAAGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGTAATTTCTCAAGAAACTTTAAGTAAAGAATTaaa                      <  2:234715/126‑1 (MQ=255)
                                                                                                ttATGCTTTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATTAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaaa                     <  1:22617/129‑1 (MQ=255)
                                                                                                     cTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAATTTAAGTAAAGAATTAAAAGTGTc                >  1:97884/1‑129 (MQ=255)
                                                                                                                  atactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTAtt   <  1:51620/129‑1 (MQ=255)
                                                                                                                   tactacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTa  >  2:144587/1‑129 (MQ=255)
                                                                                                                   tactacAAATAAAAAATATACATCTAAAAGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAAAAAAAAAAATTTATCAAAAAACTTTAAGTAAAGAATTAAAAGTATAAAGAAAAACTATTa  >  2:216981/1‑129 (MQ=255)
                                                                                                                             |                                                                                                                       
CACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTA  >  CP000731/25538‑25782

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