Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 911,610 G→A intergenic (‑221/‑107) USA300HOU_0869 ← / → sufC hypothetical membrane protein/ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730911,6100GA97.0% 106.0 / ‑3.7 33intergenic (‑221/‑107)USA300HOU_0869/sufChypothetical membrane protein/ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (17/15);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAGTAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGAAATCAAAG  >  CP000730/911475‑911741
                                                                                                                                       |                                                                                                                                   
tttaaTTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATaaataagt                                                                                                                                <  2:317779/138‑1 (MQ=255)
ttAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATaaataagt                                                                                                                                <  1:339471/141‑1 (MQ=255)
  aaaTTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTta                                                                                                                              <  1:314508/141‑1 (MQ=255)
      tCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTGTATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAAt                                                                                                                          <  2:95540/141‑1 (MQ=255)
               ccATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAAc                                                                                                                 >  1:100836/1‑141 (MQ=255)
               ccATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAAc                                                                                                                 >  1:113681/1‑141 (MQ=255)
                    gagaAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCatt                                                                                                            <  2:353560/141‑1 (MQ=255)
                      gaAATTATGTGCTCTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCattat                                                                                                          <  2:370206/141‑1 (MQ=255)
                          ttATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACa                                                                                                      >  2:198055/1‑141 (MQ=255)
                           tATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACaa                                                                                                     <  2:156601/141‑1 (MQ=255)
                                      attatCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGa                                                                                          <  1:177402/141‑1 (MQ=255)
                                             tttATATTTTTAAAGAGAGCGGATTAGAGTTGGACGTTAGGGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACaaa                                                                                   >  2:200679/1‑141 (MQ=255)
                                                              agCGGATTAGAGTTGTACTTTAAAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACACATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAAc                                                                  <  2:193477/141‑1 (MQ=255)
                                                                     tAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTaa                                                           <  2:327597/141‑1 (MQ=255)
                                                                        agTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGaa                                                        >  1:50143/1‑141 (MQ=255)
                                                                                         attaAAGATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAAGTTTATATTTTACgg                                       >  2:303663/1‑141 (MQ=255)
                                                                                         attaAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAGTAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACgg                                       >  1:45543/1‑141 (MQ=255)
                                                                                         attaAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACgg                                       <  1:350013/141‑1 (MQ=255)
                                                                                         attaAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACgg                                       >  1:214792/1‑141 (MQ=255)
                                                                                              aaaTAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGa                                  >  1:289651/1‑141 (MQ=255)
                                                                                                ataataAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTAGAAAGTTTTTTCGACCCCAATTTAATATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAat                                >  1:34311/1‑141 (MQ=255)
                                                                                                    taaAAAGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATaa                            >  2:222811/1‑141 (MQ=255)
                                                                                                        aaGGGTTTAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATg                        <  2:138463/141‑1 (MQ=255)
                                                                                                             tttAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATc                   >  1:277974/1‑141 (MQ=255)
                                                                                                             tttAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATc                   <  1:200679/141‑1 (MQ=255)
                                                                                                             tttAAGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATc                   <  1:13695/141‑1 (MQ=255)
                                                                                                                aaGTTGTTTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAAc                >  1:8689/1‑141 (MQ=255)
                                                                                                                  gttgttTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACAt              >  1:206344/1‑141 (MQ=255)
                                                                                                                    tgttTATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTa            <  1:98902/141‑1 (MQ=255)
                                                                                                                      tttATATTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGa          >  2:332104/1‑141 (MQ=255)
                                                                                                                      ttcatacTAAGGTATAAATAAGTTAAAATTAACTGAACGCATTATTACAAAGTCTTTTTGAGTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTAATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGa          <  2:194666/138‑1 (MQ=255)
                                                                                                                         atatTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGAAAt       >  2:168607/1‑141 (MQ=255)
                                                                                                                         atagTACGGTATAAATAACTTATAATTAACTGCACGTATTATTACAAAGTCTTTTTGGCTATATATTAAAATTATTCTAAACAAGTTAATAAACCTTTATATTTTACGGAGAGCATATAAAATGGCATCAACATTAAAAAt       <  2:169609/141‑1 (MQ=255)
                                                                                                                           atTAAGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGAAATCa     >  1:339188/1‑141 (MQ=255)
                                                                                                                              aaGGTATAAATAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGAAATCAAAg  >  2:345769/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                   
TTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAAGAGAGCGGATTAGAGTTGTACTTTAGAGTTATTAAAAATAATAAAAAGGGTTTAAGTTGTTTATATTAAGGTATAAGTAAGTTATAATTAACTGAACGCATTATTACAAAGTCTTTTTGACTACAAATTAAAATTATTATAAACTAGTTAAGAAAACTTTATATTTTACGGAGGGAATATAAAATGGCATCAACATTAGAAATCAAAG  >  CP000730/911475‑911741

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