Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,629,481	G→A	*64* (TAG→TAA)	yfgG →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,629,481	0	G	A	98.4%	195.7 / ‑6.4	64	*64* (TAG→TAA)	yfgG	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (27/35);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTT  >  NC_000913/2629348‑2629608                                                                                                                                      |                                                                                                                                tCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGtt                                                                                                                            <  2:292943/139‑1 (MQ=255)  cGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTa                                                                                                                           <  2:340811/139‑1 (MQ=255)    tACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTAcc                                                                                                                         <  2:925188/139‑1 (MQ=255)           catcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCa                                                                                                                  >  2:1598555/1‑139 (MQ=255)              cacatttAGCTGCGTCTTTAGCCGTTTTATCTACTCGTCCGTGGGTGCCTCGCAGCACCACCTTAACAAAAGAGAAGCTCAGCAATCCACACTCTCCGTCCAATCACCGTTACAACGTTAACGGTTACCTTCTCCACtt                                                                                                               <  2:713102/135‑1 (MQ=255)              caGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACtt                                                                                                               >  1:1336671/1‑139 (MQ=255)                gCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTc                                                                                                             <  1:1493231/139‑1 (MQ=255)                   ttATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTcaca                                                                                                          <  1:926584/139‑1 (MQ=255)                      tcttcttctTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGaa                                                                                                       <  1:1043524/139‑1 (MQ=255)                                                tCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGc                                                                             >  2:893413/1‑139 (MQ=255)                                                tCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGc                                                                             >  1:959664/1‑139 (MQ=255)                                                   tccgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCttt                                                                          >  1:489746/1‑139 (MQ=255)                                                     cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                        <  2:717703/139‑1 (MQ=255)                                                     cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                        >  2:1181946/1‑139 (MQ=255)                                                     cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                        <  1:223928/139‑1 (MQ=255)                                                             gcTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCGCACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGTATGCTTTTGCTGTTTgg                                                                >  1:1017113/2‑139 (MQ=255)                                                                     caCCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCa                                                        >  1:1563031/1‑139 (MQ=255)                                                                     caCCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCa                                                        >  2:547658/1‑139 (MQ=255)                                                                     caCCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACAGTTACAACGTTAACGGTTACCTTCTGCACTTTCCCAGCACATCACGCCACTTCGCTGTCGGATGCTTTTGCTGTTTGCGATTATCa                                                        >  2:1013993/1‑139 (MQ=255)                                                                          tCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCg                                                                                  <  1:300348/108‑1 (MQ=255)                                                                          tCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCg                                                                                  >  2:300348/1‑108 (MQ=255)                                                                             gagCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCa                                                <  2:1336671/139‑1 (MQ=255)                                                                                cAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGat                                             >  1:417407/1‑139 (MQ=255)                                                                                 aaaaaaGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTcaca                                                                                                          >  1:206592/1‑77 (MQ=255)                                                                                 aaaaaaGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTcaca                                                                                                          <  2:206592/77‑1 (MQ=255)                                                                                   aaaaGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGttt                                                                  >  2:1236748/1‑115 (MQ=255)                                                                                   aaaaGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGttt                                                                  <  1:1236748/115‑1 (MQ=255)                                                                                       gaagCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGTAGATATTCgtt                                      <  2:959664/139‑1 (MQ=255)                                                                                          gCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                            >  1:1543266/1‑66 (MQ=255)                                                                                          gCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                            <  2:1543266/66‑1 (MQ=255)                                                                                               gCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTaa                              >  2:697261/1‑139 (MQ=255)                                                                                                  aTCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGc                                                    >  1:740024/1‑114 (MQ=255)                                                                                                  aTCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGc                                                    <  2:740024/114‑1 (MQ=255)                                                                                                     cacaCTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTa                        <  2:417407/139‑1 (MQ=255)                                                                                                     cacaCTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTa                        <  2:1017113/139‑1 (MQ=255)                                                                                                     cacaCTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTa                        <  2:1563031/139‑1 (MQ=255)                                                                                                        acTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGa                                              <  1:918925/114‑1 (MQ=255)                                                                                                        acTCTCCGTCGAATCACCAGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGa                                              >  2:918925/1‑114 (MQ=255)                                                                                                          tctcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCt                   <  1:24252/139‑1 (MQ=255)                                                                                                          tctcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCt                   >  2:24252/1‑139 (MQ=255)                                                                                                            tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACt                                                                                           >  2:542356/1‑65 (MQ=255)                                                                                                            tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACt                                                                                           <  1:542356/65‑1 (MQ=255)                                                                                                            tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAAC‑GCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTa                >  1:783629/1‑139 (MQ=255)                                                                                                              tgTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAg                                                                                                         >  2:1034747/2‑49 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTATCTTCTCCACTTTCACAg                                                                                                         <  1:1034747/49‑1 (MQ=38)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                            <  2:1646422/46‑1 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                            >  1:1646422/1‑46 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCa                                                                                             >  2:1463583/1‑61 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCa                                                                                             <  1:1463583/61‑1 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                        <  2:908767/82‑1 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                        >  1:908767/1‑82 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCt                   >  2:1357972/1‑135 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCt                   <  1:1357972/135‑1 (MQ=255)                                                                                                              cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTAt               >  2:435867/1‑139 (MQ=255)                                                                                                                    aTCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAAC‑GCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact        <  2:783629/139‑1 (MQ=255)                                                                                                                     tCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact        <  1:435867/139‑1 (MQ=255)                                                                                                                     tCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact        <  1:1181946/139‑1 (MQ=255)                                                                                                                     tCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact        <  1:547658/139‑1 (MQ=255)                                                                                                                     tCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact        <  1:563483/139‑1 (MQ=255)                                                                                                                      cACCGGTACAACGTTAACGGTTACCTTCTCCACtt                                                                                                               <  2:1091397/35‑1 (MQ=255)                                                                                                                      cACCGGTACAACGTTAACAGTTACCTTCTCCACtt                                                                                                               >  1:1091397/1‑35 (MQ=255)                                                                                                                        ccGGTACAACGTTATCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGCATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGAt     <  1:1013993/139‑1 (MQ=255)                                                                                                                           gTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGtt  <  2:489746/139‑1 (MQ=255)                                                                                                                           gTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGtt  <  1:893413/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                TCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTT  >  NC_000913/2629348‑2629608

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