Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,456,471	T→G	F32C (TTC→TGC)	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,456,471	0	T	G	100.0%	178.7 / NA	47	F32C (TTC→TGC)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/29);  total (18/29)

CCGCTTGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTTCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTACCAGTTCTTTTTA  >  NC_000913/3456324‑3456616                                                                                                                                                    |                                                                                                                                                  ccGCTTGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCaa                                                                                                                                                <  1:626051/151‑1 (MQ=255) ccGCTTGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCaa                                                                                                                                                <  1:593065/151‑1 (MQ=255)   gCTTGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAAt                                                                                                                                               >  1:451265/1‑150 (MQ=255)    cTTGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATa                                                                                                                                              >  2:516648/1‑150 (MQ=255)      tGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATaa                                                                                                                                             <  1:793310/149‑1 (MQ=255)      tGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACg                                                                                                                                           >  1:230952/1‑151 (MQ=255)       gATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACg                                                                                                                                           <  2:230947/150‑1 (MQ=255)          aaaTGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCg                                                                                                                                        <  1:231135/150‑1 (MQ=255)           aaTGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGat                                                                                                                                      <  2:176022/151‑1 (MQ=255)                ttCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCGCCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATcc                                                                                                                                  <  2:488395/150‑1 (MQ=255)                  ccGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGcc                                                                                                                               <  2:806712/151‑1 (MQ=255)                  ccGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGc                                                                                                                                <  2:591739/150‑1 (MQ=255)                     tCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAg                                                                                                                             >  1:586913/1‑150 (MQ=255)                        gggAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCg                                                                                                                         <  1:579063/151‑1 (MQ=255)                        gggAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCg                                                                                                                         >  2:286840/1‑151 (MQ=255)                                        gaatgGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCa                                                                                                         <  1:7060/151‑1 (MQ=255)                                         aatgGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATTCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCa                                                                                                         <  2:594362/150‑1 (MQ=255)                                             gACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATc                                                                                                    <  2:554166/151‑1 (MQ=255)                                                  cGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGc                                                                                               >  2:216702/1‑151 (MQ=255)                                                     cATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCaaa                                                                                            <  1:461446/151‑1 (MQ=255)                                                        gAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACg                                                                                         <  2:689078/151‑1 (MQ=255)                                                         aaaGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGa                                                                                        <  2:40274/151‑1 (MQ=255)                                                                tCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGAt                                                                                 >  1:829289/1‑151 (MQ=255)                                                                tCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGAt                                                                                 >  2:63552/1‑151 (MQ=255)                                                                 cAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATc                                                                                >  1:372464/1‑151 (MQ=255)                                                                     aaaaTCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGAcc                                                                            <  2:469335/151‑1 (MQ=255)                                                                               gctgctTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTAc                                                                   >  2:830326/1‑150 (MQ=255)                                                                                  tcttgctGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAggg                                                               >  1:50970/2‑151 (MQ=255)                                                                                     tgctGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGa                                                              <  2:240788/149‑1 (MQ=255)                                                                                         ggCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCAtt                                                        <  1:903849/151‑1 (MQ=255)                                                                                          gcagcaCTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCAttt                                                       <  2:151472/151‑1 (MQ=255)                                                                                              caCTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGt                                                   <  1:244802/151‑1 (MQ=255)                                                                                                 tactCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACg                                                >  2:420564/1‑151 (MQ=255)                                                                                                   ctCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGc                                               >  2:39400/1‑150 (MQ=255)                                                                                                     cATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCa                                              <  1:425541/149‑1 (MQ=255)                                                                                                       tGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTa                                           <  2:329016/150‑1 (MQ=255)                                                                                                         ccTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATg                                        <  2:417841/151‑1 (MQ=255)                                                                                                             gtgCAGGACACGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATAc                                    >  1:867514/1‑151 (MQ=255)                                                                                                                    acacGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGc                             >  1:118680/1‑151 (MQ=255)                                                                                                                      acGCTGAGAACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGcc                            >  2:626207/1‑150 (MQ=255)                                                                                                                           gagaACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAgg                      >  2:895990/1‑151 (MQ=255)                                                                                                                           gagaACGAACAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAg                       <  2:218466/150‑1 (MQ=255)                                                                                                                                    cAATACGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTAcc             >  1:110883/1‑151 (MQ=255)                                                                                                                                        aCGGCGCGAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTACCAGt          <  1:441261/150‑1 (MQ=255)                                                                                                                                            cgcgAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTACCAGTTCttt     <  2:586898/151‑1 (MQ=255)                                                                                                                                              cgAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTACCAGTTCttttt   <  2:919773/151‑1 (MQ=255)                                                                                                                                               gAACTGCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTACCAGTTCTTTTTa  <  1:423676/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCGCTTGATAAATGTTTCCGTCGGGGAACTTACAGGAATGAATGGACTGCGTCATGAAAGGACTCAATAAAATCACCTGCTGCTTGCTGGCAGCACTACTCATGCCTTGTGCAGGACACGCTGAGAACGAACAATACGGCGCGAACTTCAATAACGCCGATATCCGCCAGTTCGTGGAAATAGTGGGTCAGCATCTTGGCAAAACGATCCTGATCGACCCTTCGGTACAGGGAACCATTTCCGTACGCAGTAATGATACGTTTAGCCAACAGGAGTACTACCAGTTCTTTTTA  >  NC_000913/3456324‑3456616

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