Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,787,410	A→G	T318A (ACT→GCT)	aroH →	3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,787,410	0	A	G	100.0%	195.4 / NA	62	T318A (ACT→GCT)	aroH	3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/52);  total (10/52)

CCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGC  >  NC_000913/1787322‑1787499                                                                                         |                                                                                          ccGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTa                                                                                       <  1:518502/93‑1 (MQ=255)  cGCAATGGCTCTACGGCGACTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTAc                                                                                      >  2:4566272/1‑93 (MQ=255)    cAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACgg                                                                                    <  1:994170/93‑1 (MQ=255)         gCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAAt                                                                               <  1:3474744/93‑1 (MQ=255)          ctctACGGCGATTGCTGAAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATc                                                                              <  1:3461896/93‑1 (MQ=255)           tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                             >  2:3133405/1‑93 (MQ=255)           tctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATcc                                                                             <  1:854717/93‑1 (MQ=255)            ctACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCa                                                                            <  1:4629132/93‑1 (MQ=255)              aCGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCAtt                                                                          <  2:4553788/93‑1 (MQ=255)                ggCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAc                                                                        <  2:3289724/93‑1 (MQ=255)                 gCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAcc                                                                       <  2:3986553/93‑1 (MQ=255)                 gCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTAcc                                                                       <  2:2963706/93‑1 (MQ=255)                    aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                    <  2:3655169/93‑1 (MQ=255)                    aTTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAc                                                                    <  2:3500127/93‑1 (MQ=255)                      tGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGAccc                                                                  <  2:2081972/93‑1 (MQ=255)                        cTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgt                                                                <  2:4357116/93‑1 (MQ=255)                         tGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtg                                                               >  1:109089/1‑93 (MQ=255)                         tGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCgtg                                                               <  2:2229960/93‑1 (MQ=255)                           gAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTc                                                             <  2:745141/93‑1 (MQ=255)                           gAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTc                                                             <  2:2295375/93‑1 (MQ=255)                           gAATTATGGAGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTc                                                             <  1:2853428/93‑1 (MQ=255)                           gAATTATGACTGAAAGTTGCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTc                                                             <  2:1908445/93‑1 (MQ=255)                            aaTTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTct                                                            >  1:354245/1‑93 (MQ=255)                             aTTAGGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctg                                                           <  2:3384891/93‑1 (MQ=255)                              ttATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgg                                                          <  2:1467688/93‑1 (MQ=255)                                 tGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggct                                                       >  2:3433986/1‑93 (MQ=255)                                 tGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTctgggct                                                       <  2:458956/93‑1 (MQ=255)                                   gcgcAAAGTTGCCTGCGCGAAGAAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGGCAATCCATTACCGACCCGTGTctgggctgg                                                     <  2:570814/93‑1 (MQ=255)                                     ggAAAGTTTCCTGGGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGgga                                                   <  2:1531324/93‑1 (MQ=255)                                      gAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGggag                                                  >  2:687339/1‑93 (MQ=255)                                         aGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGAt                                               <  2:1749236/93‑1 (MQ=255)                                         aGTTTCCCGCGCGAAGGAACGCAAAAAATCGTCGGCCGTCCTCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGAt                                               <  1:1159809/93‑1 (MQ=255)                                           tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  1:4430877/93‑1 (MQ=255)                                           tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  1:4275260/93‑1 (MQ=255)                                           tttCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  1:969412/93‑1 (MQ=255)                                           tttCCTGCGCAAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATAc                                             <  2:4101461/93‑1 (MQ=255)                                               cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                         <  2:2412860/93‑1 (MQ=255)                                               cTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGaa                                         >  2:2194133/1‑93 (MQ=255)                                                tGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAAc                                        <  1:2325997/93‑1 (MQ=255)                                                tGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAAc                                        <  2:4294106/93‑1 (MQ=255)                                                   gcgAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGcc                                     <  1:2680927/93‑1 (MQ=255)                                                     gAAGGAACGCAAAAGATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  2:4076625/93‑1 (MQ=255)                                                     gAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTg                                   <  1:4389831/93‑1 (MQ=255)                                                      aaGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTgg                                  >  1:268680/1‑93 (MQ=255)                                                       aGGAACGCAAAAAATCGTTGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGGCTGGGCTGGGAGGATACCGAACGCCTGGt                                 >  1:4363650/1‑93 (MQ=255)                                                       aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCTACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 <  2:608420/93‑1 (MQ=255)                                                       aGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGt                                 <  2:4543960/93‑1 (MQ=255)                                                             gCAAAAAATCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGaaaa                           >  2:210804/1‑93 (MQ=255)                                                               aaaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAAc                         <  2:2463071/93‑1 (MQ=255)                                                                aaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACt                        <  2:1861617/93‑1 (MQ=255)                                                                aaaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACt                        <  2:3148044/93‑1 (MQ=255)                                                                 aaaaTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTc                       <  2:1306979/93‑1 (MQ=255)                                                                    aTCGTCGGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGcc                    <  2:3922488/93‑1 (MQ=255)                                                                     tcgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGctt                   <  1:4531403/93‑3 (MQ=255)                                                                      cgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  2:2435644/93‑1 (MQ=255)                                                                      cgtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  1:1753906/93‑1 (MQ=255)                                                                      cgtcgGCAGTCAGCCACTCGTTGCCGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctc                  <  2:2235077/93‑1 (MQ=255)                                                                       gtcgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCctct                 <  2:1148279/93‑1 (MQ=255)                                                                         cgGCAGTCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGc               <  1:1495032/93‑1 (MQ=255)                                                                               tCAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGa         <  1:1641436/93‑1 (MQ=255)                                                                                cAGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGAt        <  2:5030/93‑1 (MQ=255)                                                                                 aGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATa       <  1:938784/93‑1 (MQ=255)                                                                                 aGCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATa       <  2:4134409/93‑1 (MQ=255)                                                                                  gCCGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAc      <  1:2229403/93‑1 (MQ=255)                                                                                   gcgctcgctTACGGTCAATCCATTACCGACCCGTGTCTGTGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAcc     <  1:1490739/92‑1 (MQ=255)                                                                                   ccGCTCGCTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAcc     <  1:894676/93‑1 (MQ=255)                                                                                    cgctcgctTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATAccc    <  2:1658250/93‑1 (MQ=255)                                                                                      cTCGCTTACGATCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGc  <  2:1926639/93‑1 (MQ=255)                                                                                         |                                                                                          CCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAACGCAAAAAATCGTCGGCAGTCAGCCGCTCACTTACGGTCAATCCATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGC  >  NC_000913/1787322‑1787499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑