Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,066,319	Δ1 bp	10.0%	coding (855/897 nt)	ygfI ←	putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,066,319	0	C	.	10.0%	165.2 / 12.3	50	coding (855/897 nt)	ygfI	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base . (2/3);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGC  >  NC_000913/3066179‑3066464                                                                                                                                               |                                                                                                                                                  cGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTga                                                                                                                                            <  3:64590/148‑2 (MQ=38)    tATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACggg                                                                                                                                        <  1:44511/149‑2 (MQ=37)     aTCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGt                                                                                                                                        >  4:136080/1‑148 (MQ=255)         ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATaa                                                                                                                                    >  4:82788/1‑148 (MQ=255)          gAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                                                                   >  4:41232/1‑148 (MQ=255)               tGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatgat                                                                                                                             <  2:40044/149‑1 (MQ=255)                   ccTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACAAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                                          >  3:43635/1‑148 (MQ=255)                     tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                                                                        <  3:379739/148‑1 (MQ=255)                     tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                                                                       <  1:86725/149‑1 (MQ=255)                          gACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                                                                   >  4:228317/1‑148 (MQ=255)                                     cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATAtt                                                                                                        <  3:82788/148‑1 (MQ=255)                                      acGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                                       >  4:32807/1‑148 (MQ=255)                                        gATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                                                                    <  1:113978/149‑1 (MQ=255)                                          tCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                                                  >  2:279509/1‑149 (MQ=255)                                             tAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                                                                >  3:49307/1‑148 (MQ=255)                                                   cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                                          <  3:320548/148‑1 (MQ=255)                                                             ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                                                                >  3:393954/1‑148 (MQ=255)                                                              tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                                                               >  4:344853/1‑148 (MQ=255)                                                                  ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                                                           <  4:393954/148‑1 (MQ=255)                                                                  ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                                                                          >  1:314259/1‑149 (MQ=255)                                                                         ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                                                             <  4:362930/139‑1 (MQ=255)                                                                         ttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                                                             >  3:362930/1‑139 (MQ=255)                                                                             ttGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGt                                                                                                                                              >  4:354143/1‑70 (MQ=255)                                                                             ttGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGt                                                                                                                                              <  3:354143/70‑1 (MQ=255)                                                                              tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                                                              <  1:36566/149‑1 (MQ=255)                                                                                cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg                                                             <  3:254566/148‑1 (MQ=255)                                                                                   ccGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGt                                                          <  3:122663/148‑1 (MQ=255)                                                                                      gCTCGATGCTGGGTTAATGACTCGTCCCGTCACAAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGc                                                       <  4:43635/148‑1 (MQ=255)                                                                                         cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCTTATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg                                                    >  3:37847/1‑148 (MQ=255)                                                                                          gATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt                                                  <  1:303814/149‑1 (MQ=255)                                                                                            tGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTc                                                <  2:314259/149‑1 (MQ=255)                                                                                              cTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                                              <  1:66892/149‑1 (MQ=255)                                                                                                gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCa                                            >  1:360226/1‑149 (MQ=255)                                                                                                gggTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCAGGATTCAcc                                             >  4:52391/1‑148 (MQ=255)                                                                                                   ttAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTc                                         <  2:324414/149‑1 (MQ=255)                                                                                                       tGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGcc                                      <  3:32807/148‑1 (MQ=255)                                                                                                         aCTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg                                                        >  3:400201/1‑128 (MQ=255)                                                                                                         aCTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg                                                        <  4:400201/128‑1 (MQ=255)                                                                                                                ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATg                            >  2:304955/1‑149 (MQ=255)                                                                                                                  gTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                          >  1:395076/1‑149 (MQ=255)                                                                                                                   tCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                          <  3:344853/148‑1 (MQ=255)                                                                                                                   tCACGAT‑CTTCAGCAATA‑TGCTGAAGCAACGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATCCgg                         <  1:47902/149‑1 (MQ=255)                                                                                                                    cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                         >  4:174499/1‑148 (MQ=255)                                                                                                                        atgcTTCAGCATAT‑TGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAata                                                                   >  4:10276/4‑102 (MQ=38)                                                                                                                        atgcTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAata                                                                   <  3:10276/99‑1 (MQ=38)                                                                                                                        atgcTTCAGC‑ATATTGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATc                     >  4:253471/4‑148 (MQ=255)                                                                                                                                    ta‑tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGaaaa         >  1:38890/1‑149 (MQ=255)                                                                                                                                           gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGTTCACCGTCAGAAAATGTcc    >  4:97004/1‑148 (MQ=255)                                                                                                                                           gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCcg   >  2:250293/1‑149 (MQ=255)                                                                                                                                             aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCcgc  >  3:190384/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                  CGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGC  >  NC_000913/3066179‑3066464

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