Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,874,078	C→T	*35* (TAG→TAA)	yoaI ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,874,078	0	C	T	100.0%	172.0 / NA	54	*35* (TAG→TAA)	yoaI	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/31);  total (23/31)

CTCTACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAACTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAC  >  NC_000913/1873947‑1874212                                                                                                                                    |                                                                                                                                       ctctACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCgg                                                                                                                                 >  1:250167/1‑139 (MQ=255)         ttCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCg                                                                                                                         >  1:1182862/1‑139 (MQ=255)          tCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGa                                                                                                                        >  2:1261413/1‑139 (MQ=255)          tCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGa                                                                                                                        >  1:57007/1‑139 (MQ=255)          tCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGa                                                                                                                        >  1:830277/1‑139 (MQ=255)           cGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGAt                                                                                                                       >  2:45908/1‑139 (MQ=255)           cGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCG‑CGCACCGCGTAATTAGCCGGTCCGGTCGATa                                                                                                                      >  1:1062474/1‑139 (MQ=255)                 aTTTGATGAAATCGCCTTCACTGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGt                                                                                                                                >  1:890829/1‑124 (MQ=255)                 aTTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGt                                                                                                                                <  2:890829/124‑1 (MQ=255)                   ttGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAgaagaa                                                                                                               <  2:250167/139‑1 (MQ=255)                             cGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAc                                                                                                     <  1:438760/139‑1 (MQ=255)                                          ggCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACaaa                                                                                                       >  2:1292541/1‑124 (MQ=255)                                          ggCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACaaa                                                                                                       <  1:1292541/124‑1 (MQ=255)                                          ggCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCa                                                                                        >  1:1058279/1‑139 (MQ=255)                                              tGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGt                                                                                                                            >  2:144031/1‑99 (MQ=255)                                              tGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGt                                                                                                                            <  1:144031/99‑1 (MQ=255)                                                cTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCa                                                                                        >  2:221784/1‑133 (MQ=255)                                                cTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCa                                                                                        <  1:221784/133‑1 (MQ=255)                                                        aTCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGa                                                                          >  1:633189/1‑139 (MQ=255)                                                          ccGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATg                                                                        <  2:1045305/139‑1 (MQ=255)                                                          ccGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATg                                                                        <  1:45908/139‑1 (MQ=255)                                                          ccGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATg                                                                        <  1:1261413/139‑1 (MQ=255)                                                           cGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGa                                                                       <  1:695588/139‑1 (MQ=255)                                                            gTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGAc                                                                      <  1:519389/139‑1 (MQ=255)                                                            gTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGAc                                                                      <  2:830277/139‑1 (MQ=255)                                                            gTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGAc                                                                      <  1:1036765/139‑1 (MQ=255)                                                             ttATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAc                                                                                                     >  2:459203/1‑107 (MQ=255)                                                             ttATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAc                                                                                                     <  1:459203/107‑1 (MQ=255)                                                                        aCCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAgt                                                          <  2:57007/139‑1 (MQ=255)                                                                          cgcgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAgtgt                                                        >  1:213595/1‑139 (MQ=255)                                                                          cgcgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAgtgt                                                        <  2:773303/139‑1 (MQ=255)                                                                           gcgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGtc                                                       >  2:255723/1‑139 (MQ=255)                                                                            cgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGAc                                                                                                          <  1:877421/87‑1 (MQ=255)                                                                            cgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGAc                                                                                                          >  2:877421/1‑87 (MQ=255)                                                                                         cattgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGTGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTg                                         >  2:993270/1‑139 (MQ=255)                                                                                         aattgccaaACCA‑AACTGGGCAAATTCCCCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGa                                        <  2:1062474/138‑1 (MQ=255)                                                                                          attgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGa                                        <  2:213595/139‑1 (MQ=255)                                                                                          attgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGa                                        <  1:1182410/139‑1 (MQ=255)                                                                                           ttgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAcc                                                                                                    <  1:822933/78‑1 (MQ=255)                                                                                           ttgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAcc                                                                                                    >  2:822933/1‑78 (MQ=255)                                                                                            tgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATc                                      >  1:161423/1‑139 (MQ=255)                                                                                                 aaCCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATaa                                                                                                                     >  1:799077/1‑55 (MQ=255)                                                                                                 aaCCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATaa                                                                                                                     <  2:799077/55‑1 (MQ=255)                                                                                                      gAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGAc                                                                      >  1:665881/1‑97 (MQ=255)                                                                                                      gAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGAc                                                                      <  2:665881/97‑1 (MQ=255)                                                                                                            gTCAATTTCCGCCGCACTGTGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGcccc                      <  1:993270/139‑1 (MQ=255)                                                                                                            gTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGcccc                      <  2:1182862/139‑1 (MQ=255)                                                                                                              cAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCg                    >  1:1008594/1‑139 (MQ=255)                                                                                                               aaTTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCgg                   <  2:510131/139‑1 (MQ=255)                                                                                                                 tttCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGAc                                                  >  2:150747/1‑106 (MQ=255)                                                                                                                 tttCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGAc                                                  <  1:150747/106‑1 (MQ=255)                                                                                                                      gccgcACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAg            <  1:255723/139‑1 (MQ=255)                                                                                                                      gccgcACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAg            <  2:633189/139‑1 (MQ=255)                                                                                                                      gccgcACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAg            <  2:161423/139‑1 (MQ=255)                                                                                                                                gTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAc  <  1:1243042/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       CTCTACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAACTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAC  >  NC_000913/1873947‑1874212

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