Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 230,497 A→G 11.0% H323R (CAC→CGC)  U712_01065 → Fatty‑acid peroxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881230,4970AG11.0% 88.3 / 5.1 38H323R (CAC→CGC) U712_01065Fatty‑acid peroxygenase
Reads supporting (aligned to +/- strand):  ref base A (22/10);  new base G (2/2);  total (25/13)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

TGTTTGTGCAGGAGGTCCGCAGATATTATCCGTTCGGCCCGTTTTTAGGGGCGCTTGTCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAGGGGAAGGCATTA  >  CP006881/230362‑230634
                                                                                                                                       |                                                                                                                                         
tGTTTGTGCAGGAGGTCCGCAGATATTATCCGTTCGGCCCGTTTTTAGGGGCGCTTGTCAAAAAAGATTTTGTATGGAATAACTGTTAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACgc                                                                                                                                        >  8:78224/1‑138 (MQ=255)
    tgtgCAGGAGGTCCGCAGATATTATCCGTTCGGCCCGTTTTTAGGGGCGCTTGTCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCt                                                                                                                                    >  5:72832/1‑139 (MQ=255)
                         ttatCCGTTCGGCCCGTTTTTAGGGGCGCTTGGCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGACCCACGACCCTCGTCTATGGGATCATCCCGAt                                                                                                               >  2:88160/1‑139 (MQ=255)
                         ttatCCGTTCGGCCCGTTTTTAGGGGCGATTGTCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGACGGGCACAGCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGAt                                                                                                               <  7:228574/139‑1 (MQ=255)
                              cGTTCGGCCCGTTTTTCGGGGCGCTTGTCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGACCGACCCACGACCCTCGTCTATGGGATCATCCCGATGAAtt                                                                                                          <  3:86373/139‑1 (MQ=255)
                                                   cgcTTGTCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCg                                                                                       >  5:177292/1‑137 (MQ=255)
                                                                      tGTTTGCAAAAACTTTGATTTTAAGACGGGCACATGGGTGCTGCTTTATTTATATGCAACCACCCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGtt                                                                   <  1:237246/138‑1 (MQ=255)
                                                                            gAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGatat                                                             >  6:295788/1‑138 (MQ=255)
                                                                             aaTAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGCGCGGGAAGAAAATCTGTTTGATATGa                                                           >  7:297015/1‑139 (MQ=255)
                                                                             aaTAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGa                                                           <  6:72832/139‑1 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGTTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAgg                                                 >  7:181386/1‑139 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCGCGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCg                                                                                  >  3:314433/1‑106 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCGCGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCg                                                                                  >  3:319398/1‑106 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCGCGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCg                                                                                  <  4:314433/106‑1 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCGCGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCg                                                                                  <  4:319398/106‑1 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAgg                                                 >  2:212153/1‑139 (MQ=255)
                                                                                       aGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAgg                                                 >  6:146286/1‑139 (MQ=255)
                                                                                             agaagGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCggggg                                           >  3:337322/1‑139 (MQ=255)
                                                                                                 gggCACATCGGTGCTGCTTGATTTATATGGAACGAACCCCTACCCTCGTCTATGGGATCATCACGAGTAATTCCGGCCGGAACGATTTGCGGAGAGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCAc                                       >  6:137829/1‑139 (MQ=255)
                                                                                                  ggCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACg                                      >  2:318441/1‑139 (MQ=255)
                                                                                                  ggCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACg                                      >  4:261583/1‑139 (MQ=255)
                                                                                                    cacaTCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACTcc                                    >  1:199574/1‑139 (MQ=255)
                                                                                                        ccGGTGCGGCTTGATTATTTTGCACCCCCCCCCGCCCCTCGTCTATGGGAGCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCgaga                                <  4:307957/138‑1 (MQ=255)
                                                                                                              ctgctTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCg                                                                                  >  8:87427/1‑83 (MQ=255)
                                                                                                              ctgctTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCg                                                                                  <  7:87427/83‑1 (MQ=255)
                                                                                                                 ctTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCAcc                       >  6:313139/1‑139 (MQ=255)
                                                                                                                    gATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCg                      <  5:146286/137‑1 (MQ=255)
                                                                                                                      tttATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGt                  >  1:330292/1‑139 (MQ=255)
                                                                                                                        tataTGGAACGACCCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTc                 <  1:212153/138‑1 (MQ=255)
                                                                                                                           atGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCACGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAgg             >  4:258428/1‑139 (MQ=255)
                                                                                                                           atGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAgg             >  6:114423/1‑139 (MQ=255)
                                                                                                                           atGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAgg             >  7:94923/1‑139 (MQ=255)
                                                                                                                               aacgaacCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAgggg           >  4:139399/1‑137 (MQ=255)
                                                                                                                               aacgaacCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAgggg           >  7:60290/1‑137 (MQ=255)
                                                                                                                                acgaacCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAGGGGa          <  8:60290/137‑1 (MQ=255)
                                                                                                                                 cgaacCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAggg            >  2:60433/1‑134 (MQ=255)
                                                                                                                                 cgaacCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAggg            <  1:60433/134‑1 (MQ=255)
                                                                                                                                      cACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAGGGGAAGGCATTa  <  8:297015/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
TGTTTGTGCAGGAGGTCCGCAGATATTATCCGTTCGGCCCGTTTTTAGGGGCGCTTGTCAAAAAAGATTTTGTATGGAATAACTGTGAGTTTAAGAAGGGCACATCGGTGCTGCTTGATTTATATGGAACGAACCACGACCCTCGTCTATGGGATCATCCCGATGAATTCCGGCCGGAACGATTTGCGGAGCGGGAAGAAAATCTGTTTGATATGATTCCTCAAGGCGGGGGGCACGCCGAGAAAGGCCACCGCTGTCCAGGGGAAGGCATTA  >  CP006881/230362‑230634

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: