Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,021,817 C→T 100% V139I (GTT→ATT)  USA300HOU_1873 ← acyl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,021,8170CT91.3% 53.9 / ‑0.8 23V139I (GTT→ATT) USA300HOU_1873acyl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (12/9);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAACAATCTGTCTTTATATTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGACCAC  >  CP000730/2021680‑2021952
                                                                                                                                         |                                                                                                                                       
tAACAATCTGGCTTTATATTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAatat                                                                                                                                       <  1:130570/140‑1 (MQ=255)
   cAATCTGTCTTTATATTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAAtatcat                                                                                                                                    <  2:171142/140‑1 (MQ=255)
          tCTTTATATTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTgg                                                                                                                             >  1:75121/1‑140 (MQ=255)
                tatTTTTTGTCTGAGATATGTTGCACTTGAGCTAAATATGCAACTTGTTGTTTACTTTCAATGCGACGTGCTTTTCTTCTTTCAACATGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATac                                                                                                                       <  2:73523/140‑1 (MQ=255)
                tatTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATac                                                                                                                       <  1:269189/140‑1 (MQ=255)
                      tAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAgg                                                                                                                 <  1:31433/140‑1 (MQ=255)
                                 atGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGATACATGTTTAg                                                                                                      >  1:277047/1‑140 (MQ=255)
                                            gAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTtca                                                                                           >  1:272230/1‑140 (MQ=255)
                                            gAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTtca                                                                                           >  2:90945/1‑140 (MQ=255)
                                             aGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATAATCTTCTGGATATACACAAGGTACTTGTTTAGGTTTGCAtttct                                                                                          >  2:28641/1‑137 (MQ=255)
                                                     atTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTCTAGGTTTGCTTTCATCATCTAAc                                                                                  <  2:170151/140‑1 (MQ=255)
                                                          aTTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGcaac                                                                             >  1:176748/1‑140 (MQ=255)
                                                                tttttACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGt                                                                       <  2:32386/138‑1 (MQ=255)
                                                                   ttACTTTCAATGCGACGTGCTTTTCTTCGTTCACCACGTTGTGGTGCTGTATAATAGAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTGTCGGTTTGCCTTCATAATCTAACGCAACAAATGTTaa                                                                    >  2:278132/1‑140 (MQ=255)
                                                                       tttCAATGCGGCGTCCGTTGTTGCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTaa                                                                <  1:254412/140‑1 (MQ=255)
                                                                        ttCAATGCGACGTGCTTTTCTTCCTTCAACTCGTTGTCTCTCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCACGTACATTTTTAGGTTTGCCTTCATCATCTAACGCATCAAATGTTAAGTAAc                                                               <  2:109939/140‑1 (MQ=255)
                                                                                gACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATTCACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAAGGCa                                                       >  2:216575/1‑140 (MQ=255)
                                                                                 aCGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCac                                                      >  2:205593/1‑139 (MQ=255)
                                                                                     gCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATACACCAGGTAAATGTTTAGGTTTGCCTTCATCGTCTAACGCAACAAATGTTAAGTAACATAATGCAACTaa                                                  >  2:25026/1‑140 (MQ=255)
                                                                                     gCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTAACTCTGATATTCTGGATACACCGCATGTACAAGTGAAGCAGAGCCGTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTaa                                                  >  1:129108/1‑140 (MQ=255)
                                                                                              cGTTCATCATGTTGTGGTGCTGTATCATTAAACCTTTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAACGTTAAGTAACTTAATGCAGCTAAATCATGCtt                                         <  1:65332/140‑1 (MQ=255)
                                                                                                     cacGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGttatta                                  >  2:286738/1‑140 (MQ=255)
                                                                                                            tggtgCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATATTAACTAACTTAGTGCAGCTAAATCCTGCTTGTTATTAAATACAt                           >  1:88164/1‑140 (MQ=255)
                                                                                                                                   tttCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGAcc    <  1:267494/140‑1 (MQ=255)
                                                                                                                                     tCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTCCCTCAACCAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGACCAc  >  1:241933/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TAACAATCTGTCTTTATATTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGACCAC  >  CP000730/2021680‑2021952

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