Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,749,745	A→G	C202R (TGT→CGT)	ydhU ←	putative cytochrome YdhU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,749,745	0	A	G	97.4%	125.8 / ‑4.8	38	C202R (TGT→CGT)	ydhU	putative cytochrome YdhU
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/20);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TACGCCCCGTGGTGCAAAGATAAAGATGACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACACAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTCCCCCTGCATAATGCCAA  >  NC_000913/1749613‑1749876                                                                                                                                     |                                                                                                                                    tACGCCCCGTGGTGCAAAGATAAAGATGACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCagcagc                                                                                                                              <  1:415327/140‑1 (MQ=255)  aCGCCCCGTGGTGCAAAGATAAAGATGACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGcc                                                                                                                             <  2:696080/140‑1 (MQ=255)        cGTGGTGCAAAGATAAAGATGACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCa                                                                                                                       <  1:175491/140‑1 (MQ=255)                       aagatGACCGAAGATCAAAAAGTGGATTACAAATGCCTGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGgca                                                                                                        <  1:213274/140‑1 (MQ=255)                       aagatGACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGgca                                                                                                        >  1:375838/1‑140 (MQ=255)                            gACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACACAGCAGCCCCGTCAATAGTAACAACGGCAGCAAt                                                                                                   <  2:159776/140‑1 (MQ=255)                              ccGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGAGGCAGCAGGCCCTTCAATAGTATCAACGGCAGCAATcc                                                                                                 >  1:192483/1‑140 (MQ=255)                              ccGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATcc                                                                                                 >  1:668384/1‑140 (MQ=255)                              ccGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATcc                                                                                                 >  2:487397/1‑140 (MQ=255)                                    aTAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACAt                                                                                           >  1:538240/1‑140 (MQ=255)                                       aaaaaGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGac                                                                                        <  2:668384/140‑1 (MQ=255)                                       aaaaaGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGac                                                                                        <  1:590830/140‑1 (MQ=255)                                               gCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACAt                                                                                >  1:28046/1‑140 (MQ=255)                                                      aaTGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGAc                                                                         <  2:397716/140‑1 (MQ=255)                                                       aTGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGAcc                                                                        <  2:140902/140‑1 (MQ=255)                                                       aTGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGAcc                                                                        <  1:668342/140‑1 (MQ=255)                                                       aTGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGAcc                                                                        >  1:427107/1‑140 (MQ=255)                                                       aTGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGAcc                                                                        >  1:184871/1‑140 (MQ=255)                                                           cAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGTAATCCATACATGACACCAACATAGGCGACctgct                                                                    >  1:107008/1‑140 (MQ=255)                                                              agCAAAATGTGTCTGCAATAACCAGTATCTTCCGCAAGGTAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTa                                                                 <  2:192483/140‑1 (MQ=255)                                                                        gtCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTaaa                                                       >  1:72250/1‑140 (MQ=255)                                                                         tCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAAt                                                      <  2:72250/140‑1 (MQ=255)                                                                               aTAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGAc                                                <  2:184871/140‑1 (MQ=255)                                                                               aTAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGAc                                                <  2:375838/140‑1 (MQ=255)                                                                                tAACCAGTATCTTACGCCAGGAAACACATTTCCCACGCCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACt                                               <  1:312140/140‑1 (MQ=255)                                                                                   ccAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTggg                                            >  1:91630/1‑140 (MQ=255)                                                                                        aTCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGgttgtt                                       <  2:49183/140‑1 (MQ=255)                                                                                              cGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCggg                                 >  2:413998/1‑140 (MQ=255)                                                                                              cGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCggg                                 >  1:243444/1‑140 (MQ=255)                                                                                               gCCAGGAAACACATCTCCCACGACTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGa                                <  2:243444/140‑1 (MQ=255)                                                                                                 cAGGAAACACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATATGCGACCTGATGTAAGGGATTAAATTGAGACAGGGTTGTTGCCGGGaaa                              >  2:202624/1‑140 (MQ=255)                                                                                                     aaaCACATCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGAt                          <  2:538240/140‑1 (MQ=255)                                                                                                          caTCTCCCACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGttct                     <  1:413998/140‑1 (MQ=255)                                                                                                               cccACGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTcccc                <  2:427107/140‑1 (MQ=255)                                                                                                                   cGGCTTGCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTCCCCCTGc            >  2:33858/1‑140 (MQ=255)                                                                                                                         gCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTCCCCCTGCATAATg      >  1:612342/1‑140 (MQ=255)                                                                                                                         gCGGATAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTCCCCCTGCATAATg      >  1:526337/1‑140 (MQ=255)                                                                                                                             aTAGAGACGCAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTCCCCCTGCATAATGCcaa  <  1:16245/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    TACGCCCCGTGGTGCAAAGATAAAGATGACCGAAGATAAAAAAGAGGCTTATAAATGCCAGAGCAAAATGTGTCTGCAATAACCAGTATCTTACGCCAGGAAACACATCTCCCACGGCTTGCGGATAGAGACACAGCAGCCCCGTCAATAGTAACAACGGCAGCAATCCATACATGACACCAACATAGGCGACCTGCTGTAAGGGATTAAATTTAGACTGGGTTGTTGCCGGGAAAGGATGTTCTTCCCCCTGCATAATGCCAA  >  NC_000913/1749613‑1749876

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