Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,377,215	A→G	L69P (CTG→CCG)	sspA ←	stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,377,215	0	A	G	97.8%	146.6 / NA	45	L69P (CTG→CCG)	sspA	stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/22);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAG  >  NC_000913/3377068‑3377361                                                                                                                                                    |                                                                                                                                                   atgatgGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACgggg                                                                                                                                                 >  2:116986/1‑151 (MQ=255)    atgGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCa                                                                                                                                              <  1:428850/151‑1 (MQ=255)          ttCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCAc                                                                                                                                        >  2:4071/1‑151 (MQ=255)           tcatcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCAc                                                                                                                                        >  1:441950/1‑150 (MQ=255)           tcatcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCAc                                                                                                                                        <  2:24498/150‑1 (MQ=255)                   gtgtACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATccacca                                                                                                                               <  2:454658/151‑1 (MQ=255)                         cAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATTCGAGATTCCCACGGGGTCAGCTCACGATCCACCCGGGTCg                                                                                                                         >  1:226402/1‑151 (MQ=255)                         cAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCg                                                                                                                         <  2:297340/151‑1 (MQ=255)                             cTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGaa                                                                                                                      >  1:426142/1‑150 (MQ=255)                              tttttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAAc                                                                                                                     >  2:484434/1‑150 (MQ=255)                               ttttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAAc                                                                                                                     <  2:394298/149‑1 (MQ=255)                                            gcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGtt                                                                                                      >  2:536212/1‑151 (MQ=255)                                            gcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGtt                                                                                                      <  1:182476/151‑1 (MQ=255)                                              aTGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCCTCCAGATATTCCATAATGATGCGCGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCggg                                                                                                        >  1:416241/1‑147 (MQ=255)                                                gTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAg                                                                                                   <  2:437937/150‑1 (MQ=255)                                                  aCAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTc                                                                                                >  2:94842/1‑151 (MQ=255)                                                   cAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCa                                                                                               <  2:647064/151‑1 (MQ=255)                                                                 cgcgAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAgg                                                                                 <  2:300853/151‑1 (MQ=255)                                                                    gAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                                                              >  2:200419/1‑151 (MQ=255)                                                                         aCCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                           <  2:492133/149‑1 (MQ=255)                                                                         aCCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTg                                                                          <  2:493900/150‑1 (MQ=255)                                                                          ccGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTc                                                                        <  2:212042/151‑1 (MQ=255)                                                                                aaaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTc                                                                  <  1:291505/151‑1 (MQ=255)                                                                                aaaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTc                                                                  <  2:358311/151‑1 (MQ=255)                                                                                 aaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTcc                                                                 >  2:172719/1‑151 (MQ=255)                                                                                    aGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACg                                                              >  2:117948/1‑151 (MQ=255)                                                                                      gCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgtg                                                            <  1:451406/151‑1 (MQ=255)                                                                                       cATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgtgt                                                           >  2:9698/1‑151 (MQ=255)                                                                                         tCAGTGGCGGATGCGGGGAACGCCCATCCAGATATTCCATAATGATGCGAGATTCCCCCGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTTCACgtgt                                                           >  1:228080/1‑149 (MQ=255)                                                                                                  gATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGa                                                 <  1:401913/150‑1 (MQ=255)                                                                                                    tGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGaa                                                >  1:663451/1‑149 (MQ=255)                                                                                                       gggAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt                                            >  2:358096/1‑150 (MQ=255)                                                                                                           aaCGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACAcc                                       >  1:325239/1‑151 (MQ=255)                                                                                                             cGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCtt                                     >  2:31541/1‑151 (MQ=255)                                                                                                              ggtCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACCCCtct                                    <  1:290976/149‑3 (MQ=255)                                                                                                                 cATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtc                                   <  1:354501/149‑1 (MQ=255)                                                                                                                  aTCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTca                                <  2:561074/151‑1 (MQ=255)                                                                                                                        atatTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagccagc                          <  2:246739/151‑1 (MQ=255)                                                                                                                               aTAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCATGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGc                   >  1:253632/1‑151 (MQ=255)                                                                                                                                tAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCg                  >  1:82785/1‑151 (MQ=255)                                                                                                                                   tgatgCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGAc               >  2:597965/1‑151 (MQ=255)                                                                                                                                     atgCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCt             <  1:178198/151‑1 (MQ=255)                                                                                                                                      tgCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTg            >  2:530948/1‑151 (MQ=255)                                                                                                                                       gCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCATTTCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACCTACCCCTTTCTCAGCCAGCACAATGCGGACCtt            <  1:410758/150‑2 (MQ=255)                                                                                                                                         gagaTTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATg         >  1:393230/1‑151 (MQ=255)                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAg  <  1:598396/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAG  >  NC_000913/3377068‑3377361

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