Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,988,911 Δ1 bp 40.0% intergenic (+205/‑341) 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,9110A.40.0% 92.9 / 58.3 40intergenic (+205/‑341)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (16/8);  new base . (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 = 4.35e-01

TCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCATTTCCTACGTTCACACCAATCCCAAGGA  >  AE015451/5988772‑5989138
                                                                                                                                           |                                                                                                                                                                                                                                   
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)
                                                                                                                                      tGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCATTTCCTACGTTCACACCAATCCCAAGGa  <  2:267014/233‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                                                                                                   
TCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGAACCATTTCCTACGTTCACACCAATCCCAAGGA  >  AE015451/5988772‑5989138

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: ‑