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,2980AT99.4% 57.3 / ‑6.9 24T15T (ACA→ACTUSA300HOU_2392FemAB family peptidoglycan biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (9/14);  minor base C (0/1);  total (9/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 = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCAC  >  CP000730/2532162‑2532434
                                                                                                                                        |                                                                                                                                        
gtgtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGCGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTActacaa                                                                                                                                      <  1:179123/141‑1 (MQ=255)
 tgtATGTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGAATGAGTTCAATGCGTTTACTACaaa                                                                                                                                     >  1:78172/1‑141 (MQ=255)
     tttGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAt                                                                                                                                 >  2:214342/1‑141 (MQ=255)
              aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACAt                                                                                                                        >  1:92861/1‑141 (MQ=255)
                ttttttATAAAATAAAAAAATAGCAGTATTGATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATATTTCACATTa                                                                                                                      <  2:59975/141‑1 (MQ=255)
                                aaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACTTTACACACAATCAGCTAtt                                                                                                      >  1:249100/1‑141 (MQ=255)
                                 aaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCGGCTATTc                                                                                                     <  1:135756/141‑1 (MQ=255)
                                    tAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCCAGGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAtt                                                                                                  <  1:35556/141‑1 (MQ=255)
                                          tattaATAAATTAAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATc                                                                                            <  1:73304/141‑1 (MQ=255)
                                                     tCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAt                                                                                 <  2:10586/141‑1 (MQ=255)
                                                            ctataGTTCAATTCAGAGAGGGATGGGATAAAATGAATTTTCTAACGTTGACTTCGGATGAGTTCAATGCGCTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATGCATTACAATCATAGAGTTGATTTAAAAg                                                                          <  2:210846/140‑1 (MQ=255)
                                                                 gTTCAATTAAGAGACGGATGGGATAAAATGAACTTTGTAACGTGGACTCCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAt                                                                     >  2:110996/1‑141 (MQ=255)
                                                                  ttCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtg                                                                    >  1:195617/1‑141 (MQ=255)
                                                                    aaaTTCGAAGAGGGATGGAATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgtg                                                                  <  2:189413/140‑1 (MQ=255)
                                                                         cAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                             <  2:92861/141‑1 (MQ=255)
                                                                                         aaaaTGAACTTCGTCCCGTTGACTTCGGCCGAGCTCACTCAGTTTAACAAAACGCTGTTTTCACCTTACCAATGAGCTGCTATTAGTTAGATGTTTAGAGTTGATTTAAAAGGCGACGTGCATCTTGTAGGTGTTAAAGAt                                             <  1:182563/141‑1 (MQ=255)
                                                                                            aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATTCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                          <  2:216977/141‑1 (MQ=255)
                                                                                               aaCTTTGTAAAGTTGACTTCGGATGAGTTCAGTGCGTTTATTAGAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGTTGACAAt                                       >  2:146096/1‑141 (MQ=255)
                                                                                                aCTTTGTAACGTTGACTTCGGATGAGTTCACTGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGAATATTGTAGGGGTTAAAGATGACAATg                                      >  1:59283/1‑141 (MQ=255)
                                                                                                aCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTGGGGGTTAAAGATGGCAATg                                      >  1:83237/1‑141 (MQ=255)
                                                                                                    tGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCTTTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCa                                  <  2:199220/141‑1 (MQ=255)
                                                                                                      tcccGTTGACTCCGGATTAATCCAATGCGTTTTCTAATAAGTATTTTTCACATAAGACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAg                                <  1:130647/138‑1 (MQ=255)
                                                                                                        aCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTg                              <  1:174026/141‑1 (MQ=255)
                                                                                                                                gCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGaa      <  1:249019/141‑1 (MQ=255)
                                                                                                                                    ttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACATAAgcac  <  2:138117/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCAC  >  CP000730/2532162‑2532434

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: