Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 264,098 G→T L50L (CTC→CTA ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA264,0980GT96.4% 74.9 / ‑4.9 28L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (15/12);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_pae_tpiA/263963‑264243
                                                                                                                                       |                                                                                                                                                 
gcgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATGGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGa                                                                                                                                      <  2:9084/149‑1 (MQ=255)
          aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                            <  2:227700/149‑1 (MQ=255)
           tctcCCGCTGCGGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTGCGGCGATTTCTGGGACg                                                                                                                           >  1:266693/1‑149 (MQ=38)
                    gCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCGCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                  <  1:304650/149‑1 (MQ=255)
                                tAGGGTTGTGCATTACAGGTCTCTCCTTTTGTCAGTTCAGGGTGTGGCTGATCAGGCAGGTATAGGGCCAATTGCGGTCGTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAg                                                                                                      <  2:363797/149‑1 (MQ=21)
                                        cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                           >  2:55309/1‑136 (MQ=255)
                                        cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                           <  1:55309/136‑1 (MQ=255)
                                         tCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCg                                                                                             >  2:399441/1‑149 (MQ=255)
                                                    ctcCACCATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                  <  2:325277/149‑1 (MQ=255)
                                                        aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTg                                                                              >  2:19710/1‑149 (MQ=255)
                                                         ctcttGTCAGTTCTGGGTGATGCTCATCACGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACTTCCCTGCCGGTGTAATGAGAGACGCTACCGTGCCACTGGTATTTGc                                                                             <  2:77308/146‑1 (MQ=21)
                                                              gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                                        >  1:419188/1‑149 (MQ=255)
                                                                                tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCCCTGGTATTTGCCGGGGCAGTAGAGAAAGATACgg                                                      >  1:264815/1‑149 (MQ=25)
                                                                                tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATCTCCGGAACGTCCCTGCCGGTGGAATGACAGACGCTACCGTGCCCCTGGTATTTGCCGGTGCACTCGCGAAAGATCCgg                                                      >  2:2048/1‑149 (MQ=255)
                                                                                                gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                      >  2:22697/1‑149 (MQ=255)
                                                                                                 ggCCCATGGCGGTCCTGGCGGCGTTCGGCTTATACCGCTAGAACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                     <  2:266693/149‑1 (MQ=255)
                                                                                                   ccaaTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCCTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGt                                   <  1:85441/146‑1 (MQ=255)
                                                                                                             tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                                    <  1:64357/90‑1 (MQ=255)
                                                                                                             tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                                    >  2:64357/1‑90 (MQ=255)
                                                                                                                tGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGTCAAAGATTCGGGCATAACTATTCGTGCGGTATATAGTCATTGc                      >  2:175183/1‑149 (MQ=17)
                                                                                                                 ggcggcGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGcc                     >  1:30533/1‑149 (MQ=255)
                                                                                                                  gcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTCGCGAAAGATTCGGGGCTCAGGATGCTGTCGGTATATCGCTATTGccc                    >  2:49628/1‑149 (MQ=14)
                                                                                                                  gcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCggg                                                     >  2:61514/1‑116 (MQ=255)
                                                                                                                  gcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCggg                                                     <  1:61514/116‑1 (MQ=255)
                                                                                                                            gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTAAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTaat          >  1:301278/1‑147 (MQ=255)
                                                                                                                              gTAAACCGCTAGGAGTCCGGCAATGTCCGGAACGTCCTGGCCGGTGTAATGAGAGACGCTAGCGTGCCACTGGTATTTGCCAGTGCAGTAGCGAAAGATTCGGCACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGa        <  2:282825/149‑1 (MQ=14)
                                                                                                                                ataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATa      >  1:103125/1‑149 (MQ=255)
                                                                                                                                   ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGGATATCGTCATTGCCCTCCGTTTACTTATAAtt   >  1:326994/1‑149 (MQ=255)
                                                                                                                                    cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAttt  >  1:340511/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
GCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_pae_tpiA/263963‑264243

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