Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,988,889 C→A 27.2% intergenic (+183/‑363) 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,8890CA27.2% 58.5 / 19.6 33intergenic (+183/‑363)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (14/10);  new base A (3/6);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCG  >  AE015451/5988705‑5989103
                                                                                                                                                                                        |                                                                                                                                                                                                                      
gAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCCCGCCTGGCACAGGAGGGATTGGCCAGa                                                                                                                                                                                                <  1:329054/209‑1 (MQ=255)
                ggCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGtt                                                                                                                                                                            <  1:234742/212‑1 (MQ=255)
                     ggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTggccga                                                                                                                                                                                                 <  2:505559‑M2/148‑3 (MQ=255)
                                        gagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTg                                                                                                                                                                                                                     >  1:300811‑M2/21‑148 (MQ=255)
                                        gagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTg                                                                                                                                                                                                                     <  2:300807‑M2/128‑1 (MQ=255)
                                                gcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAtt                                                                                                                                                                                                       <  2:446633‑M2/142‑1 (MQ=255)
                                                gcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAtt                                                                                                                                                                                                       >  1:446637‑M2/13‑154 (MQ=255)
                                                  gccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGa                                                                                                                                                                                                <  1:435474‑M2/149‑1 (MQ=255)
                                                                                     cTACGTTTGTTTCTGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCt                                                                                                                                                                                                                      >  2:205147/1‑102 (MQ=255)
                                                                                     cTACGTTTGTTTCTGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCt                                                                                                                                                                                                                      <  1:205151/102‑1 (MQ=255)
                                                                                                tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatct                                                                                                                                         <  2:296451‑M2/168‑34 (MQ=255)
                                                                                                tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggc                                                                                                                                                  >  2:472032‑M2/1‑135 (MQ=255)
                                                                                                tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggc                                                                                                                                                  <  1:472036‑M2/159‑25 (MQ=255)
                                                                                                                     gACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccg                                                                                                                                                                      <  1:9521‑M2/118‑5 (MQ=255)
                                                                                                                      aCGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACaaa                                               >  1:129182/1‑236 (MQ=255)
                                                                                                                            aCGCGAGCCCCCTGTCGCACGCCTCGATATCGCGGCGTACAAACAAGGCCGTCGCGCGCGCCTgg                                                                                                                                                                                                                    >  2:68961/1‑65 (MQ=255)
                                                                                                                                             catgtctcgatATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaaca                        <  1:441568‑M2/230‑147 (MQ=255)
                                                                                                                                                  cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTgg                                                                                                       <  1:289836/152‑1 (MQ=255)
                                                                                                                                                  cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTgg                                                                                                       >  2:289832/1‑152 (MQ=255)
                                                                                                                                                   tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGCAGGGATTGGcc                                                                                                                                                                                                   >  2:19136/1‑59 (MQ=255)
                                                                                                                                                   tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggg                                                                                                                                  <  1:90020‑M2/123‑41 (MQ=255)
                                                                                                                                                     gATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGa                                                                                                      >  1:32735/1‑150 (MQ=255)
                                                                                                                                                          cgcgTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTa                                                                                                                                                                                      <  2:96556/64‑1 (MQ=14)
                                                                                                                                                               cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACATACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGTCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAg                                                                                                               >  1:198526/1‑131 (MQ=255)
                                                                                                                                                               cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTc      >  1:403450/1‑236 (MQ=255)
                                                                                                                                                                  acaaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaat                                                       >  1:408609‑M2/1‑68 (MQ=255)
                                                                                                                                                                    aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCcg  >  1:435321/1‑235 (MQ=255)
                                                                                                                                                                              gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggc                                                                                                                                                  >  1:517030‑M2/1‑57 (MQ=255)
                                                                                                                                                                              gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggc                                                                                                                                                  <  2:517026‑M2/81‑25 (MQ=255)
                                                                                                                                                                              gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAAtt                         >  2:142697/1‑202 (MQ=255)
                                                                                                                                                                              gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATGTCCTGACGATCAc                                                                          >  2:88597/1‑153 (MQ=255)
                                                                                                                                                                              gTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccggcca                                                                                     <  2:408605‑M2/141‑86 (MQ=255)
                                                                                                                                                                              gTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCcgc                                                                                                                                                                     >  2:146275/1‑61 (MQ=18)
                                                                                                                                                                                        |                                                                                                                                                                                                                      
GAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCG  >  AE015451/5988705‑5989103

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