Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,532,298 A→T 100% T15T (ACA→ACT USA300HOU_2392 → FemAB family peptidoglycan biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,532,2980AT96.1% 64.1 / ‑5.3 23T15T (ACA→ACTUSA300HOU_2392FemAB family peptidoglycan biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (4/18);  minor base C (0/1);  total (4/19)
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.39e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532162‑2532433
                                                                                                                                        |                                                                                                                                       
gtgtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTActacaa                                                                                                                                     >  1:88738/1‑141 (MQ=255)
          aaTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTc                                                                                                                           <  1:323798/141‑1 (MQ=255)
              agattttGATAACATAAAAAAATCGCATTATTACTAAATTCAAAAATAATATATACAAACCGAGAGGGATGGGATAAAAAGAACTAAGTAACGTTGACTTCGGATGGGTTCAATGCGTTTACTACAAAGCATTTTTCAAAt                                                                                                                       >  2:87276/4‑141 (MQ=255)
                         aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATc                                                                                                            <  2:172140/141‑1 (MQ=255)
                            taaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                         <  2:282099/141‑1 (MQ=255)
                            taaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                         <  1:67189/141‑1 (MQ=255)
                                 aaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAAGTCAGATAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACATACAATCTGCTATTc                                                                                                    <  2:289679/141‑1 (MQ=255)
                                  aaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGATTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCa                                                                                                   <  1:324692/141‑1 (MQ=255)
                                    tAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAtt                                                                                                 <  2:319572/141‑1 (MQ=255)
                                                 aaaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGt                                                                                    <  1:176870/141‑1 (MQ=255)
                                                  aactCAAAAAGTATAGTTTAATTTAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGTATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                   <  2:223229/138‑1 (MQ=255)
                                                  aaTTCAATAATTATAGTCAAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                   <  2:58397/141‑1 (MQ=255)
                                                                  ttCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtg                                                                   <  2:175659/141‑1 (MQ=255)
                                                                           gagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTg                                                          <  2:101875/141‑1 (MQ=255)
                                                                               gtgatgggatAAAAATAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAgg                                                      <  1:231270/139‑1 (MQ=255)
                                                                                   agggctACAATGAACTTTGTAACGTGGACTTCGGAATAGTTCAATGCGTTTACTACAAAGCAATTGTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAACGCGATGTGCATCTTGTAGGGGtt                                                  <  2:88709/140‑1 (MQ=255)
                                                                                     ggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGGTGTGCATCTTGTAGGGGttga                                                >  1:265552/1‑139 (MQ=255)
                                                                                               aaCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGCTTTAAAAGGCGCTGTGCATCTTGTAGGGGTTAAATATGACAAt                                      >  2:102048/1‑141 (MQ=255)
                                                                                                 cTTTGTAACGGTGACGTCGGATGAGTTCAATGTGTTTACTACAAGGCTTTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATgg                                    <  1:87276/141‑1 (MQ=255)
                                                                                                              aCTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCa                       >  1:276139/1‑141 (MQ=255)
                                                                                                                  cGGATGAGTTCAATGCGTCTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt                   <  1:90748/141‑1 (MQ=255)
                                                                                                                                  gTTTACCACAAAGCATTTTTCACATTACACACAATCAGATGTTCATTAAGAACTGAGACTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATTGTCAAGTGATTGCAGGATGCTTATTGACAGAAgc   <  1:9789/141‑1 (MQ=255)
                                                                                                                                   tttACTACAAAGCTTTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGGTCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  <  2:102008/141‑1 (MQ=255)
                                                                                                                                   tttACTACAAAGCATTTTACACATTACACACAATCAACTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  <  1:322145/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532162‑2532433

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: