Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,988,848 C→T 30.5% intergenic (+142/‑404) 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,8480CT30.5% 57.5 / 26.0 36intergenic (+142/‑404)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (9/16);  new base T (3/8);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATT  >  AE015451/5988628‑5989080
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        
tgaCCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATcgcg                                                                                                                                                                                                                            <  1:891672/235‑1 (MQ=255)
                        gtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                                                                                                                                   <  1:562563‑M2/123‑1 (MQ=255)
                         gCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                                                                                                                                   <  1:1118235/235‑1 (MQ=255)
                                               gggagcGGTGGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTgg                                                                                                                                                                                             <  1:843150/218‑1 (MQ=255)
                                                           ttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTaca                                                                                                                                                                                                                     >  2:1013476‑M2/79‑183 (MQ=255)
                                                           ttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                                                                                                                                   <  1:1098194‑M2/123‑1 (MQ=255)
                                                            tGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgcg                                                                                                                                                                                                  <  2:1031696/201‑1 (MQ=255)
                                                              aTTCGGCTGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTa                                                                                                                                                               <  1:487629/234‑1 (MQ=255)
                                                                 gggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAg                                                                                                                                                              <  1:1013477‑M2/159‑1 (MQ=255)
                                                                  gggtGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgc                                                                                                                                                                                                     <  1:477333/189‑1 (MQ=255)
                                                                         ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc                                                                                                                                                                                                   <  2:799278‑M2/123‑1 (MQ=255)
                                                                             gAGGCCGTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACg                                                                                                                                                                                                                                         <  2:984018/145‑1 (MQ=255)
                                                                                 ccTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCGCGCGAGCCCCCTGGCGCACg                                                                                                                                                                                                                                         <  2:339862/141‑1 (MQ=255)
                                                                                       ggCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTg                                                                                                                                                                                              <  1:770800/178‑1 (MQ=255)
                                                                                              gCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTCCGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                                                                                                                                               <  1:237073/138‑1 (MQ=255)
                                                                                              gCCTTGAGGGTTGTAGTGCCTGTGAGATTAGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcg                                                                                                                                                                                                        <  1:317833/161‑1 (MQ=255)
                                                                                                         ttcagcgCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGaaacaa                                                                                                                                                                    <  1:952802/179‑1 (MQ=2)
                                                                                                                  tgggagatcgagcgccgcccgcgCGCCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTg                                                                                                                                                                                              <  1:391296‑M2/128‑1 (MQ=255)
                                                                                                                        aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                                                                                                                                               >  1:264607/1‑112 (MQ=255)
                                                                                                                        aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                                                                                                                                               <  2:264606/112‑1 (MQ=255)
                                                                                                                                              cTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCg                                                                                                                                                                                                                                    <  2:890695/85‑1 (MQ=255)
                                                                                                                                                                 ccTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaa                                                                                                                                                                        <  1:952/126‑1 (MQ=255)
                                                                                                                                                                   tACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACag                                                                                                                                                                                        >  2:29530/1‑108 (MQ=255)
                                                                                                                                                                             tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAACAACTATCGCCAGGCa                                                                                     >  2:930578/1‑197 (MQ=255)
                                                                                                                                                                             tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTcgc                                              >  1:983013/1‑236 (MQ=255)
                                                                                                                                                                               gggCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgat                                                                                                                    <  2:904564‑M2/164‑32 (MQ=255)
                                                                                                                                                                               gggCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGGTTACCCGCGAAgcgc                                                                                                                                       >  2:1058992/1‑144 (MQ=2)
                                                                                                                                                                                                  gatttgcgcgcgaccgCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtcccctcaa                                                         <  1:88567‑M2/198‑91 (MQ=255)
                                                                                                                                                                                                       gcgcgcGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaa                                                                                        >  1:453382‑M2/1‑108 (MQ=255)
                                                                                                                                                                                                       gcgcgcGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaa                                                                                        <  2:453381‑M2/167‑60 (MQ=255)
                                                                                                                                                                                                       gcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATGCTGCAATTTc                                                              >  1:910230/1‑194 (MQ=255)
                                                                                                                                                                                                         aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAAcc                    >  1:965376/1‑234 (MQ=255)
                                                                                                                                                                                                                  cccTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATc          >  1:708067/1‑235 (MQ=255)
                                                                                                                                                                                                                        cgcACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTcgccg                                            >  1:463318/1‑195 (MQ=255)
                                                                                                                                                                                                                         gaaCGCCTTGATAGCGCGTCGTAAAACCAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGACAGAAACAAACGTa                                                                                                                                                               <  2:929403/77‑1 (MQ=255)
                                                                                                                                                                                                                          cACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAAtt  >  1:1030782/1‑235 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        
TGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATT  >  AE015451/5988628‑5989080

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