Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,961 T→C 34.8% intergenic (+199/‑347) kefB‑III → / → PP_mr64 glutathione‑regulated potassium/H+ antiporter/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,9610TC34.8% 18.3 / 16.0 23intergenic (+199/‑347)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (6/9);  new base C (2/6);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

CGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCA  >  NC_002947/5988813‑5989097
                                                                                                                                                    |                                                                                                                                        
cgcccgcGCGGCGCTTCGCGGGTAAACCAGACACTCCGTTTGTTTCGGGCCAGTTATTCCTATGGGACGCGCACGCGAGCCCCCTGGCCGACGCCTCGATATCGCGTCGTATAAACAAGGCGGTCGCGCGCGCGTGGCACAGGAGGCAtt                                                                                                                                         <  2:483742/150‑1 (MQ=25)
     gcgcgGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCa                                                                                                                                    <  2:288482/150‑1 (MQ=255)
          gcgcTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaac                                                                                                                               <  2:189660/150‑1 (MQ=255)
             cTTCGCGGGTAATCCCGCTCCTTCGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGagacaaa                                                                                                                            <  2:453455/150‑1 (MQ=255)
              ttCGCGCGTGAACACGCTCCTACGTTTGGTTCGGGACAGTTATTCCTGTGCGACGCGCACGCGAGCCCTGTGGCGCACGCCGCTATAGCGCGTCGTACAAACAAGGCGGTCCCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                           <  2:126670/150‑1 (MQ=255)
                cgcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGt                                                                                                                         <  2:72221/150‑1 (MQ=255)
                 gcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAg                                                                                                                       <  1:308643/150‑1 (MQ=9)
                   gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGACCGAAACAAACGTAGGa                                                                                                                     >  1:294359/1‑150 (MQ=2)
                     gTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGTCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAg                                                                                                                    >  1:77552/1‑150 (MQ=255)
                       aaaCCCGCTCCTACGTTTGTTTCGGGCTGTTTATTCCTGCGGGGTTTGCGCGTGACCGCCTTGGGGCATGTCTCGCTATCGCGTCGTACCAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAATGTAGGAGCgg                                                                                                                 <  2:191503/150‑1 (MQ=255)
                       aaaCCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCCCGACTCGATATCGCGTCGTACAAACAAGGCGGTTGGGCGCGCCTTGCAGATGAGGGATTGGCCAGAAACAAACGTAGGAGCg                                                                                                                  <  2:424967/150‑1 (MQ=37)
                           ccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGtt                                                                                                              >  2:466497/1‑150 (MQ=255)
                                cctcagtttgtttcgggcCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAccc                                                                                                         >  2:132307‑M2/6‑147 (MQ=255)
                                          tttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgcc                                                                                               <  2:176947‑M2/150‑14 (MQ=255)
                                                                                 cccTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTa                                                        >  1:119984/1‑150 (MQ=255)
                                                                                            gtctcGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacc                                            <  1:393353‑M2/148‑65 (MQ=255)
                                                                                                            gtACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAAAGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgcc                            <  2:10374‑M2/150‑81 (MQ=255)
                                                                                                             tACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCa                            >  1:221284/1‑150 (MQ=255)
                                                                                                                aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg                                                                      >  1:48713/1‑105 (MQ=255)
                                                                                                                aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg                                                                      <  2:48713/105‑1 (MQ=255)
                                                                                                                aaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtcc                        <  1:420105‑M2/150‑85 (MQ=255)
                                                                                                                                       ggCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCa  <  1:176947/150‑1 (MQ=255)
                                                                                                                                              gaggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgcc                            >  1:438477‑M2/1‑36 (MQ=255)
                                                                                                                                              gaggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgcc                            <  2:438477‑M2/116‑81 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
CGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCA  >  NC_002947/5988813‑5989097

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