Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,317,518 C→T 7.8% V507I (GTC→ATC)  treSA ← trehalose synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,317,5180CT 7.8% 144.2 / 2.7 51V507I (GTC→ATC) treSAtrehalose synthase A
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

GCACAGAAGCCTCTGCCTGTGGGTCAAGCCCGGCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCGGTGGTGCAGGCAATGCCGT  >  NC_002947/3317371‑3317661
                                                                                                                                                     |                                                                                                                                               
gCACAGAAGCCTCTGCCTGTGGGTCAAGCCCGGCCAGGTCATAGCCGCCCCTGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑Gacc                                                                                                                                                <  8:342846/149‑1 (MQ=255)
    aGAAGCCTCTGCCTGTGGGTCAAGCCCGGCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCAccc                                                                                                                                            >  6:224684/1‑149 (MQ=255)
       aGCCTCTGCCTGTGGGTCAAGCCCGGCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCg                                                                                                                                         >  4:86404/1‑149 (MQ=255)
                       tCAAGCCCGGCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTa                                                                                                                         <  8:285745/149‑1 (MQ=255)
                           gCCCGGCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACAt                                                                                                                      >  5:29943/1‑148 (MQ=255)
                               ggCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGAc                                                                                                                   >  1:272926/1‑147 (MQ=255)
                                   aGGTCATAGCCGCCCCGGTGAATAAAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTACGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCCCCCCCGGCTGCATGGCGTTGTAAATGACcagcag                                                                                                             >  2:36207/1‑149 (MQ=255)
                                      tCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACcagcagcag                                                                                                          >  5:283849/1‑149 (MQ=255)
                                      tCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACcagcagcag                                                                                                          <  3:86404/149‑1 (MQ=255)
                                       aaTACCCCCCCCGGTGAATCCAGAGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCATCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGCC‑GACCACCCCAGCCTGCATCGCGTTGTACATGAc                                                                                                                   >  1:176127/2‑139 (MQ=255)
                                           gccgccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGc                                                                                                      <  5:224684/148‑1 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GATCACCCCCGGCTGCATGGCGTTGTACATGACcagcag                                                                                                             >  7:303304/1‑138 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GATCACCCCCGGCTGCATGGCGTTGTACATGACcagcag                                                                                                             <  4:186465/138‑1 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GATCACCCCCGGCTGCATGGCGTTGTACATGACcagcag                                                                                                             >  3:186465/1‑138 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GATCACCCCCGGCTGCATGGCGTTGTACATGACcagcag                                                                                                             <  8:303304/138‑1 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGc                                                                                                                               >  1:111164/1‑120 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGc                                                                                                                               <  2:111164/120‑1 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTg                                                                                                       >  7:21296/1‑144 (MQ=255)
                                              gccCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTg                                                                                                       <  8:21296/144‑1 (MQ=255)
                                                cccTGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCTGCGGGCAAGGGCAGTG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATTGCTTTGTACATGACCAGCAGCAGGTGCACcttc                                                                                                <  1:36207/149‑1 (MQ=255)
                                                cccGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACcttc                                                                                                <  1:145460/149‑1 (MQ=255)
                                                  cGGTGAATCCAGCGGGGATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGGTGTACATGACCAGCAGCAGGTGCACcttctt                                                                                              <  2:46058/149‑1 (MQ=255)
                                                            aGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTc                                                                                     >  6:101332/1‑148 (MQ=255)
                                                               gggTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACa                                                                                                                       <  2:163274/111‑1 (MQ=255)
                                                                   aTCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCg                                                                             >  5:92892/1‑149 (MQ=255)
                                                                      gCCATCGCTCATGCGTTCGGGAACCGGCTCGGGGGGCAAGGGCATGG‑GAACGACCAAGTCCCATCCGGTAAAGGG‑GTCCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCt                                                                           <  7:59375/148‑1 (MQ=255)
                                                                       ccATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTac                                                                         >  7:170045/1‑149 (MQ=255)
                                                                       ccATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTa                                                                          >  5:310209/1‑148 (MQ=255)
                                                                       ccATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTa                                                                          >  5:315480/1‑148 (MQ=255)
                                                                        cATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTaca                                                                        >  2:120059/1‑149 (MQ=255)
                                                                         aTCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGGCCACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATc                                                                               >  1:269363/1‑142 (MQ=255)
                                                                         aTCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATc                                                                               <  2:269363/141‑1 (MQ=255)
                                                                         aTCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTacac                                                                       >  2:44703/1‑149 (MQ=255)
                                                                                 tGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGt                                                               <  1:44703/149‑1 (MQ=255)
                                                                                   cGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTc                                                             <  8:170045/149‑1 (MQ=255)
                                                                                   cGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTc                                                             <  2:272926/149‑1 (MQ=255)
                                                                                              cGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGaa                                                  <  6:283849/149‑1 (MQ=255)
                                                                                              cGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGAACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGaa                                                  >  4:253912/1‑149 (MQ=255)
                                                                                                      cggGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGCGGACCACCCCCgg                                                                                                                                        >  1:264100/1‑56 (MQ=255)
                                                                                                      cggGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCgg                                                                                                                                        <  2:264100/55‑1 (MQ=255)
                                                                                                             gggcaggt‑caCCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGGTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGcagcag                                   <  7:133182/149‑1 (MQ=255)
                                                                                                                      caCCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAgg                           <  5:101332/149‑1 (MQ=255)
                                                                                                                       aCCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGc                          >  8:355814/1‑149 (MQ=255)
                                                                                                                               ggTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCggtg                  <  6:315480/149‑1 (MQ=255)
                                                                                                                               ggTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCggtg                  <  6:310209/149‑1 (MQ=255)
                                                                                                                               ggTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCggtg                  >  1:260106/1‑149 (MQ=255)
                                                                                                                                     aGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCGGTGGTGCAg            <  1:120059/149‑1 (MQ=255)
                                                                                                                                       ccGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTgg                       <  4:10928/136‑1 (MQ=255)
                                                                                                                                       ccGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTgg                       >  3:10928/1‑136 (MQ=255)
                                                                                                                                               ggC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTgcagca                                    >  7:176972/1‑115 (MQ=255)
                                                                                                                                               ggC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTgcagca                                    <  8:176972/115‑1 (MQ=255)
                                                                                                                                               ggC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCGGTGGTGCAGGCAATGCCGt  >  7:102644/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCACAGAAGCCTCTGCCTGTGGGTCAAGCCCGGCCAGGTCATAGCCGCCCCGGTGAATCCAGCGGGTATCGCCATCGAGCATGCGTTCGGCAACCGCTTCGGCGGGCAAGGGCAGGG‑CACCGACCAGGTCCCAGCCGGACAAGGC‑GACCACCCCCGGCTGCATGGCGTTGTACATGACCAGCAGCAGGTGCACCTTCTTGATCAGTTCGATATCCGCTACACCAATCTGTTCGAGGTCGCGAATACCCAGTGCAGCAGCGATCAGGCTGGCGGTGGTGCAGGCAATGCCGT  >  NC_002947/3317371‑3317661

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