Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,756,899 T→C 7.5% intergenic (‑1/+99) USA300HOU_2612 ← / ← nrdG hypothetical protein/anaerobic ribonucleotide reductase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,756,8990TC 7.5% 152.1 / 3.3 53intergenic (‑1/+99)USA300HOU_2612/nrdGhypothetical protein/anaerobic ribonucleotide reductase small subunit
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base C (2/2);  total (23/30)
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.26e-01

TTTATTTATTGTAAAACGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGA  >  CP000730/2756759‑2757039
                                                                                                                                            |                                                                                                                                            
tttatttattGTAAAACGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTCTTGGTATCAAAATGATTGTTAGGTCATATTCTATCAat                                                                                                                                              >  2:218273/1‑141 (MQ=255)
  tatttattGTAAAACGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTATCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAatat                                                                                                                                            >  2:278767/1‑141 (MQ=255)
  tatttattGTAAAACGCTACCAAAAAGTGCGACGTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATCCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATAATCTCTCAatat                                                                                                                                            >  2:232122/1‑141 (MQ=255)
    tttattGTAAAACGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAatatat                                                                                                                                          >  1:114525/1‑141 (MQ=255)
       attGTAAAACGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATAtttt                                                                                                                                       <  2:356181/141‑1 (MQ=255)
             aaaCGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATaaa                                                                                                                                 >  1:143931/1‑141 (MQ=255)
             aaaCGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATaaa                                                                                                                                 <  1:185565/141‑1 (MQ=255)
              aaCGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCCCTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAg                                                                                                                                <  1:174321/141‑1 (MQ=255)
              aaCGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAg                                                                                                                                <  1:174363/141‑1 (MQ=255)
              aaCGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAg                                                                                                                                <  2:404383/141‑1 (MQ=255)
              aaCGCTACCAAAAAGTTAGACTTCCGTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAg                                                                                                                                <  1:305000/141‑1 (MQ=255)
                 gCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTTTAAAGAAt                                                                                                                             >  2:12187/1‑141 (MQ=255)
                  cTACCAAAAAGTTAGATTTCCTTCCCACTAAAATACCACTTTTATCTTTCAAATTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAGAGAAtt                                                                                                                            >  1:228350/1‑141 (MQ=255)
                  cTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAAtt                                                                                                                            <  1:272965/141‑1 (MQ=255)
                         aaaGTTAGACTTCCTTCCCTCTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTtta                                                                                                                     <  1:295080/141‑1 (MQ=255)
                          aaGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTttat                                                                                                                    >  1:445679/1‑141 (MQ=255)
                                 aCTTCCTTCCCACTAAAATTCCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAGATTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACttt                                                                                                             <  2:331559/141‑1 (MQ=255)
                                          ccACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGt                                                                                                    <  1:64874/141‑1 (MQ=255)
                                              tAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTacct                                                                                                >  1:19735/1‑141 (MQ=255)
                                               aaaaTACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTaccta                                                                                               >  1:25974/1‑141 (MQ=255)
                                                  aTACCCCCTGTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTacctaacc                                                                                            <  1:75487/141‑1 (MQ=255)
                                                    aCCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTacctaaccta                                                                                          <  1:138279/141‑1 (MQ=255)
                                                       aCTTTTCTCTTTCAACTTTTTTAAAAAACGCATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTACCATTAAATTATGTACCTAACCTacaa                                                                                       >  1:223748/1‑141 (MQ=255)
                                                        cTTTTCTCTTTCAACTGTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTaaaaa                                                                                      <  1:171523/141‑1 (MQ=255)
                                                            tctcTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAACATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAaagaag                                                                                  >  1:28999/1‑141 (MQ=255)
                                                            tctcTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAACATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAaagaag                                                                                  >  1:73681/1‑141 (MQ=255)
                                                              tctTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGcc                                                                                >  1:301728/1‑141 (MQ=255)
                                                               ctTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCa                                                                               >  1:311419/1‑141 (MQ=255)
                                                                         ttttAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGaa                                                                     <  1:372474/141‑1 (MQ=255)
                                                                         ttatAAAAAAAGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGAAAAAAAAAAAAAAAGAAGCCAAGGCAACGaa                                                                     <  2:291140/138‑1 (MQ=255)
                                                                         tataaaacaaCGGATATTCAACCTTTAGTAATGGTATCAAAATGCTTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAACTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGaa                                                                     <  1:217494/137‑1 (MQ=255)
                                                                            tAAAAAACGGATATGCCACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGt                                                                  <  2:6791/141‑1 (MQ=255)
                                                                            tAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGt                                                                  <  1:77109/141‑1 (MQ=255)
                                                                                  aCGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGTCAAGGTAAcattcgttgcctt                                                            <  1:151450/141‑13 (MQ=255)
                                                                                     gATATGCAACTTTTAGTATTGGTATCAAAATAATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAACGAAGCCAAAGCCACGAATGTTACCttcgc                                                         >  1:423963/1‑138 (MQ=255)
                                                                                       tatGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACtt                                                       >  2:1893/1‑141 (MQ=255)
                                                                                         tGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGCCTTCt                                                     >  1:103256/1‑141 (MQ=255)
                                                                                                 ttAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACata                                             <  2:397268/141‑1 (MQ=255)
                                                                                                 ttAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAACATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACata                                             <  2:73681/141‑1 (MQ=255)
                                                                                                 ttAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAACATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACata                                             <  2:28999/141‑1 (MQ=255)
                                                                                                 atatcagtGGTATCAAAATGATTGATAGGTCATATTCTATCAATATATTTTTACAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACata                                             <  2:451819/134‑1 (MQ=255)
                                                                                                      aTTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCaa                                        <  1:186807/141‑1 (MQ=255)
                                                                                                        tGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCaact                                      >  1:453052/1‑141 (MQ=255)
                                                                                                            aTCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCaactaact                                  <  1:299424/141‑1 (MQ=255)
                                                                                                              cAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACtat                                <  1:190561/141‑1 (MQ=255)
                                                                                                              cAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACtat                                <  1:191979/141‑1 (MQ=255)
                                                                                                                aaaTGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACtatat                              >  2:393154/1‑141 (MQ=255)
                                                                                                                  aTGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTACCTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAactaaatacttat                            >  2:296850/1‑136 (MQ=255)
                                                                                                                    gATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTc                          >  2:453124/1‑141 (MQ=255)
                                                                                                                         ttAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACGTATTCAACTAACTATATATTCAATCa                     >  2:263512/1‑141 (MQ=255)
                                                                                                                          tAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCAt                    >  2:193337/1‑141 (MQ=255)
                                                                                                                              tCATATTCTATCAATATATTTTTATAAAGAATGGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATAcgc                <  2:114525/141‑1 (MQ=255)
                                                                                                                               cATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATAcgcg               >  1:420664/1‑141 (MQ=255)
                                                                                                                                  atTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCAt            >  1:408703/1‑141 (MQ=255)
                                                                                                                                      tATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCga        >  2:69456/1‑141 (MQ=255)
                                                                                                                                        tCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCgaga      <  2:302802/141‑1 (MQ=255)
                                                                                                                                         cAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGagag     <  2:312814/141‑1 (MQ=255)
                                                                                                                                            tatatTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGa  >  1:398683/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTTATTTATTGTAAAACGCTACCAAAAAGTTAGACTTCCTTCCCACTAAAATACCACTTTTCTCTTTCAACTTTTTTAAAAAACGGATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGA  >  CP000730/2756759‑2757039

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: ‑