Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,020,512	A→C	13.8%	intergenic (+23/‑28)	ssnA → / → ygfM	putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,020,512	0	A	C	13.8%	142.7 / 19.4	53	intergenic (+23/‑28)	ssnA/ygfM	putative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base C (3/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

AGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACC  >  NC_000913/3020364‑3020652                                                                                                                                                      |                                                                                                                                             aGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCa                                                                                                                                              <  2:101127/149‑1 (MQ=255) aGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCa                                                                                                                                              <  2:75790/149‑1 (MQ=255)         ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAagag                                                                                                                                      <  2:313214/149‑1 (MQ=255)         ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                      <  2:9468/149‑1 (MQ=255)          gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAgg                                                                                                                                     <  3:174305/149‑1 (MQ=255)            cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGt                                                                                                                                   <  4:100259/149‑1 (MQ=255)               tGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATACCCTCTTCCC‑TGCGGGAAGAgggcatc                                                                                                                                <  3:132386/149‑5 (MQ=255)                   gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATaa                                                                                                                            <  3:602251/149‑1 (MQ=255)                   gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATaa                                                                                                                            <  3:305045/149‑1 (MQ=255)                   gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATaa                                                                                                                            >  2:177130/1‑149 (MQ=255)                            tATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                   >  2:307920/1‑149 (MQ=255)                                      gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAAc                                                                                                         >  4:317017/1‑149 (MQ=255)                                         aGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAt                                                                                                      <  3:187726/149‑1 (MQ=255)                                               aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                              <  4:99271/135‑1 (MQ=255)                                               aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                              >  3:99271/1‑135 (MQ=255)                                                 ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTaaa                                                                                                                               >  4:595346/1‑116 (MQ=255)                                                 ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTaaa                                                                                                                               <  3:595346/116‑1 (MQ=255)                                                          gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGC‑GGGAAGAgg                                                                                                                                     >  1:560150/1‑101 (MQ=255)                                                          gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                     <  2:560150/101‑1 (MQ=255)                                                                 tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                              >  3:206896/1‑149 (MQ=255)                                                                      gCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGgcg                                                                         >  3:293920/1‑149 (MQ=255)                                                                        acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                       >  2:275513/1‑149 (MQ=255)                                                                          aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                      >  3:130751/1‑84 (MQ=255)                                                                          aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                      <  4:130751/84‑1 (MQ=255)                                                                            gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                   >  2:45876/1‑149 (MQ=255)                                                                             ccAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAAc                                                                  >  4:440011/1‑149 (MQ=255)                                                                                 aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                              <  4:293920/149‑1 (MQ=255)                                                                                              aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAg                                                 <  1:89608/149‑1 (MQ=255)                                                                                                tATGTGGCGTCGAATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGa                                               <  3:121809/149‑1 (MQ=255)                                                                                                    tgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                           <  1:177130/149‑1 (MQ=255)                                                                                                        gTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCg                                       >  4:137388/1‑149 (MQ=255)                                                                                                           ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGC‑GGGAagag                                                                                                                                      <  1:266264/51‑1 (MQ=25)                                                                                                           ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAagag                                                                                                                                      >  2:266264/1‑51 (MQ=25)                                                                                                           ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAggg                                                                                                                                           <  4:51424/46‑1 (MQ=255)                                                                                                           ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAggg                                                                                                                                           >  3:51424/1‑46 (MQ=255)                                                                                                                gatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTc                               <  3:440011/149‑1 (MQ=255)                                                                                                                       tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                        >  2:416370/1‑149 (MQ=255)                                                                                                                       tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                        >  2:21918/1‑149 (MQ=255)                                                                                                                          cATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgca                     <  1:416370/149‑1 (MQ=255)                                                                                                                           aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                    >  4:162852/1‑149 (MQ=255)                                                                                                                               aTGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCaaa                >  2:134034/1‑149 (MQ=255)                                                                                                                                    aGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa           >  4:348758/1‑149 (MQ=255)                                                                                                                                     gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt          <  3:317017/149‑1 (MQ=255)                                                                                                                                          ccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                           <  4:574961/111‑1 (MQ=255)                                                                                                                                          ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                           >  3:574961/1‑111 (MQ=255)                                                                                                                                          ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                           >  3:108234/1‑111 (MQ=255)                                                                                                                                          ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                           <  4:108234/111‑1 (MQ=255)                                                                                                                                          ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                             <  1:405037/77‑1 (MQ=255)                                                                                                                                          ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                             >  2:405037/1‑77 (MQ=255)                                                                                                                                          ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac     >  3:101501/1‑149 (MQ=255)                                                                                                                                          ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGACGGGGGCAGCAAACTCAATGCTac     >  2:406927/1‑149 (MQ=255)                                                                                                                                           ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca    >  1:72419/1‑149 (MQ=255)                                                                                                                                           ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca    >  3:603491/1‑149 (MQ=255)                                                                                                                                           ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCC‑GGGGCAGCAAACTCAATGCTacac   <  3:429885/149‑1 (MQ=255)                                                                                                                                             ctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACAcc  <  1:275513/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             AGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACC  >  NC_000913/3020364‑3020652

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: ‑