Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,317,517 A→G I300T (ATC→ACC)  trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,317,5170AG88.5% 171.9 / 9.2 61I300T (ATC→ACC) trpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (27/27);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

CCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTC  >  NC_000913/1317385‑1317650
                                                                                                                                    |                                                                                                                                     
ccTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGt                                                                                                                                 >  1:1694527/1‑139 (MQ=255)
       gCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGAttc                                                                                                                          >  1:1633199/1‑139 (MQ=255)
       gCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGAttc                                                                                                                          >  2:1296684/1‑139 (MQ=255)
                gCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTg                                                                                                                 >  1:731117/1‑139 (MQ=255)
                gCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTg                                                                                                                 >  2:1289409/1‑139 (MQ=255)
                 cGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCa                                                                                                                      <  1:1099964/133‑1 (MQ=255)
                 cGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCa                                                                                                                      >  2:1099964/1‑133 (MQ=255)
                      tGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTc                                                                                                           >  1:1559725/1‑139 (MQ=255)
                        aaGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCtt                                                                                                         <  2:1694527/139‑1 (MQ=255)
                        aaGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCtt                                                                                                         <  2:1559725/139‑1 (MQ=255)
                             ttCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGAttc                                                                                                                          >  1:1348037/1‑117 (MQ=255)
                             ttCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGAttc                                                                                                                          <  2:1348037/117‑1 (MQ=255)
                                     cTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGt                                                                                                            >  2:857035/1‑123 (MQ=255)
                                     cTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGt                                                                                                            <  1:857035/123‑1 (MQ=255)
                                        catcatCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         <  2:938186/139‑1 (MQ=255)
                                             tcGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         >  1:696115/1‑134 (MQ=255)
                                             tcGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         <  2:696115/134‑1 (MQ=255)
                                                gTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCa                                                                                 <  1:1289409/139‑1 (MQ=255)
                                                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCa                                                                                                                      >  2:1355956/1‑91 (MQ=255)
                                                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCa                                                                                                                      <  1:1355956/91‑1 (MQ=255)
                                                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAg                                                                      >  2:860709/1‑139 (MQ=255)
                                                           gTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAg                                                                      >  2:952698/1‑139 (MQ=255)
                                                              aTCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGtt                                                                                                    <  2:1385762/106‑1 (MQ=255)
                                                              aTCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGtt                                                                                                    >  1:1385762/1‑106 (MQ=255)
                                                                  gcgcGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCa                                                               >  1:1341943/1‑139 (MQ=255)
                                                                       tCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGa                                                          <  1:952698/139‑1 (MQ=255)
                                                                             gCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGa                                                                           >  2:932905/1‑116 (MQ=255)
                                                                             gCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGa                                                                           <  1:932905/116‑1 (MQ=255)
                                                                                         cgcgTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCt                                                                                                          >  1:1294793/1‑73 (MQ=255)
                                                                                         cgcgTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCt                                                                                                          <  2:1294793/73‑1 (MQ=255)
                                                                                           cgTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCa                                                                                 >  2:1400099/1‑96 (MQ=255)
                                                                                           cgTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCa                                                                                 <  1:1400099/96‑1 (MQ=255)
                                                                                             tgttgtGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         <  1:1203765/86‑1 (MQ=255)
                                                                                             tgttgtGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         >  2:1203765/1‑86 (MQ=255)
                                                                                              gttgtGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTg                                   <  2:731117/139‑1 (MQ=255)
                                                                                                tgtgGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCa                                                                                 <  2:124543/91‑1 (MQ=255)
                                                                                                tgtgGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCa                                                                                 >  1:124543/1‑91 (MQ=255)
                                                                                                 gtgGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCg                                                                            >  2:189527/1‑95 (MQ=255)
                                                                                                 gtgGGCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCg                                                                            <  1:189527/95‑1 (MQ=255)
                                                                                                   gggCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         >  1:1520265/1‑80 (MQ=255)
                                                                                                   gggCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCgg                                                                                         <  2:1520265/80‑1 (MQ=255)
                                                                                                    ggCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGcc                             <  2:1633199/139‑1 (MQ=255)
                                                                                                    ggCCGACAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGcc                             >  1:1108936/1‑139 (MQ=255)
                                                                                                         aCAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATa                                                                               >  1:1317532/1‑84 (MQ=255)
                                                                                                         aCAGACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATa                                                                               <  2:1317532/84‑1 (MQ=255)
                                                                                                            gACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTg                                   <  1:1051941/125‑1 (MQ=255)
                                                                                                            gACGGGAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTg                                   >  2:1051941/1‑125 (MQ=255)
                                                                                                               gggAAATCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATAc                  <  2:102695/139‑1 (MQ=255)
                                                                                                                    aTCCAGTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGa             <  2:1108936/139‑1 (MQ=255)
                                                                                                                       cAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGa             >  2:311446/1‑136 (MQ=255)
                                                                                                                       cAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGa             <  1:311446/136‑1 (MQ=255)
                                                                                                                         gTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCttt                                                                                   <  2:53870/64‑1 (MQ=255)
                                                                                                                         gTCCGGCGGAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCttt                                                                                   >  1:53870/1‑64 (MQ=255)
                                                                                                                           ccGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCa                                                       <  2:1538790/90‑1 (MQ=255)
                                                                                                                           ccGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCa                                                       >  1:1538790/1‑90 (MQ=255)
                                                                                                                              gcggAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCt   <  1:1513687/139‑1 (MQ=255)
                                                                                                                               cggAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATAcc                                                                             <  2:246743/64‑1 (MQ=255)
                                                                                                                               cggAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATAcc                                                                             >  1:246743/1‑64 (MQ=255)
                                                                                                                               cggAGGTGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTc  >  2:1667586/1‑139 (MQ=255)
                                                                                                                               cggAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCg                                                                            <  1:928378/65‑1 (MQ=255)
                                                                                                                               cggAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCg                                                                            >  2:928378/1‑65 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATGGAGTAAGATTCTTCAATCTGCCCGTCTTCGGTTTGCATCATCGGCGCTTTCATACCGAAATAGATACCCACGCGACCATGTTTTAGCGGTGCGCCGTGCTCGCCAGTTTCGATACCGTGACCACCTGGCTC  >  NC_000913/1317385‑1317650

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