Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,503,491 C→G 32.3% T672R (ACG→AGG)  USA300HOU_2364 → AraC transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,503,4910CG32.3% 33.4 / 31.5 38T672R (ACG→AGG) USA300HOU_2364AraC transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (7/18);  new base G (4/8);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CCTCAAGACTCTTTGATGATTATTAAAACATTGAATCATGAACATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTAAAGGCTATTTA  >  CP000730/2503352‑2503630
                                                                                                                                            |                                                                                                                                           
ccTCAAGACTCTTTGATGATTATTAAAACATTGAATCATGAACATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAc                                                                                                                                             >  1:327809/1‑140 (MQ=255)
  ttaaGACTCTTTGATGCTTATTAAAACATTGAATCATGAACATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACgg                                                                                                                                           <  1:10823/138‑1 (MQ=255)
  tCAAGGCTCTTTGATGATTATTAAAACATTGAATCATGAACATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAggg                                                                                                                                           >  2:308612/1‑140 (MQ=255)
         tctTTGATGATTATTAAAACATTGAATCATGAACATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGc                                                                                                                                    <  1:107445/140‑1 (MQ=255)
                                          cATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACt                                                                                                   >  2:354839/1‑140 (MQ=255)
                                          cATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACt                                                                                                   >  1:327807/1‑140 (MQ=255)
                                                     tcaaCATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATc                                                                                        <  1:148489/140‑1 (MQ=255)
                                                         cATTTGCTTCCAATACGCGATCAACTTGTTTATAAAGAAAAAAAAATTATAGATGAATTCGTC‑ATAATGAATGAGCCTAAAAGCGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTAAAATCacga                                                                                    >  1:266412/1‑140 (MQ=255)
                                                                    aaTAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCATACATTTGAACTCAAGTTAAATCACGACGAGGTTAAat                                                                         >  2:56407/1‑140 (MQ=255)
                                                                     aaaaacgctCAACTTGTTTATATAAAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAACTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGACGTTAAata                                                                        <  2:196324/132‑1 (MQ=255)
                                                                     aTAAGAGATCCACTTGTTTATCTAGAAAAAGAAGTAATCAGTCGA‑TAGACTATTTCGAATGACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAata                                                                        <  1:257543/140‑1 (MQ=255)
                                                                         gcgctgatcgtGTTTATATAGTAAAAGAAATTGTAGATGAATTAGAC‑AAAACGAATGACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGCCGAGGTTAAATAGAtt                                                                    <  2:41421/130‑1 (MQ=255)
                                                                         gCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATAtt                                                                    <  2:142248/140‑1 (MQ=255)
                                                                         gCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATAtt                                                                    <  1:233390/140‑1 (MQ=255)
                                                                           gATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTg                                                                  <  2:239855/140‑1 (MQ=255)
                                                                                aCTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTa                                                             <  1:123336/140‑1 (MQ=255)
                                                                                    gTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAAcc                                                         <  1:267100/140‑1 (MQ=255)
                                                                                      ttATATTGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACGGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCaa                                                       <  1:23471/140‑1 (MQ=255)
                                                                                              taaaaTAAATTTTAGATGAATTAGAT‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATa                                               <  2:24990/139‑1 (MQ=255)
                                                                                              aaaaaGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTGCTCCAAGCTaactc                                               <  2:165146/140‑4 (MQ=255)
                                                                                                 aaGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTa                                            <  2:198305/140‑1 (MQ=255)
                                                                                                     aaTTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAg                                        >  1:124964/1‑140 (MQ=255)
                                                                                                     aaTTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAg                                        <  1:151819/140‑1 (MQ=255)
                                                                                                      aTTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGt                                       >  1:252268/1‑140 (MQ=255)
                                                                                                        tttAGATGAATTATAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTcc                                     <  1:90860/140‑1 (MQ=255)
                                                                                                        tttAGATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTTAAGTGAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTcc                                     <  1:315760/140‑1 (MQ=255)
                                                                                                        tgtAGATCAATTCGAC‑AAAACGAGTTACCGTAAAACGGAGCTTGCAGTTCAAGCAGGAACTGGTCGTACATTTGAACTCAAGTTAAATCTCGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACGAACAGTcc                                     <  2:262762/138‑1 (MQ=255)
                                                                                                           aGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCctct                                  >  1:274042/1‑140 (MQ=255)
                                                                                                            gATGAATTAGAC‑AAAACGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCtt                                 <  1:308612/140‑1 (MQ=255)
                                                                                                               gAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTtgtg                              <  1:193217/140‑1 (MQ=255)
                                                                                                                caTTAGAC‑AAGACGAACTCACCTTGAACGTCGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTgtgt                             <  2:140840/139‑1 (MQ=255)
                                                                                                                    aGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGATAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAg                         >  1:23866/1‑140 (MQ=255)
                                                                                                                      aC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGTTTTAGACTAAGCTAAATACTAACAGTCCTCTTGTGATTAGtt                       >  1:47308/1‑140 (MQ=255)
                                                                                                                         aaaaCGAATTACCCTAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAACTACTAACAGTCCTCTTGTGTTTAGTTTc                     >  2:316051/1‑140 (MQ=255)
                                                                                                                          aaaCGAATTACCCGAAAAGGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCt                    <  2:49798/140‑1 (MQ=255)
                                                                                                                                 ttACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACTTTTGAACTCAAGTTATATCTCGTCGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTa             <  1:36076/140‑1 (MQ=255)
                                                                                                                                 ttACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTa             <  1:36104/140‑1 (MQ=255)
                                                                                                                                  tACCCTAAAAAGAGGCTTTCAGTTCAATAAGAAACTGGTCGTGCATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTaa            >  2:208419/1‑140 (MQ=255)
                                                                                                                                   aCCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTaaa           >  1:245645/1‑140 (MQ=255)
                                                                                                                                       tAAAACGGAGCTTGCAGTTCAAGAAGCAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTAAAGGCt       <  2:193809/140‑1 (MQ=255)
                                                                                                                                            cGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTAAAGGCTATTta  >  2:165463/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CCTCAAGACTCTTTGATGATTATTAAAACATTGAATCATGAACATGGTTCAATTCAACATTTGCTTCCAATAAGCGATCAACTTGTTTATATAGAAAAAGAAATTTTAGATGAATTAGAC‑AAAACGAATTACCCTAAAACGGAGCTTGCAGTTCAAGAAGAAACTGGTCGTACATTTGAACTCAAGTTAAATCACGACGAGGTTAAATATATTTGCTTTAAACCAAGCTAAATACTAACAGTCCTCTTGTGTTTAGTTTCTTACGTTAAAGGCTATTTA  >  CP000730/2503352‑2503630

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