Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 195,286 T→A intergenic (‑97/‑82) USA300HOU_0183 ← / → USA300HOU_0184 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730195,2860TA95.7% 71.1 / NA 23intergenic (‑97/‑82)USA300HOU_0183/USA300HOU_0184hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (7/15);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGA  >  CP000730/195148‑195416
                                                                                                                                          |                                                                                                                                  
tttGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTaaa                                                                                                                                  >  2:79350/1‑141 (MQ=255)
  tGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAaaaa                                                                                                                                >  2:264420/1‑141 (MQ=255)
    tttCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAaaaaaa                                                                                                                              <  2:175976/141‑1 (MQ=255)
     tttttttttaaCTAGATAAGAGAGAATCTTTATCCATATTGATACCGGCCTATATATATGTTTCTTTACTTATTCCCTTTTTTGTTTTTTATTCTCTTTGAAAATATATTTTTCAGTATTTTTTCTACATATTTaaaacaa                                                                                                                             >  2:417106/10‑141 (MQ=255)
                                         tAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATtatcg                                                                                         >  2:49176/1‑139 (MQ=255)
                                           aTGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTa                                                                                       >  2:145327/1‑141 (MQ=255)
                                                   ccccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATAttttt                                                                               <  2:418777/141‑1 (MQ=255)
                                                      cTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCAAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCa                                                                            <  2:329890/141‑1 (MQ=255)
                                                             tatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTg                                                                     >  2:214220/1‑141 (MQ=255)
                                                                      ttatttatACTCTAACATTTTTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATCTCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATaggag                                                            <  2:399900/141‑1 (MQ=255)
                                                                       ctctcataCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATaggagg                                                           <  2:283870/136‑1 (MQ=255)
                                                                        aCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGt                                                          >  2:239980/1‑141 (MQ=255)
                                                                                 cTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAAt                                                 <  2:410675/141‑1 (MQ=255)
                                                                                   aaCATTATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGa                                               <  2:54808/141‑1 (MQ=255)
                                                                                      aTGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACat                                            >  2:322348/1‑141 (MQ=255)
                                                                                        gATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACatta                                          <  2:289252/141‑1 (MQ=255)
                                                                                        caTTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATAAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCTTTGTTATAGGAGGTCTTATTAATGATatta                                          <  2:301929/140‑1 (MQ=255)
                                                                                             ttATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTAttttt                                     <  2:343146/141‑1 (MQ=255)
                                                                                                        gAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTa                          <  2:341695/141‑1 (MQ=255)
                                                                                                                   ttACAGAATTTTAAAAAAAAAAAAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAt               <  2:293271/141‑1 (MQ=255)
                                                                                                                    tACAGAATTTTATCTAAAAATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAtt              <  2:1559/141‑1 (MQ=255)
                                                                                                                          aTTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCAt        <  2:444892/141‑1 (MQ=255)
                                                                                                                                tCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa  <  2:291060/141‑1 (MQ=255)
                                                                                                                                tCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAACCTAACAATCTTAATTTTGCATcgacga  >  2:260780/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                  
TTTGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGA  >  CP000730/195148‑195416

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: