Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,317,648 C→A 11.2% E256D (GAG→GAT)  trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,317,6480CA11.2% 156.5 / 3.2 56E256D (GAG→GAT) trpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base A (2/5);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGC  >  NC_000913/1317505‑1317790
                                                                                                                                               |                                                                                                                                              
aGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCcac                                                                                                                                           <  7:341310/149‑1 (MQ=255)
aGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCcac                                                                                                                                           <  7:255860/149‑1 (MQ=255)
 gTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCcaca                                                                                                                                          <  5:165868/149‑1 (MQ=255)
               gAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCCCACCAATCAgg                                                                                                                                 >  3:393061/1‑144 (MQ=255)
                aGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTcc                                                                                                                                             >  5:69834/1‑131 (MQ=255)
                aGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTcc                                                                                                                                             <  6:69834/131‑1 (MQ=255)
                    aGATTCTTCTCGCAGCCCGTCCTCTGTTTGACTCATCGGCGCTTTCATACCTAAATCGAAACCCACGCTACCATGTTTTAGAGTTTAGATTTTTTTTTTTTTTTTGATACGGTGACCACCTTGATCTACACCAATCAGGCCGACGTTg                                                                                                                        <  2:497084/148‑1 (MQ=255)
                          ttcAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCAGTTTTTTTTTTTTTTGATACGGTGACCACCTTGATCTTCACCAATCAGGc                                                                                                                                <  2:113584/134‑1 (MQ=255)
                          ttcAATCTGCCCGTCTTCGGTATGCATCATCGGCGCTTTCAAACAGAAATAGATACCCACGCGACCATGTATTAGAGATGCGCCGTGAACGACAGTTTCGATACCGTGACCACCAGTCTCCACACCAATCAGGc                                                                                                                                >  1:113584/1‑134 (MQ=255)
                             aaTCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACACACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGATCCACACCAATCAGGCAGACGTTGGTTTCATtga                                                                                                              >  3:115963/1‑149 (MQ=255)
                              aTCTGACCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGAACATGTTTTAGCGCTGCGCCGTGCTCGCCAGTTTCGAAACCGTGACCACCTGGATCCACACCAATCAGGCCGACGTAGGTTTCaaagat                                                                                                             >  2:122932/1‑144 (MQ=255)
                                          tCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCaaa                                                                                                 >  2:97965/1‑149 (MQ=255)
                                             gTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACAt                                                                                              >  6:465834/1‑149 (MQ=255)
                                                   atcatcGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGAt                                                                                        >  4:188975/1‑149 (MQ=255)
                                                    tcatcGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGt                                                                                                                      >  6:151369/1‑118 (MQ=255)
                                                    tcatcGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGt                                                                                                                      <  5:151369/118‑1 (MQ=255)
                                                    tcatcGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATg                                                                                       >  6:344161/1‑149 (MQ=255)
                                                    tcatcGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATg                                                                                       >  7:104377/1‑149 (MQ=255)
                                                                 cATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgc                                                                          >  4:149076/1‑149 (MQ=255)
                                                                  aTACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGTTTTTTTTTTTTGATACGGTGACCACCTTGATCTACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgcc                                                                         <  4:115963/149‑1 (MQ=255)
                                                                   tACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgccg                                                                        >  2:240825/1‑149 (MQ=255)
                                                                    aCCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgc                                                                          <  7:381585/146‑1 (MQ=255)
                                                                    aCCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgc                                                                          >  8:381585/1‑146 (MQ=255)
                                                                     ccGAAATAGATACCCACGCGACCATGTTTTAGCTGTGCGCCGTTCTCGCCAGTTTCGATACCGTGACCACCTGTCTACACACCAATCAGGCAGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACcgccgc                                                                       <  1:122932/148‑1 (MQ=255)
                                                                       gAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTTTTTTTTTTTTTTGATACGGTGACCACCTTGATCTACACCAAt                                                                                                                                     <  4:130457/84‑1 (MQ=255)
                                                                       gAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAAt                                                                                                                                     >  3:130457/1‑84 (MQ=255)
                                                                        aaaTAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                                                    >  4:189875/1‑148 (MQ=255)
                                                                         aaTAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAaca                                                                  <  1:97965/149‑1 (MQ=255)
                                                                              aTACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGc                                                             <  7:274041/149‑1 (MQ=255)
                                                                                     cgcgACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACg                                                      <  3:189875/149‑1 (MQ=255)
                                                                                        gACCATGTTTTAGCGGTGCGCCGTGTTTTTTATTTTTGATACGGTGACCACCTTGATCTACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgcc                                                                         <  2:90344/127‑1 (MQ=255)
                                                                                        gACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAAccgcc                                                                         >  1:90344/1‑127 (MQ=255)
                                                                                        gACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCa                                                   <  5:344161/149‑1 (MQ=255)
                                                                                        gACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACAAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCa                                                   >  1:317544/1‑149 (MQ=255)
                                                                                          ccATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATc                                                 >  6:12191/1‑149 (MQ=255)
                                                                                               ttttAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCgg                                              >  5:144583/1‑147 (MQ=255)
                                                                                                 ttAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAg                                           <  3:459616/148‑1 (MQ=255)
                                                                                                            ccGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTctct                               <  3:188975/149‑1 (MQ=255)
                                                                                                                gCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCTGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTc                            >  3:302787/1‑148 (MQ=255)
                                                                                                                  tCGTCAGTTTCGATACGGTGACCACCTTGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAg                                                                                                     <  6:105845/73‑1 (MQ=255)
                                                                                                                      cAGTTTCGATACCGTGACCACCTTGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAac                                                                   <  4:263442/103‑1 (MQ=255)
                                                                                                                      cAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAac                                                                   >  3:263442/1‑103 (MQ=255)
                                                                                                                         tttCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTgc                  >  6:461920/1‑149 (MQ=255)
                                                                                                                           tCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTgcgc                >  2:324223/1‑149 (MQ=255)
                                                                                                                             gATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCtt              >  2:205947/1‑149 (MQ=255)
                                                                                                                               tACCGTGACCACCTGGCTCCACACAAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTg            <  2:317544/149‑1 (MQ=255)
                                                                                                                                  cGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGtt         <  7:480796/149‑1 (MQ=255)
                                                                                                                                  cGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGtt         >  8:396491/1‑149 (MQ=255)
                                                                                                                                    tGACCACCTGGCTCCACACCAATCAGGACGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGAATTCGAAACGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCAAAAATCTGAGCTATGGTatc       >  8:434652/1‑149 (MQ=255)
                                                                                                                                       ccaccTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTttcttc    <  1:324223/149‑1 (MQ=255)
                                                                                                                                         accTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGc  <  3:149076/149‑1 (MQ=255)
                                                                                                                                         accTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGAAGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGc  >  7:116050/1‑149 (MQ=255)
                                                                                                                                           cTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACg                                                      >  2:161644/1‑95 (MQ=255)
                                                                                                                                           cTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACg                                                      <  1:161644/95‑1 (MQ=255)
                                                                                                                                             ggCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGt          <  7:432244/137‑1 (MQ=255)
                                                                                                                                             ggCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGt          >  8:432244/1‑137 (MQ=255)
                                                                                                                                               |                                                                                                                                              
AGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGC  >  NC_000913/1317505‑1317790

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