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%	172.9 / NA	53	A414T (GCC→ACC)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/27);  total (26/27)

AAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCA  >  NC_000913/3457490‑3457751                                                                                                                               |                                                                                                                                        aaaaaCGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAatggca                                                                                                                             <  1:1738652/139‑1 (MQ=255)         tGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCa                                                                                                                     <  1:309293/139‑1 (MQ=255)                          ccAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    <  2:496324/139‑1 (MQ=255)                          ccAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    <  1:68248/139‑1 (MQ=255)                           cAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    >  1:764979/1‑138 (MQ=255)                           cAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    <  2:764979/138‑1 (MQ=255)                             aTACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAACAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTg                                                                                                 >  2:1159842/1‑139 (MQ=255)                             aTACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTg                                                                                                 >  1:1149270/1‑139 (MQ=255)                                 cGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCGGGCTTCTTCAATGGCGACTGGGGa                                                                                             <  2:1149270/139‑1 (MQ=255)                                   gATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                                  <  2:668140/116‑1 (MQ=255)                                   gATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                                  >  1:668140/1‑116 (MQ=255)                                   gATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGt                                                                                           >  2:1625406/1‑139 (MQ=255)                                           aTTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTa                                                                                   <  1:1625406/139‑1 (MQ=255)                                                        cgcAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAg                                                                      >  1:1571653/1‑139 (MQ=255)                                                         gcAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                                     <  1:235565/139‑1 (MQ=255)                                                                  ggCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                               <  2:498101/104‑1 (MQ=255)                                                                  ggCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                               >  1:498101/1‑104 (MQ=255)                                                                  ggCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaa                                                            >  2:265684/1‑139 (MQ=255)                                                                  ggCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTAAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaa                                                            >  1:381286/1‑139 (MQ=255)                                                                         tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACa                                                                                                                                   <  2:501868/61‑1 (MQ=255)                                                                         tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACa                                                                                                                                   >  1:501868/1‑61 (MQ=255)                                                                         tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCgcgc                                                                             >  1:216339/1‑115 (MQ=255)                                                                         tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCgcgc                                                                             <  2:216339/115‑1 (MQ=255)                                                                         tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                        >  2:967312/1‑136 (MQ=255)                                                                         tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                        <  1:967312/136‑1 (MQ=255)                                                                                 gAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGcc                                             <  1:265684/139‑1 (MQ=255)                                                                                   aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    >  2:838356/1‑82 (MQ=255)                                                                                   aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    <  1:838356/82‑1 (MQ=255)                                                                                         ggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAg                                                                                            >  1:1212265/1‑84 (MQ=255)                                                                                         ggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAg                                                                                            <  2:1212265/84‑1 (MQ=255)                                                                                         ggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAg                                     >  2:1677831/1‑139 (MQ=255)                                                                                            tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                                  >  2:187695/1‑139 (MQ=255)                                                                                              aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCg                                >  1:1477193/1‑139 (MQ=255)                                                                                                 aGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACaa                                                                                                                                  <  1:1351786/38‑1 (MQ=255)                                                                                                 aGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACaa                                                                                                                                  >  2:1351786/1‑38 (MQ=255)                                                                                                 aGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTaa                             >  1:909342/1‑139 (MQ=255)                                                                                                   cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTg                                                                                                 <  2:1676841/69‑1 (MQ=255)                                                                                                   cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTg                                                                                                 >  1:1676841/1‑69 (MQ=255)                                                                                                   cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCgcgc                                                                             >  2:1734484/1‑89 (MQ=255)                                                                                                   cgcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCgcgc                                                                             <  1:1734484/89‑1 (MQ=255)                                                                                                   cgcgAATCCTGCCTGGGATATGTTAAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                           <  2:381286/139‑1 (MQ=255)                                                                                                      gAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACGTCCTCGCCACCCCAAGCATCGTAACGCTg                        <  2:783239/139‑1 (MQ=255)                                                                                                      gAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTg                        <  2:1477193/139‑1 (MQ=255)                                                                                                        aTCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCg                                               >  2:1068997/1‑114 (MQ=255)                                                                                                        aTCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCg                                               <  1:1068997/114‑1 (MQ=255)                                                                                                        aTCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGa                      >  1:1088216/1‑139 (MQ=255)                                                                                                           ctgcctgGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataa                   <  2:1088216/139‑1 (MQ=255)                                                                                                             gcctgGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataata                 >  2:1488357/1‑139 (MQ=255)                                                                                                                        gTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTcc      >  2:434288/1‑139 (MQ=255)                                                                                                                          ttAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTGACCTCGCTAGCCAGTAACAATAAAAATGACATCCTCGCGCCCCCAAGCATCCTAACGCTGGATAATAAACCCGCGTCCtt    <  1:1528490/139‑1 (MQ=255)                                                                                                                           tAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAAT‑AACTCGCGTCCTTCa  <  2:1571653/139‑1 (MQ=255)                                                                                                                            aGCACCTACAATGGCATGGCCGCAGGCTTCTTCAAt                                                                                                         <  2:512236/36‑1 (MQ=255)                                                                                                                            aGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGa                                                                                             >  1:503855/1‑48 (MQ=255)                                                                                                                               |                                                                                                                                        AAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCA  >  NC_000913/3457490‑3457751

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