Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,205,731	G→A	*61* (TAG→TAA)  pseudogene (1/412 nt)	ymfR → beeE →	e14 prophage, uncharacterized protein pseudogene, portal protein family, e14 prophage,Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,205,731	0	G	A	97.9%	144.5 / ‑4.3	47	*61* (TAG→TAA)  pseudogene (1/412 nt)	ymfR beeE	e14 prophage, uncharacterized protein pseudogene, portal protein family, e14 prophage,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/24);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGG  >  NC_000913/1205596‑1205854                                                                                                                                        |                                                                                                                            tttGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCACTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAatgt                                                                                                                          <  1:1422173/139‑1 (MQ=255)  ttGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGtt                                                                                                                         <  2:997237/139‑1 (MQ=255)  ttGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGtt                                                                                                                         <  2:1544216/139‑1 (MQ=255)        acGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTc                                                                                                                   >  1:1173850/2‑139 (MQ=255)          gTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCgg                                                                                                                 >  1:652029/1‑139 (MQ=255)           tATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCggg                                                                                                                <  2:972839/139‑1 (MQ=255)                  ccTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTAttt                                                                                                         >  1:880901/1‑139 (MQ=255)                                    cggcggGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTAttt                                                                                                         >  2:1132769/1‑121 (MQ=255)                                    cggcggGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTAttt                                                                                                         <  1:1132769/121‑1 (MQ=255)                                      gcggGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCAcc                                                                                     >  1:1339540/1‑139 (MQ=255)                                       cggGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCg                                                                                    <  2:1578715/139‑1 (MQ=255)                                            tttGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGAc                                                                               <  1:1103075/139‑1 (MQ=255)                                             ttgttgttgCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGAcc                                                                              >  1:1690067/1‑139 (MQ=255)                                                        ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCcag                                                                       <  2:1855141/135‑1 (MQ=255)                                                        ggggCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCcag                                                                       >  1:1855141/1‑135 (MQ=255)                                                            cgcTGTTCGTGTTCTGGTCGTGGCTAGTGGGGCGATATCTCCACCGTACACAGGCGTCTGTCCGCGGAGGTAAATAATGTTCTTGTCGGGATTACTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGCGCTg                                                               <  1:1875089/139‑1 (MQ=255)                                                                 tgCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGa                                                          <  2:1339540/139‑1 (MQ=255)                                                                           ggTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                <  2:1690067/139‑1 (MQ=255)                                                                            gTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGc                                               >  1:34304/1‑139 (MQ=255)                                                                            gTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGc                                               <  2:1173850/139‑1 (MQ=255)                                                                                 ggCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGa                                                          >  2:609245/1‑123 (MQ=255)                                                                                 ggCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGa                                                          <  1:609245/123‑1 (MQ=255)                                                                                 ggCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                <  2:606382/133‑1 (MQ=255)                                                                                 ggCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                >  1:606382/1‑133 (MQ=255)                                                                                 ggCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCg                                          >  1:887557/1‑139 (MQ=255)                                                                                  gCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGt                                         <  2:1318925/139‑1 (MQ=255)                                                                                    tggtggCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTAt                                       <  2:34304/139‑1 (MQ=255)                                                                                       tggCGTGATATCTCGACCGTCCACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGac                                    <  2:1240985/139‑1 (MQ=255)                                                                                            cgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                <  1:1396333/122‑1 (MQ=255)                                                                                            cgATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                >  2:1396333/1‑122 (MQ=255)                                                                                             gATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCAc                                                                            <  1:875530/93‑1 (MQ=255)                                                                                             gATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCAc                                                                            >  2:875530/1‑93 (MQ=255)                                                                                                      aCCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACg                                                                                         <  1:1351633/71‑1 (MQ=255)                                                                                                      aCCGTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACg                                                                                         >  2:1351633/1‑71 (MQ=255)                                                                                                         gTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAg                  >  1:884363/1‑139 (MQ=255)                                                                                                         gTACACAGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAg                  >  2:770262/1‑139 (MQ=255)                                                                                                               aGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                <  2:731624/103‑1 (MQ=255)                                                                                                               aGTCGTCTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCggg                                                >  1:731624/1‑103 (MQ=255)                                                                                                                 tcgtcTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGacac                                  >  2:1184263/1‑115 (MQ=255)                                                                                                                 tcgtcTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATcagcag          >  2:1428262/1‑139 (MQ=255)                                                                                                                 tcgtcTGTCGGCGGAGCTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGacac                                  <  1:1184263/115‑1 (MQ=255)                                                                                                                    tcTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGCCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCa       >  1:296630/1‑139 (MQ=255)                                                                                                                     cTGTCGGCGGAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAg      <  1:1428262/139‑1 (MQ=255)                                                                                                                         cggcggAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCggg  >  2:1310638/1‑139 (MQ=255)                                                                                                                           gcggAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGacac                                  >  2:481316/1‑105 (MQ=255)                                                                                                                           gcggAGGTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGacac                                  <  1:481316/105‑1 (MQ=255)                                                                                                                                ggTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCg                                          <  1:130950/92‑1 (MQ=255)                                                                                                                                ggTAAATAATGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCg                                          >  2:130950/1‑92 (MQ=255)                                                                                                                                        |                                                                                                                            TTTGCTGGCGTATGGTGCCTGGCTGATTTATCCCCCGGCGGGTTTTGTTGTTGCCGGGGCGCTGTGCCTGTTCTGGTCGTGGCTGGTGGCGCGATATCTCGACCGTACACAGTCGTCTGTCGGCGGAGGTAAATAGTGTTCTTTTCGGGATTATTTCAACGAAAAAGTGACGCACCGGTGACCACGCCAGCAGAGCTGGCGGATGCCATCGGGCTGTCGTATGACACCTATACCGGAAAGCAGATCAGCAGTCAGCGGG  >  NC_000913/1205596‑1205854

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