Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,419,527 G→T P7T (CCT→ACT)  USA300HOU_1317 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,419,5270GT96.9% 92.1 / ‑2.9 32P7T (CCT→ACT) USA300HOU_1317hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (25/6);  total (25/7)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TTGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGGCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGAAAAATCATT  >  CP000730/1419395‑1419663
                                                                                                                                    |                                                                                                                                        
ttGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGa                                                                                                                                  >  2:167712/1‑141 (MQ=255)
 tGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaa                                                                                                                                 >  2:24424/1‑141 (MQ=255)
   gAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaata                                                                                                                               <  1:294190/141‑1 (MQ=255)
               cAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATattca                                                                                                                   >  2:47924/1‑141 (MQ=255)
                             ataACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTAttt                                                                                                     >  2:325925/1‑141 (MQ=255)
                                 cAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCAtt                                                                                                 >  2:367482/1‑141 (MQ=255)
                                 cAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGGCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCAtt                                                                                                 <  1:208393/141‑1 (MQ=255)
                                      gTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACATAATATTTTTAATAAAATCATTATTAATGTTTTGATCCATAAATAATTAGTCAGATGCAATAATGTTATATTCATTTCCCTTATATTTCATTACttt                                                                                            >  2:155334/1‑141 (MQ=255)
                                            aCTTCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATTGAATAATGTCATATTCATTATCATTTTATTTCATTACTTTCTTTGt                                                                                      >  2:249127/1‑141 (MQ=255)
                                             cTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGta                                                                                     <  1:170978/141‑1 (MQ=255)
                                             cTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGta                                                                                     >  2:59400/1‑141 (MQ=255)
                                                        aaCCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTttatt                                                                          >  2:143634/1‑141 (MQ=255)
                                                        aaCCATGTATTTAAAACTAAGCAACAGAATATTGGTGATAAGATCATTATTAATGTGGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTAATTTCATTACGTTGGTTTTATACGTTttatt                                                                          >  2:34521/1‑141 (MQ=255)
                                                                 tttaAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAAtttt                                                                 >  2:249731/1‑141 (MQ=255)
                                                                 tttaAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAAtttt                                                                 >  2:139646/1‑141 (MQ=255)
                                                                  ttaAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATCATGACATTTTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAAttttt                                                                >  1:190569/1‑141 (MQ=255)
                                                                        cTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAAc                                                          <  2:80353/141‑1 (MQ=255)
                                                                                        tGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGcagtc                                          <  2:129599/141‑1 (MQ=255)
                                                                                         gTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGcagtcg                                         >  2:301904/1‑140 (MQ=255)
                                                                                            aataaAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGtt                                      >  1:23508/1‑141 (MQ=255)
                                                                                                 aaTCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAg                                 >  2:180328/1‑141 (MQ=255)
                                                                                                    cATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCaa                              >  2:186767/1‑141 (MQ=255)
                                                                                                            aTGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATtctat                      >  2:274651/1‑141 (MQ=255)
                                                                                                             tgttgtGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATtctatc                     >  1:107116/1‑141 (MQ=255)
                                                                                                                tgtgATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATctac                  >  2:103849/1‑141 (MQ=255)
                                                                                                                 gtgATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATctact                 >  1:271239/1‑141 (MQ=255)
                                                                                                                 gtgATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATctact                 >  2:66225/1‑141 (MQ=255)
                                                                                                                  tgATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATctacta                >  1:216401/1‑141 (MQ=255)
                                                                                                                       cATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGa           >  1:222021/1‑141 (MQ=255)
                                                                                                                        ataaataaTTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGAATACATTTTATTATAAAGTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGaa          >  2:278400/1‑141 (MQ=255)
                                                                                                                           aataatTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGaaaaa       <  1:183296/141‑1 (MQ=255)
                                                                                                                                ttAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGAAAAATCAtt  <  1:167712/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                        
TTGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGGCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGAAAAATCATT  >  CP000730/1419395‑1419663

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