Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,189,827	G→A	6.8%	P14P (CCG→CCA)	yjaZ →	conserved protein YjaZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,189,827	0	G	A	6.8%	159.9 / 2.3	59	P14P (CCG→CCA)	yjaZ	conserved protein YjaZ
Reads supporting (aligned to +/- strand):  ref base G (18/37);  new base A (2/2);  total (20/39)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACAAAAGATTATGAAC  >  NC_000913/4189698‑4189963                                                                                                                                  |                                                                                                                                         aCTAATATCGTGAAACTGACTGAACTAATTGAGTCAAACGCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTc                                                                                                                             <  2:138380/143‑1 (MQ=255) aCTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTc                                                                                                                             >  2:22330/1‑143 (MQ=255)   tAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGt                                                                                                                           <  1:806356/143‑1 (MQ=255)     atatCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGa                                                                                                                         <  1:886937/143‑1 (MQ=255)     atatCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGa                                                                                                                         <  1:874899/143‑1 (MQ=255)          cTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCAGATGGTGATATTCGTGATTTAt                                                                                                                    >  1:699863/1‑143 (MQ=255)          cTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCAGATGGTGATATTCGTGATTTAt                                                                                                                    >  1:200757/1‑143 (MQ=255)               actgactgAACTTATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCAttt                                                                                                               <  2:274170/143‑1 (MQ=255)                    ctggACTAATGGAGTCAATCTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGa                                                                                                          <  2:869027/143‑1 (MQ=255)                    ctgaACTATTTGTGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGa                                                                                                          <  2:519967/143‑1 (MQ=255)                      gaACTACTTGATGCAACTTCGGCAATGACTACGTACTATTCCTGTGTAACTTTCTTAAGGAACGAGAACGAAACTGGACGTGGAAAAGTAGCGACCTTTTGGACATCCGGATGGTGATATTCGTCATTTATCATTTCTTGATg                                                                                                        <  2:184904/143‑1 (MQ=255)                          tAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCtca                                                                                                    >  1:987302/1‑143 (MQ=255)                          tAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCtca                                                                                                    >  2:630572/1‑143 (MQ=255)                          caaTGGAGTCAAACTCGGAAAGGATTCGATACTATTCCTGTGTAACTTTCTTAATGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCtca                                                                                                    <  2:116021/142‑1 (MQ=255)                                  tCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAGAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGc                                                                                                       >  2:78989/1‑132 (MQ=255)                                  tCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAGAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGc                                                                                                       <  1:78989/132‑1 (MQ=255)                                   ccaaCTCGGAAGGGATTCGGTACTATTCCTGTGTAACTTTCTTAAGGGACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGt                                                                                           <  2:85998/141‑1 (MQ=255)                                       cTCGGCAAGGATTCGATTCTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTAc                                                                                       <  1:130893/143‑1 (MQ=255)                                          ggCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGtt                                                                                    <  1:997117/143‑1 (MQ=255)                                          ggCAAGGATTCGATACGATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTATTCGtt                                                                                    <  2:676630/143‑1 (MQ=255)                                            caagtgttCGATAGTAGTCCTGTGTGACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCa                                                                                  <  2:621822/137‑1 (MQ=255)                                             aaGGATTCGATACTATTCCTGTGTAACTTTCTTCAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGATATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAg                                                                                 <  1:952946/143‑1 (MQ=255)                                             aaGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAg                                                                                 <  2:45748/143‑1 (MQ=255)                                             aaGGATTCGATACTAGTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAg                                                                                 <  1:941478/143‑1 (MQ=255)                                                gATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGcat                                                                              >  1:612245/1‑143 (MQ=255)                                                gATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGAGATTTTTTTGAGTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGcat                                                                              >  1:1005882/1‑143 (MQ=255)                                                   tCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGCTATTCGTGATTTATCATTTCTTGCTGCTCATCCGGCTGCCTACGTTCAGcatcat                                                                           >  1:259742/1‑143 (MQ=255)                                                   tCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGcatcat                                                                           >  1:776861/1‑143 (MQ=255)                                                   tCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGcatcat                                                                           <  2:91158/143‑1 (MQ=255)                                                        aCTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACAGCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCGTCAGGAggg                                                                      >  1:724192/1‑143 (MQ=255)                                                            ttCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCaaa                                                                  <  1:326245/143‑1 (MQ=255)                                                                   ttgtgtttCTTAGGGATGGAGGATTGAACGGGAAGTGGAAAAGTGGGGACCTTTTGGACATCCGGATGGTGATCTTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGTATCATGAGGGCAAAGAGCCtt                                                           <  2:888560/138‑1 (MQ=255)                                                                    tAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTGGTGATTTATCATTGCTTGATGCGCATGAGGCTGTCTACGTTCAGGATCGTGAGGGCGACGAGACttt                                                          >  1:480954/1‑143 (MQ=255)                                                                    tAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCttt                                                          >  1:815332/1‑143 (MQ=255)                                                                       cTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCTTGAGGGCAAAGAGACTTTagc                                                       >  2:382779/1‑142 (MQ=255)                                                                        tttCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTagag                                                      >  1:201737/1‑143 (MQ=255)                                                                            ttcaGGAACGAGGATGAAACAGGAAGTGGAAGAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATc                                                  <  2:878230/140‑1 (MQ=255)                                                                            ttAAGGAACGGGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTTGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATc                                                  <  2:191860/143‑1 (MQ=255)                                                                                    cGAGGATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTgg                                          <  1:926400/143‑1 (MQ=255)                                                                                    cGAGAATGGAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTgg                                          <  2:945699/143‑1 (MQ=255)                                                                                     gagaATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACTTCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTggg                                         <  1:480881/143‑1 (MQ=255)                                                                                         aTGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTGGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTAc                                     <  2:32055/143‑1 (MQ=255)                                                                                                     gTGGAAAAGTGGCTGACTGTGGGACATCCGGATGGTGATATTCGTGATTTATCATTTCGTGATGCTCATCAGGCTGTCTACGTTCCGCATAATGTTGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCa                         <  2:217288/143‑1 (MQ=255)                                                                                                        gAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTg                      >  1:822020/1‑143 (MQ=255)                                                                                                          aaaGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCt                    <  1:22330/143‑1 (MQ=255)                                                                                                           aaGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCtt                   <  1:269599/143‑1 (MQ=255)                                                                                                             gTGGCGACCTTTTGGTTATCCGGATGGTGGTATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGCCTGTGTTCAGCATCATGAGGGCAAAGAGCCTTGAGAGTAGCGCTTTGGGGTAACCTACTCTCGTCACTGCTTca                 <  2:1006003/143‑1 (MQ=255)                                                                                                                gCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaa              >  1:864221/1‑143 (MQ=255)                                                                                                                gCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGCGTTACCTACTCTCTTCACTGCTCCACaa              >  1:540741/1‑143 (MQ=255)                                                                                                                 cGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaaa             >  1:16628/1‑143 (MQ=255)                                                                                                                  gtccTTTTGGACATCCAGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaaaa            <  2:200757/141‑1 (MQ=255)                                                                                                                  gatcTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATCCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGCGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaaaa            <  2:251173/140‑1 (MQ=255)                                                                                                                  gACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaaaa            >  1:489427/1‑143 (MQ=255)                                                                                                                  gACCTTTTGGACATCCAGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaaaa            <  2:699863/143‑1 (MQ=255)                                                                                                                  gACCTTGTGGACATCCGGATTGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACaaaa            <  1:240806/143‑1 (MQ=255)                                                                                                                    gcTTTTGGACATTCGGATGGTGATATTCGTGATTGATCAGTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACAAAAGa          <  1:408157/142‑1 (MQ=255)                                                                                                                    ccTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACAAAAGa          <  1:624594/143‑1 (MQ=255)                                                                                                                        ttGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACAAAAGATTAt      <  1:499776/143‑1 (MQ=255)                                                                                                                            aCATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACAAAAGATTATGAAc  <  1:267157/143‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         ACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATACTATTCCTGTGTAACTTTCTTAAGGAACGAGAATGAAACAGGAAGTGGAAAAGTGGCGACCTTTTGGACATCCGGATGGTGATATTCGTGATTTATCATTTCTTGATGCTCATCAGGCTGTCTACGTTCAGCATCATGAGGGCAAAGAGCCTTTAGAGTATCGCTTTTGGGTTACCTACTCTCTTCACTGCTTCACAAAAGATTATGAAC  >  NC_000913/4189698‑4189963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑