Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,457,616	G→A	A414T (GCC→ACC)	gspD →	general secretory pathway component, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,457,616	0	G	A	100.0%	193.4 / NA	60	A414T (GCC→ACC)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/31);  total (29/31)

GGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCG  >  NC_000913/3457483‑3457744                                                                                                                                      |                                                                                                                                 ggCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGTAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCGTGCCTGGGATATGTTTAGCACCTAc                                                                                                                             <  2:1441536/139‑1 (MQ=255) ggCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTAc                                                                                                                             >  1:1322614/1‑139 (MQ=255) ggCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTAc                                                                                                                             >  1:1441536/1‑139 (MQ=255)                   cgcACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                           <  1:513916/139‑1 (MQ=255)                   cgcACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                           <  1:268232/139‑1 (MQ=255)                   cgcACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGACTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                           <  2:1153680/139‑1 (MQ=255)                            aTTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAAt                                                                                                  >  2:1299325/1‑139 (MQ=255)                            aTTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAAt                                                                                                  <  2:1322614/139‑1 (MQ=255)                             tttACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATg                                                                                                 >  2:925245/1‑139 (MQ=255)                             tttACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATg                                                                                                 >  2:257012/1‑139 (MQ=255)                               tACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGc                                                                                               <  1:925245/139‑1 (MQ=255)                                    aTACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTg                                                                                          >  1:924960/1‑139 (MQ=255)                                      aCCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  1:1320319/139‑1 (MQ=255)                                      aCCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  1:257012/139‑1 (MQ=255)                                          gATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                           >  1:790077/1‑116 (MQ=255)                                          gATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                           <  2:790077/116‑1 (MQ=255)                                           aTTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTa                                                                                   >  2:1340078/1‑139 (MQ=255)                                             tACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACt                                                                                 <  1:1340078/139‑1 (MQ=255)                                                  aTTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTa                                                                            <  1:1299325/139‑1 (MQ=255)                                                             tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAatggca                                                                                                                      <  1:43268/86‑1 (MQ=255)                                                             tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAatggca                                                                                                                      >  2:43268/1‑86 (MQ=255)                                                                         ggCTGATTATAAAAAGAATGGTGGGATCACCAGCTCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGc                                                                                               <  2:985344/97‑1 (MQ=255)                                                                         ggCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGc                                                                                               >  1:985344/1‑97 (MQ=255)                                                                             gATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                 <  2:924960/139‑1 (MQ=255)                                                                              aTTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                             >  2:255092/1‑94 (MQ=255)                                                                              aTTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                             <  1:255092/94‑1 (MQ=255)                                                                              aTTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGa                                                >  1:483705/1‑139 (MQ=255)                                                                               ttATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACcgcg                                                                       <  1:280419/115‑1 (MQ=255)                                                                               ttATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACcgcg                                                                       >  2:280419/1‑115 (MQ=255)                                                                                tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACa                                              >  1:951700/1‑139 (MQ=255)                                                                                                ggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        >  1:974534/1‑81 (MQ=255)                                                                                                ggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAg                              >  2:697427/1‑139 (MQ=255)                                                                                                ggATCACCAGCGCGAATCCTGCCTGGGATATATTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  2:974534/81‑1 (MQ=255)                                                                                                   tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                              >  2:476418/1‑104 (MQ=255)                                                                                                   tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                              <  1:476418/104‑1 (MQ=255)                                                                                                   tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                           >  1:700170/1‑139 (MQ=255)                                                                                                     aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCg                         >  2:575709/1‑139 (MQ=255)                                                                                                     aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCg                         >  1:379510/1‑139 (MQ=255)                                                                                                        aGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        >  2:872622/1‑73 (MQ=255)                                                                                                        aGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  1:872622/73‑1 (MQ=255)                                                                                                        aGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTCCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTaa                      <  2:378417/139‑1 (MQ=255)                                                                                                          cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        >  2:1016125/1‑71 (MQ=255)                                                                                                          cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  1:1016125/71‑1 (MQ=255)                                                                                                          cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                    <  1:697427/139‑1 (MQ=255)                                                                                                          cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                    <  2:483705/139‑1 (MQ=255)                                                                                                           gcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGGTGGCCAGTAACAATAAAAATGACATCTTCGCCACCCCAAGCAGCGTAACGc                   <  2:700170/139‑1 (MQ=255)                                                                                                           gcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGc                   <  2:951700/139‑1 (MQ=255)                                                                                                            cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGcttct                                                                                                       >  2:950403/1‑54 (MQ=255)                                                                                                            cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGcttct                                                                                                       <  1:950403/54‑1 (MQ=255)                                                                                                            cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  1:300742/69‑1 (MQ=255)                                                                                                            cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        >  2:300742/1‑69 (MQ=255)                                                                                                            cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                    <  1:225973/137‑1 (MQ=255)                                                                                                            cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                    >  2:225973/1‑137 (MQ=255)                                                                                                                    gcctgGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        >  2:781650/1‑61 (MQ=255)                                                                                                                    gcctgGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                        <  1:781650/61‑1 (MQ=255)                                                                                                                    gcctgGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTAGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataata          >  1:1154302/1‑139 (MQ=255)                                                                                                                         ggATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTc     >  1:405887/1‑139 (MQ=255)                                                                                                                            tatGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTAGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgcg  <  2:1154302/139‑1 (MQ=255)                                                                                                                               gTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGAc                                               >  2:1276815/1‑91 (MQ=255)                                                                                                                               gTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGAc                                               <  1:1276815/91‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 GGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCG  >  NC_000913/3457483‑3457744

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