Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,457,324	T→C	D316D (GAT→GAC)	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,324	0	T	C	100.0%	139.9 / NA	39	D316D (GAT→GAC)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

TGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCT  >  NC_000913/3457194‑3457470                                                                                                                                   |                                                                                                                                                   tGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGt                                                                                                                                >  2:554526/1‑151 (MQ=255)  gAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTc                                                                                                                               <  1:533408/151‑1 (MQ=255)     atatGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCAtt                                                                                                                            >  1:507893/1‑151 (MQ=255)        tGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTAcc                                                                                                                         >  1:192068/1‑151 (MQ=255)            aaaGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGc                                                                                                                       <  2:114110/149‑1 (MQ=255)            aaaGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCt                                                                                                                      <  2:396002/150‑1 (MQ=255)             aaGCCACGAATCTGGTGGAAGTGCTAACCGTTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGa                                                                                                                    <  2:424695/151‑1 (MQ=255)                      aTCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAgtccgt                                                                                                           >  1:345271/1‑151 (MQ=255)                             ggAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGaa                                                                                                    <  1:290091/151‑1 (MQ=255)                                   gCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTc                                                                                              >  2:342247/1‑151 (MQ=255)                                    cTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCg                                                                                             >  1:671216/1‑151 (MQ=255)                                                 ccGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTg                                                                                 <  2:590534/150‑1 (MQ=255)                                                    aaaaGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGcgcg                                                                             >  2:581614/1‑151 (MQ=255)                                                        gCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTg                                                                         >  1:39097/1‑151 (MQ=255)                                                         cTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTgg                                                                        >  2:18583/1‑151 (MQ=255)                                                                aTGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCg                                                                 <  1:38703/151‑1 (MQ=255)                                                                aTGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCg                                                                 >  1:280370/1‑151 (MQ=255)                                                                   aaaaaGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCg                                                              >  1:415896/1‑151 (MQ=255)                                                                    aaaaGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGt                                                             >  1:354956/1‑151 (MQ=255)                                                                     aaaGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTg                                                            <  1:92332/151‑1 (MQ=255)                                                                            aTGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGt                                                     >  1:308133/1‑151 (MQ=255)                                                                             tGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTg                                                    >  2:336924/1‑151 (MQ=255)                                                                                   tAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGtt                                              >  1:314633/1‑151 (MQ=255)                                                                                    aaGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGCGCTGGTTg                                             >  2:10377/1‑151 (MQ=255)                                                                                     aGCCCTCCTCTTCTGGCGCGATGGATACCGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGa                                            <  1:176906/151‑1 (MQ=255)                                                                                        cctcctCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGc                                         >  1:342038/1‑151 (MQ=255)                                                                                           cctcTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAat                                      <  1:554542/151‑1 (MQ=255)                                                                                                  tGGCGCGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGtt                               >  1:188749/1‑151 (MQ=255)                                                                                                    gcgcGATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTg                              >  1:668843/1‑150 (MQ=255)                                                                                                      gcgATGGATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAg                           <  1:360203/151‑1 (MQ=255)                                                                                                           ggATAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAg                      <  2:140573/151‑1 (MQ=255)                                                                                                             aTAACGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGa                    <  1:500604/151‑1 (MQ=255)                                                                                                                 cGTCGCCATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGa                <  1:99466/151‑1 (MQ=255)                                                                                                                       cATTACCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGa          <  2:463309/151‑1 (MQ=255)                                                                                                                           aCCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTaa      <  1:493315/151‑1 (MQ=255)                                                                                                                           aCCGCCGACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTa       >  1:657410/1‑150 (MQ=255)                                                                                                                             cgccgACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAAc    <  1:58408/151‑1 (MQ=255)                                                                                                                               ccgACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCt  >  2:242306/1‑151 (MQ=255)                                                                                                                               ccgACGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCt  <  2:449330/151‑1 (MQ=255)                                                                                                                                   |                                                                                                                                                   TGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCGATGGATAACGTCGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCT  >  NC_000913/3457194‑3457470

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: ‑