Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 460,284 G→T 28.2% A725S (GCA→TCA)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb460,2840GT28.2% 90.7 / 19.9 39A725S (GCA→TCA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base T (5/6);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGA  >  W3110S.gb/460217‑460348
                                                                   |                                                                
aaTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAg                                                               >  1:3456859/1‑71 (MQ=255)
                 ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcgc                                                       <  1:4893941/62‑1 (MQ=255)
                  tctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAgcgc                                                            <  1:4891845/56‑1 (MQ=255)
                  tctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCgcgg                                                     >  1:922959/1‑63 (MQ=255)
                   ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCgcgg                                                     >  1:3737840/1‑62 (MQ=255)
                   ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCgcgg                                                     >  1:2012005/1‑62 (MQ=255)
                   ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTa                                             >  1:342181/1‑70 (MQ=255)
                      cGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTaaaa                                          <  1:3371227/70‑1 (MQ=255)
                          gTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGt                                      <  1:1604257/70‑1 (MQ=255)
                             aGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTaa                                            >  1:3355987/1‑61 (MQ=255)
                                  cTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTcc                             >  1:2129965/1‑71 (MQ=255)
                                  cTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTcc                             >  1:3000488/1‑71 (MQ=255)
                                    gCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGt                           >  1:3355774/1‑71 (MQ=255)
                                     ccGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGtt                          >  1:2653575/1‑71 (MQ=255)
                                     ccGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGtt                          >  1:3440684/1‑71 (MQ=255)
                                        aTCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcgc                                                       <  1:2794313/39‑1 (MQ=255)
                                           ggtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGa                                                >  1:3830199/1‑43 (MQ=255)
                                             tggtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATa                  <  1:531189/71‑1 (MQ=255)
                                             tggtTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATa                  <  1:2794255/71‑1 (MQ=255)
                                             tggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTaa                                            <  1:1364034/45‑1 (MQ=255)
                                              ggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATa                  >  1:1075213/1‑70 (MQ=255)
                                                tttGAAAGAAAAACTCCTGGCAGCGCATCGcggcgg                                                  <  1:2690217/36‑1 (MQ=255)
                                                tttGAAAGAAAAACTCCTGGCAGCGCATCGcggcgg                                                  <  1:1647966/36‑1 (MQ=255)
                                                 ttGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAAc               >  1:4523831/1‑70 (MQ=255)
                                                  tGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAAAACAGTGCt                                   <  1:4657831/49‑1 (MQ=255)
                                                  tGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAAcgc             >  1:3502901/1‑71 (MQ=255)
                                                  tGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAAcgc             >  1:2714958/1‑71 (MQ=255)
                                                   gaaagaaaAACTCCTGTCAGCGCATCGCGGCGGGATTaaaa                                          <  1:4069815/41‑1 (MQ=255)
                                                   gaaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTaa                                            <  1:4796835/39‑1 (MQ=255)
                                                    aaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAg                                       >  1:3760302/1‑43 (MQ=255)
                                                    aaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGa           >  1:3583879/1‑71 (MQ=255)
                                                     aagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTaaaa                                          <  1:814093/39‑1 (MQ=255)
                                                     aagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTaaaa                                          <  1:1705121/39‑1 (MQ=255)
                                                     aagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTaa                                            <  1:3389364/37‑1 (MQ=255)
                                                      agaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATc         >  1:1046063/1‑71 (MQ=255)
                                                       gaaaAACTCCTGTCAGCGCATCGCGGCGGGATTaaaa                                          >  1:776515/1‑37 (MQ=255)
                                                       gaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCt                                   <  1:2265084/44‑1 (MQ=255)
                                                             cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAaga  <  1:294418/71‑1 (MQ=255)
                                                              tCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGaa                      <  1:579285/50‑1 (MQ=255)
                                                                   |                                                                
AATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGA  >  W3110S.gb/460217‑460348

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