Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,988,760 C→G 14.8% intergenic (+54/‑492) 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,7600CG14.8% 75.0 / 8.2 27intergenic (+54/‑492)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (7/16);  new base G (2/2);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-02

TGGGCCTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCG  >  AE015451/5988537‑5988957
                                                                                                                                                                                                                               |                                                                                                                                                                                                     
tGGGCCTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCT                                                                                                                                                                                            >  1:144251‑M1/1‑235 (MQ=255)
                           aCCGCATTGCACGGTTTACACAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAAcc                                                                                                                                                                               <  1:383082/221‑1 (MQ=255)
                                       ggTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAg                                                                                                                                                      <  2:144247/234‑1 (MQ=255)
                                                               tGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgac                                                                                                                             >  1:502744‑M1/1‑172 (MQ=255)
                                                               tGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAcgcg                                                                                                                              <  1:498098/234‑1 (MQ=255)
                                                                    gAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACgtttg                                                                                                                                                                 <  1:49612/194‑1 (MQ=255)
                                                                        cccAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                             <  1:399390/194‑1 (MQ=255)
                                                                                    tGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCaatccc                                                                                                                                                  <  1:256110/193‑6 (MQ=255)
                                                                                    tGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCaatccc                                                                                                                                                  >  2:256106/1‑188 (MQ=255)
                                                                                     gCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTaaa                                                                                                                                                                                 >  2:399386/1‑161 (MQ=255)
                                                                                          atgaCCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGc                                                                                                                                                                            <  2:409518/161‑1 (MQ=255)
                                                                                          atgaCCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATcgc                                                                                                  >  1:354410/1‑235 (MQ=255)
                                                                                                        gTGATGCAGACGGCCAGGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                             <  2:119679/162‑1 (MQ=255)
                                                                                                                          gaggCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                             <  2:378178/144‑1 (MQ=255)
                                                                                                                           aggCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGAc                                                                                                                                       <  1:320830/165‑1 (MQ=255)
                                                                                                                           aggCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGAc                                                                                                                                       >  2:320826/1‑165 (MQ=255)
                                                                                                                                                                        gAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCCCGCCTGGCACAGGAGGGATTGGCCAGa                                              <  1:329054/209‑1 (MQ=255)
                                                                                                                                                                              ttttGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcct                                                        <  2:430093‑M1/193‑133 (MQ=255)
                                                                                                                                                                                        ggCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGtt                          <  1:234742/212‑1 (MQ=255)
                                                                                                                                                                                          ccTTGAGGGTTGTAGTGGCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacc                                                                                                                      <  2:509154‑M1/119‑71 (MQ=255)
                                                                                                                                                                                                        gTGCCTGTGAGATTGAGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtctaattctcagctgggtcctacttccccgcgaaccatttcctacgttcacaccaatccccaaggaattatccta  >  1:193237‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                         tGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtcta                                                                         <  1:409538‑M1/149‑116 (MQ=255)
                                                                                                                                                                                                         tGCCTGTGAGATTGGGCGCCGCGCGGGCGGCGCTcgatctcatgggcgcccaaaaactatcgccaggcacctggaagccatcctgcaatttcctgacgatcacctcgccgccatgacgcgaactcacaaaaaccatcacaaatcgtcta                                                                         >  2:409534‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                           ccTGTGAGATTGGTCGCCCCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTc                                                                                                                                                             <  2:94469/63‑1 (MQ=255)
                                                                                                                                                                                                                gagaTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                   >  2:327706/1‑100 (MQ=255)
                                                                                                                                                                                                                gagaTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGgc                                                                                                                   <  1:327710/100‑1 (MQ=255)
                                                                                                                                                                                                                     tGGGCGTCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCGCCTACGTGTGTTTCGGGGCAGTTATTCATGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                    <  1:176478/110‑1 (MQ=25)
                                                                                                                                                                                                                               |                                                                                                                                                                                                     
TGGGCCTGAGCGATGAACAGGCGGCGGACCGCATTGCACGGTTTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCG  >  AE015451/5988537‑5988957

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