Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,981,890 G→A 8.5% V166V (GTC→GTT)  otsB ← trehalose‑6‑phosphate phosphatase, biosynthetic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,981,8900GA 8.5% 144.5 / 2.3 47V166V (GTC→GTT) otsBtrehalose‑6‑phosphate phosphatase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base G (26/17);  new base A (2/2);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGG  >  NC_000913/1981752‑1982036
                                                                                                                                          |                                                                                                                                                  
cAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAacacac                                                                                                                                           >  2:103428/1‑148 (MQ=255)
 aGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTACGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCAACAACACACt                                                                                                                                          >  1:308354/1‑148 (MQ=255)
    cGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTc                                                                                                                                       >  1:261233/1‑148 (MQ=255)
       ttAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCct                                                                                                                                    <  6:115216/148‑1 (MQ=255)
          aCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCctgc                                                                                                                                  <  8:306374/147‑1 (MQ=255)
                         gATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgtg                                                                                                                  >  4:165445/1‑148 (MQ=255)
                          aTTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTgg                                                                                                                 <  7:335202/148‑1 (MQ=255)
                           ttCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgt                                                                                                                   <  8:214215/145‑1 (MQ=255)
                           ttCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGATTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgt                                                                                                                   >  7:214215/1‑145 (MQ=255)
                                 ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGAt                                                                                                           <  3:165445/147‑1 (MQ=255)
                                  gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCt                                                                                                         >  4:296511/1‑148 (MQ=255)
                                  gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGGGTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCt                                                                                                         >  2:141550/1‑148 (MQ=255)
                                        tcatcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTaa                                                                                                   <  7:229584/148‑1 (MQ=255)
                                        tcatcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTaa                                                                                                   >  3:254133/1‑148 (MQ=255)
                                              cccAGAAATACGGGCGTTAGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGtt                                                                                             >  6:132426/1‑148 (MQ=255)
                                               ccAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACCCACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACTTTg                                                                                            >  1:281600/1‑148 (MQ=255)
                                                   aaaTACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACATTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgcgc                                                                                         >  7:102569/1‑147 (MQ=255)
                                                       aCGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAg                                                                                                      <  7:27551/130‑1 (MQ=255)
                                                       aCGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCAACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATg                                                                                    >  3:145524/1‑148 (MQ=255)
                                                       aCGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACATTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAg                                                                                                      >  8:27551/1‑130 (MQ=255)
                                                              ttCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTaa                                                                             <  1:103428/148‑1 (MQ=255)
                                                                         aaGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCa                                                                   <  1:67741/147‑1 (MQ=255)
                                                                                      gCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTc                                                                       >  5:97971/1‑130 (MQ=255)
                                                                                      gCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTc                                                                       <  6:97971/130‑1 (MQ=255)
                                                                                          aaaaGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGa                                                 >  7:89555/1‑148 (MQ=255)
                                                                                             aGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATaa                                              >  4:136936/1‑148 (MQ=255)
                                                                                                      tGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAgcgc                                      >  4:194787/1‑147 (MQ=255)
                                                                                                        ccTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATc                                                                                                          <  8:174941/77‑1 (MQ=255)
                                                                                                        ccTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATc                                                                                                          >  7:174941/1‑77 (MQ=255)
                                                                                                        ccTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCaaa                                   <  5:132426/148‑1 (MQ=255)
                                                                                                          tCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAg                                 >  2:217857/1‑148 (MQ=255)
                                                                                                           cACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGc                                >  1:252188/1‑148 (MQ=255)
                                                                                                            aCCTTTACTGGTACCTCTCGGTTTGATCTCAACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGc                                <  4:145524/147‑1 (MQ=255)
                                                                                                               tttACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATc                            >  3:213948/1‑148 (MQ=255)
                                                                                                                      gTACCTCTCGGTTTGATCTCAACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGAAGATAATGCAGCGCAAAAGCCATCCCttt                       <  2:308354/146‑1 (MQ=255)
                                                                                                                         ccTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCt                  >  5:331/1‑148 (MQ=255)
                                                                                                                          ctctCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTc                 <  2:261233/148‑1 (MQ=255)
                                                                                                                          ctctCGGTTTGATCTCGACAACACACTTTACCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACg                                                                                               <  5:308830/70‑1 (MQ=255)
                                                                                                                            ctcGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTc                 >  4:22152/1‑146 (MQ=255)
                                                                                                                             tcGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCTAg              >  2:54662/1‑148 (MQ=255)
                                                                                                                              cGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGc             >  7:47126/1‑148 (MQ=255)
                                                                                                                               ggTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGc             <  8:102569/147‑1 (MQ=255)
                                                                                                                                    gATCTCGACAACAAACTTTCCATGCTGTAACGCCATTTGTTGCAAGATCTGAGTAATACGTTGCTCTaa                                                                                      <  6:288755/69‑1 (MQ=255)
                                                                                                                                     aTCTCGACAACACACTTTCCCTGCTGTAACGCAATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAAt                                                                                     >  5:288755/1‑69 (MQ=255)
                                                                                                                                      tctcGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAg                                         >  7:147967/1‑112 (MQ=255)
                                                                                                                                      tctcGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAg                                         <  8:147967/112‑1 (MQ=255)
                                                                                                                                       ctcGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCg    >  6:292658/1‑148 (MQ=255)
                                                                                                                                          gACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCggg  <  3:136936/147‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATGCGTCTTCATGCTGCGGAGCCTGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGG  >  NC_000913/1981752‑1982036

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: ‑