Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,988,905 T→C 40.0% 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
*AE0154515,988,9050TC40.0% 35.2 / 47.4 40intergenic (+199/‑347)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (16/8);  new base C (9/7);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

ATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCATTT  >  AE015451/5988748‑5989113
                                                                                                                                                             |                                                                                                                                                                                                                
aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAtt                                                                                                                                                                                                                 <  1:4998/159‑1 (MQ=255)
                        tCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCa                                                                                                                                                                                                            <  2:699942/140‑1 (MQ=255)
                                                ttgtttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCCGCGAAg                                                                                                                                                                           >  2:617301/1‑148 (MQ=11)
                                                       gggCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgggcg                                                                                                                                                              >  1:158149/1‑154 (MQ=21)
                                                                          gatttgcgcgcgaccgCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGGCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgcaaatcctccggcaaacacctcacaatt                                                                                                            <  2:271381‑M2/180‑73 (MQ=255)
                                                                          gACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcg                                                                                                                                                                  <  1:724961‑M2/132‑19 (MQ=255)
                                                                          gACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCCTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcg                                                                                                                                                                  >  2:724960‑M2/1‑114 (MQ=255)
                                                                           aCGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGAt                                                                                       >  1:267015/1‑206 (MQ=255)
                                                                              cgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCGCGATCTCATGGGCGCCCCAAAACTATCGc                                                                                                                            >  2:238295/1‑166 (MQ=255)
                                                                                 gcgcgaCCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgct                                                                                                                                         <  2:197112‑M2/148‑44 (MQ=255)
                                                                                 aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGTGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTg                                                                                                      >  1:236007/1‑185 (MQ=255)
                                                                                 aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGGTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCCTCCTGCAATTTCCTGCCGATCACCTCGCCGCCATGAcgcgc                                                                  >  1:186535/1‑220 (MQ=255)
                                                                                 aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTAt                                                                                                                               >  1:493008/1‑160 (MQ=255)
                                                                                 aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTca                                                             >  1:673970/1‑226 (MQ=255)
                                                                                            caggCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGt                                         >  1:188518/3‑235 (MQ=255)
                                                                                                cgcACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAAt                                            >  2:740108/1‑228 (MQ=255)
                                                                                                       cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgagctcacaggcgctgaaaatcctccggcaaacaccacaaaattcca                                                                                                         <  2:33333‑M2/159‑76 (MQ=255)
                                                                                                        tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCt                            >  1:462049/1‑236 (MQ=255)
                                                                                                          gATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCt                            <  2:493007/234‑1 (MQ=255)
                                                                                                           atatCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCt                            <  2:462048/233‑1 (MQ=255)
                                                                                                              tCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccg                                                                                                                             >  2:158148‑M2/1‑77 (MQ=255)
                                                                                                              tCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcaca                                                                                                                                               <  1:612456‑M2/114‑38 (MQ=255)
                                                                                                              tCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcaca                                                                                                                                               >  2:612455‑M2/1‑77 (MQ=255)
                                                                                                                    cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcg                                                                                                                                                           <  2:313192‑M2/97‑26 (MQ=255)
                                                                                                                    cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATtct                                 >  2:356973/1‑219 (MQ=255)
                                                                                                                    cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTc                >  1:604101/1‑236 (MQ=255)
                                                                                                                      tACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccag             >  1:617302‑M2/1‑69 (MQ=255)
                                                                                                                        gaaacaagGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccaga            >  1:373793‑M2/2‑67 (MQ=255)
                                                                                                                         aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTcccc             >  2:736453/1‑234 (MQ=255)
                                                                                                                                   gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGAcgcg                                                                   >  1:415566/1‑170 (MQ=255)
                                                                                                                                   gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAAcc                                                     >  1:564967/1‑184 (MQ=255)
                                                                                                                                   gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCAttt  <  1:740109/235‑1 (MQ=255)
                                                                                                                                   gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCAttt  <  1:736454/235‑1 (MQ=255)
                                                                                                                                   gTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgc                                                                                                                                                          >  1:438561‑M2/1‑56 (MQ=255)
                                                                                                                                   gTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgc                                                                                                                                                          <  2:438560‑M2/82‑27 (MQ=255)
                                                                                                                                      gcgcgcgcCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATcaca                                               >  1:153365/1‑187 (MQ=255)
                                                                                                                                      gcgcgcgaCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAccc                                                                                                                                                                                 <  2:580920/56‑1 (MQ=9)
                                                                                                                                      gcgcgcgaCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAccc                                                                                                                                                                                 >  1:580921/1‑56 (MQ=17)
                                                                                                                                       cgcgcgaCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgccca                                                   <  2:373792‑M2/181‑130 (MQ=255)
                                                                                                                                         cgcgaCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctc                                              >  1:197113‑M2/1‑50 (MQ=255)
                                                                                                                                             ccTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGa                                                                  <  2:564966/161‑1 (MQ=255)
                                                                                                                                                             |                                                                                                                                                                                                                
ATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCATTT  >  AE015451/5988748‑5989113

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