Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,243,494 A→T 8.9% intergenic (+72/+37) USA300HOU_2114 → / ← pyrG hypothetical protein/CTP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,243,4940AT 8.9% 171.2 / 7.9 57intergenic (+72/+37)USA300HOU_2114/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base A (18/33);  new base T (2/3);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

ATTGAATTATGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCA  >  CP000730/2243356‑2243633
                                                                                                                                          |                                                                                                                                           
aTTGAATTATGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTATCATAATATTAAAAAAAAAAAAAAAACCGGTACTTTCAATg                                                                                                                                           <  2:175281/141‑1 (MQ=255)
  ggAATTAGGTATAGTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGtt                                                                                                                                         <  2:759061/140‑1 (MQ=255)
        aTGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                                                                   <  1:669455/141‑1 (MQ=255)
         tGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACatt                                                                                                                                  <  2:254513/141‑3 (MQ=255)
           tataCTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                                                                <  1:2041/141‑1 (MQ=255)
             taCTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                                                              >  2:222490/1‑141 (MQ=255)
               cTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGt                                                                                                                            >  2:482049/1‑141 (MQ=255)
                tATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                                                                           >  2:763960/1‑141 (MQ=255)
                   gACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                                                                        <  1:161249/141‑1 (MQ=255)
                     cATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGt                                                                                                                      >  2:366576/1‑141 (MQ=255)
                         tataCACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTtta                                                                                                                  >  1:572303/1‑141 (MQ=255)
                           taCACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatt                                                                                                                >  1:241085/1‑141 (MQ=255)
                            acacAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                                                                   >  2:758140/1‑120 (MQ=255)
                            acacAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                                                                   <  1:758140/121‑2 (MQ=255)
                              acaAATGATTGATATGACAATAGACTTAAAAATATAAAAGTCACTATTGTTATTTCATCCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTATATTAGc                                                                                                             <  2:545852/141‑1 (MQ=255)
                              acaAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACAGGTTTTATTAGc                                                                                                             >  1:768023/1‑141 (MQ=255)
                                   taaTTGAAACTACAAGAAACTTATAATGATAAAATTCACTATTATTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                                                                        <  2:467993/139‑1 (MQ=255)
                                    gATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAattt                                                                                                       <  2:696420/141‑1 (MQ=255)
                                      ttGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                                                                     <  2:547380/141‑1 (MQ=255)
                                      ttGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                                                                     <  1:52148/141‑1 (MQ=255)
                                       tGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatt                                                                                                    <  1:192936/141‑1 (MQ=255)
                                          aaTGATAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttat                                                                                                 >  2:165823/1‑141 (MQ=255)
                                          aaTGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGTAAatttatttat                                                                                                 >  2:224730/1‑141 (MQ=255)
                                                aaGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgtt                                                                                           <  1:630452/141‑1 (MQ=255)
                                                   agaCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGat                                                                                        <  1:567530/141‑1 (MQ=255)
                                                     aCTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGatat                                                                                      <  1:361025/141‑1 (MQ=255)
                                                         aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATCTTAAAATTATACGAAAAACCGGAACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTACTTTGTTGATATTTTa                                                                                  >  2:755566/1‑141 (MQ=255)
                                                         aGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                                                                  <  2:249556/141‑1 (MQ=255)
                                                            aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATg                                                                               <  2:242754/141‑1 (MQ=255)
                                                                taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                                                                           <  2:718837/141‑1 (MQ=255)
                                                                 aaaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCt                                                                          >  1:45606/1‑141 (MQ=255)
                                                                         cTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGCACTTTCACTGTTTAACAATGAAAGTACAGGCATTATTAACAAGTATATTTATTTTGTTGATAATTTAATGAAGCTTCAATaaa                                                                  >  1:508521/1‑141 (MQ=255)
                                                                                 aTTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaa                                                          >  2:204548/1‑141 (MQ=255)
                                                                                      aTACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc                                                     <  2:500513/141‑1 (MQ=255)
                                                                                        aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg                                                   >  1:152015/1‑141 (MQ=255)
                                                                                              tCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCgg                                             <  1:96442/141‑1 (MQ=255)
                                                                                               cACTAACATAATATTAAAATTATACGAACAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATAAATTAAAAATCGGGTGCGGa                                            <  2:393081/141‑1 (MQ=255)
                                                                                                  tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGa                                         <  1:595806/141‑1 (MQ=255)
                                                                                                  tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGa                                         <  2:241085/141‑1 (MQ=255)
                                                                                                    aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGAtt                                       <  2:301437/141‑1 (MQ=255)
                                                                                                       taataTTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTgg                                    <  2:145176/141‑1 (MQ=255)
                                                                                                           atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGtata                                >  1:166588/1‑141 (MQ=255)
                                                                                                                 aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAg                          <  1:251071/141‑1 (MQ=255)
                                                                                                                  ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGa                         <  2:302018/141‑1 (MQ=255)
                                                                                                                       cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCt                    >  2:98747/1‑141 (MQ=255)
                                                                                                                        gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTg                   >  1:682917/1‑141 (MQ=255)
                                                                                                                          aaaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGTAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTggg                 <  1:224730/141‑1 (MQ=255)
                                                                                                                              ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                                                                       <  2:260916/83‑1 (MQ=255)
                                                                                                                              ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                                                                       >  1:260916/1‑83 (MQ=255)
                                                                                                                              ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGa             <  2:681851/141‑1 (MQ=255)
                                                                                                                                 gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAAtt          >  1:396670/1‑141 (MQ=255)
                                                                                                                                    cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGAc       <  2:56888/141‑1 (MQ=255)
                                                                                                                                      ttCATTGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACaa     >  2:288026/1‑141 (MQ=255)
                                                                                                                                       tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAg    <  2:103745/141‑1 (MQ=255)
                                                                                                                                       tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAg    <  2:532565/141‑1 (MQ=255)
                                                                                                                                        cAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGc   <  2:768023/141‑1 (MQ=255)
                                                                                                                                        cAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGc   <  2:545849/141‑1 (MQ=255)
                                                                                                                                         aaTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCa  >  1:619178/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
ATTGAATTATGTATACTATGACATTTATACACAAATGATTGAAATGACAAGAGACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCGGACGATTTGGTCTAGATAAGAATTCTGGGTGGAATTGACAAGCA  >  CP000730/2243356‑2243633

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