Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,368,556 T→C intergenic (+257/‑104) PP_2962 → / → PP_2964 alcohol dehydrogenase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,368,5560TC94.7% 59.4 / ‑3.9 19intergenic (+257/‑104)PP_2962/PP_2964alcohol dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (7/11);  total (7/12)
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.11e-01

GGTCAAGCTGGCTGCAGTCGAAGTGGTGACTTCGTCGCTACTGTCGGGAGGAAAGGGTGAGGCAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGTTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTGGGGGATGTGTTCAATGGCATCGGT  >  NC_002947/3368423‑3368679
                                                                                                                                     |                                                                                                                           
ggTCAAGCTGGCTGCAGTCGAAGTGGTGACTTCGTCGCTACTGTCGGGAGGAAAGGGTGAGGCAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTc                                                                                                                           <  1:669550/136‑1 (MQ=255)
                       tggtgACTTCGTCGCTACTGTCGGGAGGAAAGGGTGAGGCAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGt                                                                                                    <  2:854134/136‑1 (MQ=255)
                       tggtgACTTCGTCGCTACTGTCGGGAGGAAAGGGTGAGGCAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGt                                                                                                    >  1:1081661/1‑136 (MQ=255)
                                cgtcgCAACTGTCGGGAGGAAAGGGTGCGGCAGAAAAATCTATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGCTCATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCgagg                                                                                           >  1:679082/1‑135 (MQ=255)
                                                        gTGAGGCAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCAc                                                                   >  1:1096834/1‑136 (MQ=255)
                                                              cAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTg                                                             >  2:528973/1‑136 (MQ=255)
                                                              cAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTg                                                             <  2:267768/136‑1 (MQ=255)
                                                               aGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTgg                                                            <  2:724151/136‑1 (MQ=255)
                                                               aGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTgg                                                            <  2:212150/136‑1 (MQ=255)
                                                                 aaaaaTCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCt                                                          <  2:287301/136‑1 (MQ=255)
                                                                    aaTCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGc                                                       <  2:229836/136‑1 (MQ=255)
                                                                                               ggggTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCt                            <  1:132938/136‑1 (MQ=255)
                                                                                               ggggTTATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCt                            <  2:330243/136‑1 (MQ=255)
                                                                                                   ttATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTggg                        >  2:1152601/1‑136 (MQ=255)
                                                                                                   ttATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTggg                        >  1:879227/1‑136 (MQ=255)
                                                                                                    tATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTgggg                       <  2:1096834/136‑1 (MQ=255)
                                                                                                    tATGCCGAGATAAGGCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTgggg                       >  2:934415/1‑136 (MQ=255)
                                                                                                                 ggCAAAAATTAGGACATTCGCTCTGTAAATATATGATTTAAAATGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTGGGGGATGTGTTCAATg          <  2:1109420/136‑1 (MQ=255)
                                                                                                                         ttAGGACATTCGTTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTGGGGGATGTGTTCAATGGCATCGGt  <  1:195351/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                           
GGTCAAGCTGGCTGCAGTCGAAGTGGTGACTTCGTCGCTACTGTCGGGAGGAAAGGGTGAGGCAGAAAAATCGATCATTCACGCGGGGTTATACCGGGGTTATGCCGAGATAAGGCAAAAATTAGGACATTCGTTCTGTAAATATATGATTTAAAAGGTTATTCGAGAGGCCGTGGCTTGCTGGTCATCCACCGTCTGGCTAGCCGTCTATGCTGATGCTTTTGCCTCGCTTGGGGGATGTGTTCAATGGCATCGGT  >  NC_002947/3368423‑3368679

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