Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,463,952 G→A S253F (TCC→TTC)  treC ← trehalose‑6‑phosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,463,9520GA87.9% 75.9 / 0.5 33S253F (TCC→TTC) treCtrehalose‑6‑phosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (14/15);  total (16/17)
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.59e-01

ACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGC  >  NC_000913/4463708‑4464172
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                            
aCCAAAGCGAGAAACAATGCGCGGCTGATCGTTGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTAAAGGCTGGTGGAGAACAttt                                                                                                                                                                                                                        <  1:306681/251‑1 (MQ=255)
                                ggTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAAAAATGTTTTCAAAGCCACAAAGTCAGGGTTAGCCAGCGTCCCTTTTTCACCACCGGGAAAAACGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACCGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAac                                                                                                                                                                                        >  1:183903/1‑251 (MQ=255)
                                     cacCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCATGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATc                                                                                                                                                                                   <  1:156448/251‑1 (MQ=255)
                                           aaCAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCcgcgatgtgtttacgccacgcgggt                                                                                                                                                                             <  1:326455/251‑22 (MQ=255)
                                             cAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACCCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGGAAACACATCGCGGTTcc                                                                                                                                                                           >  1:166071/1‑250 (MQ=255)
                                              aaGGCATTCCATGCTACGTTATGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGATTAGCCAGCGTCCATTTTTCACCACCGGGCTAATCGACCTTCAGGTGATGAAAATTAAAGGTCAACGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTTGTGGAGGACATTTCACCTACGGTCCTTAACCCCCGTGGCGTAAACACATCCCGGTTCAt                                                                                                                                                                          >  1:60589/1‑251 (MQ=255)
                                               aGGCATTCCATGCTTCGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCTCGTGGCGTAAACACATCGCGGTTCAtc                                                                                                                                                                         <  1:232199/251‑1 (MQ=255)
                                                   aTTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGACAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTAACCTACGGTCATTAACCCGCGGGGCGTAAACACATCGCGGTTCATCTCgg                                                                                                                                                                     >  1:381477/1‑250 (MQ=255)
                                                     tCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTg                                                                                                                                                                    <  1:196554/250‑1 (MQ=255)
                                                      ccATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCa                                                                                                                                                                  <  1:352494/251‑1 (MQ=255)
                                                                              ttgCCAGTGGCGGAACAATGTTTTCAACCCCACAAAATCAGGGTTAGCCAGCGTCCATTTTTCACCACCCGGAAAATCGACCTTCAGGGGATGACAATTAAAGGTCATCGACAATTCCCTGCCCGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCAGGTGGAGGACATTTCACCTACGGTCCTTAACCCGCCTGGCGTAAACACATCCCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCgc                                                                                                                                          >  1:68489/1‑250 (MQ=255)
                                                                                gCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCg                                                                                                                                        >  1:110075/1‑251 (MQ=255)
                                                                                         ggAACAATGTTTTCAACGCCACAAAGTGAGATTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTAGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGGGTAGa                                                                                                                                 >  1:84766/1‑249 (MQ=255)
                                                                                                        aCGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCg                                                                                                                <  1:113254/251‑1 (MQ=255)
                                                                                                           caacaAATTCACGTTTAGCCAGCGTCCATGTTGCACCACCGGGATAATCGACCGTTAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTTGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTACCCCGCGTGGCGTAAACACATCGCGGTTCAGCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTAGGTGTATAAGCGACGCCCGTCGCCGTcc                                                                                                             <  2:127870/249‑1 (MQ=255)
                                                                                                               aaaGTTAGGTTTAGCCAGCGTCAATTTTTCACCACCGGGATAATCGACATTCAGGTGATGAATATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATAGCTGGCAATGCTCCAGGCTGGTGGCGAACATTTCCACTACGGTCTTTAACCCGCGTGTCCTAAACACATCGCTGTTCATCTCGTGCATAAACTAGTGTGCTCTTAGCCCGTCGGTCGAGTAGCGACgtcatttcacggcgacgt                                                                                                         >  2:101191/1‑234 (MQ=255)
                                                                                                                              cAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAAcgcg                                                                                          >  1:206143/1‑251 (MQ=255)
                                                                                                                                 cGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGCCAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTACAGGTCTTCAGGGAAACGAGGGt                                                                                       >  1:372480/1‑251 (MQ=255)
                                                                                                                                                                       gatgaAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAAccc                                                  <  2:110075/250‑1 (MQ=255)
                                                                                                                                                                       gatgaAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCg                                                 >  1:239131/1‑251 (MQ=255)
                                                                                                                                                                             aatgaaaGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGAcc                                           <  1:245284/247‑1 (MQ=255)
                                                                                                                                                                             aaTTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGAcc                                           <  1:336863/251‑1 (MQ=255)
                                                                                                                                                                                    ggTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGAt                                    >  1:191738/1‑251 (MQ=255)
                                                                                                                                                                                           gACAATGCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGTCTGGTGGAGAACAGTTCACCTACGGTCATTAACCCGCGTGGCGTAGACACATCGCGGTTCATCTCGTGCAAAAACTCGTTTGCTCGTGGCCCGTCGGTGTAGATGCGACGCCCGTCGCCGTCCAGGTCTTCAGTGAAACGGGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCCACCTCACGAGCGGCCCa                             <  2:166071/251‑1 (MQ=255)
                                                                                                                                                                                                ttCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACt                        >  1:278463/1‑251 (MQ=255)
                                                                                                                                                                                                  cACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACACCACGATCGGCCCAGAACTc                       >  1:20798/1‑250 (MQ=255)
                                                                                                                                                                                                      gCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTTGAGAACATTTCTCCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAATCTCGTGGGCACGTGTCCCGTCGGTTTAGAATCGATGCCCGTCGCCGTCCAg                                                                                                           <  2:280333/162‑1 (MQ=255)
                                                                                                                                                                                                      gCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAg                                                                                                           >  1:280333/1‑162 (MQ=255)
                                                                                                                                                                                                       ccTGTCAGAGCGTCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCAACTACGGTCATTAACCCGGGGGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCAAAGAc                 <  2:60589/251‑1 (MQ=255)
                                                                                                                                                                                                            cAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACtttt             <  1:36826/250‑1 (MQ=255)
                                                                                                                                                                                                                    cgTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGctct    <  1:155908/251‑1 (MQ=255)
                                                                                                                                                                                                                    cgTAGCGCTGGCGATGCTCAAGGCTGGTGGAGAACATTGCGCCTACGGTCATTAACCCGCGTGGCGTAAACTCATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCCTCAGGGAAACGCGGGTCTTTGTAGATCAGATTCACCACATCCAGGCGCATCCCGTCGACCCCACGATCGGCTCAGAACTCACAGACTTTTTTCAGctct    <  2:278463/251‑1 (MQ=255)
                                                                                                                                                                                                                      tATCGCTGGCAATGCTCAAGGCTGGAGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGCGGCCCGTCGGGGCAGAAGCGACGCCCGTCGCCGTCCAGGCCTTCAGGGAAACGCGGGTCTTTGGAGCTCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAATACACAGCCTTTTTTCAGCTCTgc  >  1:39130/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                            
ACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGC  >  NC_000913/4463708‑4464172

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