Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,476,692 G→T 24.3% H19N (CAT→AAT)  U712_17855 ← putative undecaprenyl‑phosphate galactose phosphotransferase‑like Protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,476,6920GT24.3% 87.0 / 6.8 56H19N (CAT→AAT) U712_17855putative undecaprenyl‑phosphate galactose phosphotransferase‑like Protein
Reads supporting (aligned to +/- strand):  ref base G (30/11);  new base T (3/12);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 5.57e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCTCCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTAAA  >  CP006881/3476555‑3476829
                                                                                                                                         |                                                                                                                                         
cGTCTCCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGt                                                                                                                                          <  4:180939/139‑1 (MQ=255)
    tcCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATTTAAAAGAATGTAttt                                                                                                                                      >  6:388225/1‑139 (MQ=255)
             aTCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGGATTCAGATAAAGGAAAAGAATGTATTTTTTGAAcgg                                                                                                                             >  3:113075/1‑137 (MQ=255)
                aaaaTTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGGTTTCTGATAATGTAAAAAAAAGTATTTGTTGATCagag                                                                                                                           >  2:334664/1‑138 (MQ=255)
                       aaaaTACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCATATAAAGTAAAATAAATTATTTTTTGAACTGAGaaataacc                                                                                                                   >  4:80219/1‑133 (MQ=255)
                                 aaTCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTAAGCTAATTTAAAAAAATGTATTTTTTGTACTGATGATTCACGAAGCACAtt                                                                                                          >  1:405397/1‑138 (MQ=255)
                                      ttGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCt                                                                                                    <  7:372579/139‑1 (MQ=255)
                                      ttGGCAGCGCAATCGCAAATCCAATGAGCGCAAACAAAATATCCATTACGCGCTTTATTCCAAAATAACGCGCCGTCTTTTCATATAATTAAAAAGAATGTATTTGTTCAACTGACAATCCACGACTCACATTCATGCt                                                                                                    <  4:135001/139‑1 (MQ=255)
                                      agaccggcgaAATCGCAAATCAATTGACCGCAAACAAAAAATCAATTACGCGCTTTATTACAAAATACGGCGCCTTATTTTCATAAAATTTAAAAGATTTTTTTTGTGTAACTTATAATGCACGACTCACATTCATGCt                                                                                                    <  2:44864/133‑1 (MQ=255)
                                           aGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATTTAAAAGAATGTATTTGTTGTACTGAGAATTAACGAATCACATTCATGCtttta                                                                                                >  4:359220/1‑137 (MQ=255)
                                             cgcAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTTTTGAACTGAGAATTAACGACTCACATTCATGCTTTTctct                                                                                             >  8:71287/1‑139 (MQ=255)
                                             cgcAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGAGTATGTAAAAGAATGTTTTTGTTGAACTGAGAATTGACGACTCACATTCATGCTTTTatct                                                                                             >  5:412535/1‑139 (MQ=255)
                                               cAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGGTTTGAGATATTGTAAAAGAATTTATTTTTTGAACTGAGAATTCAAGCAGCACATTAATTCttttatatgg                                                                                           >  5:50135/1‑133 (MQ=255)
                                               cAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTATCTGc                                                                                           >  2:64680/1‑139 (MQ=255)
                                                actCGCAAATCCTATGCGAGAAAACAAAATTTCAATTACGAGTTTTTTAGAAAAATAACAAGCGTTCTTTTCTTATAATTTAAATTAATTTTTTTTTTTAACTTTTAATCCACGACTCACATTCATGCTTTTCTCTGca                                                                                          <  8:172592/137‑1 (MQ=255)
                                                aaccGCCAACCAAATGAGCGAAAACCAAATATCCATACCGCGCTTAATTGCCAAATAACCAGCCTTCTTTACAGAAAATTTAAAATATTTTATTTGTTCAACTTACAATCCACGACTCACATTCATGCTTTTCTCTGca                                                                                          <  3:31736/136‑1 (MQ=255)
                                                  tCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTGACGACTCACATTCATGCTTTTCTCTGcccg                                                                                        >  1:8859/1‑136 (MQ=255)
                                                    gCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCATAAAATGTAAAACAATGTATTTGTTGAACTGAGAATCCACGACTCACATTCATGCTTTTCTCTGcactca                                                                                      <  3:195095/139‑1 (MQ=255)
                                                       aatccaatGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTTTTGAACTTAGAATTCACGACTCACATACATGCTTTTATCTGcacctgtc                                                                                    >  1:44864/1‑133 (MQ=255)
                                                       aatccaatGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCGTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGATCATTCACGACTCACATTCATGCTTTTATCTGcactcact                                                                                    >  3:51279/1‑138 (MQ=255)
                                                            attgcgcgaAAACAAAAAAGCTATAACGCGATTTATTCCCAAAAAACGAGACCTCTTTTAATAAAATTTAAAATAATTAATTTTTTGTACTCCGAATGCACGACTCACATTCATGCTTTTCTCTGCACTCACTGacaca                                                                              <  4:429041/134‑1 (MQ=255)
                                                              tGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGGTTTTAGATAATTTAAAAGAAATTATTTTTTGAAATGAGTATTAAAGCCTCACATTCAAGGTTTTGTGTGCGCTCACGGacacacc                                                                            >  1:225833/1‑138 (MQ=255)
                                                              tGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTTAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGacacaca                                                                            >  7:15729/1‑139 (MQ=255)
                                                                acccaaaaCAAAAAATCAATTCCGCGTTTAATTCCAAAATACCGCGCCTTATTTTCATATAATTTAAAATATTTTATTTGTTTAAATTACAATCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGacacacac                                                                           <  5:133840/133‑1 (MQ=255)
                                                                  cgcAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTGTGTTAACATAATGTAAAATAATGTATTTTTTTAACTTAGAAATAAAGACACACAATAATGCTTTTTTCTGCACTCACTGCCCCCCAcccc                                                                        >  6:268605/1‑137 (MQ=255)
                                                                      aaCCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATTTAAAAGAATGTATTTGTTGAACTGAGTATTCACGACTCACATTCATGCGTTTCTCTGCGCTCGCTGACACCCACCTCAc                                                                      >  3:53219/1‑137 (MQ=255)
                                                                                 cATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATTTAAAATAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTggg                                                          <  7:69975/138‑1 (MQ=255)
                                                                                  aTTACGCTATTTATTGCAAAATAAGGCGCATTATTTTCATATAATTAAAAATAATTTTTTTTTTAAACTTATAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGAt                                                        <  3:80219/139‑1 (MQ=255)
                                                                                   ttccgcgCTTTATAGCCAAATCACGAGCCGTCTTTTCATATAATTAAAACTAATTTATTTGTTCAACTTACAATCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATg                                                       <  4:265771/136‑1 (MQ=255)
                                                                                   ttACGCGCTTTATTGCCAAATAACGAGCCGTGTTTTCAGATAATGTTAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATTCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATg                                                       >  3:168250/1‑139 (MQ=255)
                                                                                   ttACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCCCATTCATGCTTTTCTCTGCACTCACTGACACACCCCTCACGATTTGTTTGGGATg                                                       >  7:247441/1‑139 (MQ=255)
                                                                                        cgcTTTATTCCAAAAAAACGCGCCTTTTTTTCATATAATGAAAAAGAATGTTTTTGTTGAACTTAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTatg                                                  <  6:388086/139‑1 (MQ=255)
                                                                                        cgcTTTATATCCAAATAACGAGACTTTTTTTCAGATAATTTAAAACTATTTATTTGTTGAACGGATAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTatg                                                  <  5:111092/139‑1 (MQ=255)
                                                                                         tatttATTGCAAAAAAACGACCGTTTTTTTAAGATATTTAAAAAGATTTTTTTTGTTCAACTGATAATCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTatga                                                 <  6:50135/137‑1 (MQ=255)
                                                                                         gatttATTGAAAATTAACGAGCTGTATTTACTTATATTTAAAAATAATGTATTTGTTGAAATTACAATACACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTatga                                                 <  5:362829/137‑1 (MQ=255)
                                                                                              atcccaaaaTAACCCGCCTTATTTACAGAAAATTAAAAAGATTTTTTTTTTTCAACTGAGAAACCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGAc                                             <  1:334664/132‑1 (MQ=255)
                                                                                              aTTGCCAAATACCGAGCCGTCTTTTCATATAATTTAAAAGAATTTATTTGTTCAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGAc                                             <  5:217909/138‑1 (MQ=255)
                                                                                                tGCCAAATAACGAGCCGTCGTTTGAGATTATTTAAAATAATGTATTTGTTTAACTGATAATTAACGACTCACATTGATGAGTTTCTCGGCACAAACTGACACACACCGACCGCTTTGTTTGGGGTGTTATTATGACATa                                          >  6:5126/1‑139 (MQ=255)
                                                                                                  ccAAATAACGAGCCGTGTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCCCGACTCACATTGCTGCTTTTCTCTGCACTCACTGACACCCACCGCCCGATTTGGTTGGGATGGTTTGATGACCTATg                                        >  7:124305/1‑139 (MQ=255)
                                                                                                     aaacaCGCGCCGTCTTTTCATAAACTGAAAAAGATTGTATTTGTTGAACTGAGACATCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCca                                     <  8:247441/135‑1 (MQ=255)
                                                                                                               cGTCTTTCCATATAATTAAAAAGAATGTATTTGTTCAACTTCCAACCCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAg                           <  4:53219/139‑1 (MQ=255)
                                                                                                                 tcgtttAAGATAATGTAAAAGAATGTATTTGTTTTAATTAGAATTAACGCCTCACAATAATGAGTTTCTCTGCACTCCCGGACACACACCAACCGCATTGTTTTGGATGTTATTATTACCTATAGAACTCAGTTTATCa                         >  4:376042/4‑139 (MQ=255)
                                                                                                                 tCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCGCACGATTTGTTTGGGTTGTTATGATGACCTATGCCAGTCAGTTTAg                           >  2:37361/1‑137 (MQ=255)
                                                                                                                             ttGTAAAAGATTGTTTTTTTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTaa             <  4:51279/138‑1 (MQ=255)
                                                                                                                             aTGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTaa             >  7:267596/1‑139 (MQ=255)
                                                                                                                             aTGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGAAGACGGATGTTaa             >  3:118274/1‑139 (MQ=255)
                                                                                                                                 aaaaTAATTTATTTGTTCAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTg          <  7:71287/138‑1 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCcagt                                   <  6:197823/106‑1 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCcagt                                   >  5:197823/1‑106 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTaaa  >  5:383179/1‑139 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTaaa  >  3:73389/1‑139 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTaaa  >  7:379561/1‑139 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTaaa  >  2:183866/1‑139 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTaaa  >  1:308948/1‑139 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTaaa  >  1:308073/1‑139 (MQ=255)
                                                                                                                                        tGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTGTAGCAGACGCATGTAAATTCGATGTaaa  >  1:437268/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CGTCTCCAAACAAATCAAAATTGAAAATACCGCAATCATTGGCAGCGCAATCGCCAATCCAATGAGCGCAAACCAAATATCCATTACGCGCTTTATTGCCAAATAACGAGCCGTCTTTTCAGATAATGTAAAAGAATGTATTTGTTGAACTGAGAATTCACGACTCACATTCATGCTTTTCTCTGCACTCACTGACACACACCTCACGATTTGTTTGGGATGTTATGATGACCTATGCCAGTCAGTTTAGCAGACGGATGTTAATTGGATGTAAA  >  CP006881/3476555‑3476829

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