Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,431,682	T→G	E94A (GAA→GCA)	yrdD ←	ssDNA‑binding protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,431,682	0	T	G	97.8%	162.2 / ‑3.7	45	E94A (GAA→GCA)	yrdD	ssDNA‑binding protein, function unknown
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (25/19);  total (25/20)
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.75e-01

GTTTGAAGTTAATGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGAC  >  NC_000913/3431539‑3431829                                                                                                                                                |                                                                                                                                                    gTTTGAAGTTAATGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATg                                                                                                                                              >  1:697285/1‑151 (MQ=255)   ttGAAGTTAATGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGtt                                                                                                                                            >  2:875364/1‑151 (MQ=255)    tGAAGTTAATGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTc                                                                                                                                           >  1:251362/1‑151 (MQ=255)          tAATGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTc                                                                                                                                     >  1:449996/1‑151 (MQ=255)           aaTGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTc                                                                                                                                     >  2:718832/1‑150 (MQ=255)              ggCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGTCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAggg                                                                                                                                 <  2:326658/151‑1 (MQ=255)                cAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGAAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTa                                                                                                                               <  2:69174/151‑1 (MQ=255)                 aaaTTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGt                                                                                                                                >  1:482309/1‑149 (MQ=255)                  aaTTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGt                                                                                                                             >  2:656957/1‑151 (MQ=255)                   aTTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGtt                                                                                                                            >  2:38402/1‑151 (MQ=255)                      gACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTaa                                                                                                                          <  1:642916/150‑1 (MQ=255)                       acacTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATg                                                                                                                        <  1:457456/151‑1 (MQ=255)                        cacTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGc                                                                                                                       >  2:808159/1‑151 (MQ=255)                          cTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAAt                                                                                                                         <  2:11708/147‑1 (MQ=255)                             cGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCa                                                                                                                  <  1:169621/151‑1 (MQ=255)                             cGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCa                                                                                                                  <  2:546170/151‑1 (MQ=255)                             cGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCa                                                                                                                  <  2:197768/151‑1 (MQ=255)                                  aGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATaa                                                                                                              >  2:49530/1‑150 (MQ=255)                                      aTCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACAt                                                                                                          >  2:683452/1‑150 (MQ=255)                                       tCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATAc                                                                                                        >  2:178459/1‑151 (MQ=255)                                         acaAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATAcc                                                                                                       >  2:558452/1‑150 (MQ=255)                                            agagTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCaa                                                                                                     >  2:827191/1‑149 (MQ=255)                                              agTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCg                                                                                                 <  2:18913/151‑1 (MQ=255)                                                      gTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGc                                                                                         >  1:75515/1‑151 (MQ=255)                                                           gCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTaa                                                                                     <  2:541897/150‑1 (MQ=255)                                                                  cGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATAcc                                                                                 >  1:486363/1‑147 (MQ=255)                                                                  cGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGAt                                                                             <  2:243677/151‑1 (MQ=255)                                                                            gcTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCcac                                                                    >  2:333178/1‑150 (MQ=255)                                                                               ggACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGc                                                                <  2:807182/151‑1 (MQ=255)                                                                                   cAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCGGATTTGCGCCACATGCAggg                                                            <  2:371322/151‑1 (MQ=255)                                                                                    aGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGc                                                           <  2:163429/151‑1 (MQ=255)                                                                                             gTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGcc                                                  >  2:123182/1‑151 (MQ=255)                                                                                                cGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCt                                               >  2:482997/1‑151 (MQ=255)                                                                                                cGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCt                                               >  2:631010/1‑151 (MQ=255)                                                                                                         gggCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACtt                                      >  2:836159/1‑151 (MQ=255)                                                                                                          ggCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACttt                                     <  1:851493/151‑1 (MQ=255)                                                                                                               gTAATTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACg                                 >  1:845393/1‑150 (MQ=255)                                                                                                                 aaTTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGata                              <  2:469710/151‑1 (MQ=255)                                                                                                                 aaTTGCTGTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGat                               >  2:181860/1‑150 (MQ=255)                                                                                                                        gTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCa                        <  1:339486/150‑1 (MQ=255)                                                                                                                        gTTTCGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCAt                       <  2:350138/151‑1 (MQ=255)                                                                                                                            cGTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCg                    >  2:379845/1‑150 (MQ=255)                                                                                                                             gTCCGGTTTATCGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCt                  >  2:248482/1‑151 (MQ=255)                                                                                                                                        cGATAAGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCag       <  1:373338/151‑1 (MQ=255)                                                                                                                                             aGTGCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAggac  <  1:822423/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    GTTTGAAGTTAATGGCAAATTGACACTCCGGGTAGCGATCACAAGAGTGAAATGTTTTGCCATAACGGGAGCGGCGCTGGACCAGATGGCCCGTCCGACATTGGGGGCATGTAATTGCTGTTTCGTCCGGTTTATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGAC  >  NC_000913/3431539‑3431829

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