Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,437,911 T→A 13.2% intergenic (+17/‑113) trkA → / → mscL NAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,437,9110TA13.2% 110.8 / 8.6 38intergenic (+17/‑113)trkA/mscLNAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base A (2/3);  total (18/20)
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.92e-01

ACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTT  >  NC_000913/3437770‑3438053
                                                                                                                                             |                                                                                                                                              
acgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAATAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCt                                                                                                                                         <  4:269083/149‑1 (MQ=255)
 cgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtt                                                                                                                                        <  2:394962/149‑1 (MQ=255)
     atgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCCACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCttatta                                                                                                                                    <  3:185649/149‑1 (MQ=255)
        atTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTAttt                                                                                                                                 <  3:130125/149‑1 (MQ=255)
                tGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACCTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCtttta                                                                                                                         <  4:342562/149‑3 (MQ=255)
                           gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                       >  4:425294/1‑124 (MQ=255)
                           gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                       <  3:425294/124‑1 (MQ=255)
                                gAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTGATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatt                                                                                                         >  1:421298/1‑149 (MQ=255)
                                         gATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattcatt                                                                                                >  2:243409/1‑149 (MQ=255)
                                            cACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCt                                                                                             <  4:307325/149‑1 (MQ=255)
                                                           cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                        <  2:147019/123‑1 (MQ=255)
                                                           cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                        >  1:147019/1‑123 (MQ=255)
                                                              acaGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTaa                                                                           >  2:315735/1‑149 (MQ=255)
                                                                      aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAAt                                                                                  <  1:339621/134‑1 (MQ=255)
                                                                      aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAAt                                                                                  >  2:339621/1‑134 (MQ=255)
                                                                       aaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTGATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                  <  2:421298/149‑1 (MQ=255)
                                                                       aaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                  <  3:308443/149‑1 (MQ=255)
                                                                                  cGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                  >  4:239275/1‑138 (MQ=255)
                                                                                  cGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                  <  3:239275/138‑1 (MQ=255)
                                                                                       tCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGc                                                                                                                                           <  3:440040/60‑1 (MQ=255)
                                                                                       tCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGc                                                                                                                                           >  4:440040/1‑60 (MQ=255)
                                                                                           aaGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTc                                              <  1:243409/149‑1 (MQ=255)
                                                                                                tctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                         >  1:204483/1‑149 (MQ=255)
                                                                                                                  tcttctTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCatt                       <  1:315735/149‑1 (MQ=255)
                                                                                                                      cttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCTCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCAttatta                   >  4:276684/1‑149 (MQ=255)
                                                                                                                       ttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                     <  2:101675/66‑1 (MQ=255)
                                                                                                                       ttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                     >  1:101675/1‑66 (MQ=255)
                                                                                                                        tGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTGAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTaaa                 <  3:276684/149‑1 (MQ=255)
                                                                                                                          taattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGa               >  2:254594/1‑149 (MQ=255)
                                                                                                                               aataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc          <  2:204483/149‑1 (MQ=255)
                                                                                                                                ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATg                                                          >  2:324393/1‑100 (MQ=255)
                                                                                                                                ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATg                                                          <  1:324393/100‑1 (MQ=255)
                                                                                                                                ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAAGATGAGCATTATTAGAGAATTTcg         >  3:298094/1‑149 (MQ=255)
                                                                                                                                 taaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgc        >  3:337810/1‑149 (MQ=255)
                                                                                                                                 taaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc          <  4:9986/147‑1 (MQ=255)
                                                                                                                                 taaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc          >  3:9986/1‑147 (MQ=255)
                                                                                                                                  aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGATTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGCATTTcgcg       >  1:11898/1‑149 (MQ=255)
                                                                                                                                       gTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCGTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAAttt  >  2:375165/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                              
ACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTT  >  NC_000913/3437770‑3438053

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