Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,734,543	(C)8→9	coding (492/732 nt)	yfiH ←	UPF0124 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,734,535	1	.	C	95.2%	149.9 / ‑0.6	42	coding (500/732 nt)	yfiH	UPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (22/18);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734432‑2734663                                                                                                          |                                                                                                                                   aCGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:1175971/139‑1 (MQ=255)                          tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                                                       >  1:1389810/1‑139 (MQ=255)                             gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                                                    >  2:424999/1‑139 (MQ=255)                                aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                                                 <  2:643352/139‑1 (MQ=255)                                       tATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                          <  2:382482/139‑1 (MQ=255)                                       tATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                          <  1:863887/139‑1 (MQ=255)                                       tATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                          <  1:6755/139‑1 (MQ=255)                                        aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGc                                                                                                >  2:571119/1‑100 (MQ=255)                                        aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:571119/100‑1 (MQ=255)                                        aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCCAATTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAAAATATTTTCCGGATTATCAga                                                         >  1:1254814/1‑138 (MQ=255)                                             cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAAc                                                    >  1:764039/1‑139 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGaa                                                                                                                  >  2:719154/1‑75 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGaa                                                                                                                  <  1:719154/75‑1 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                      >  1:791461/1‑103 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                      <  2:791461/103‑1 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                  <  2:510619/107‑1 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                  >  1:510619/1‑107 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                                                  >  1:441054/1‑139 (MQ=255)                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                                                  >  2:911576/1‑139 (MQ=255)                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 <  1:999307/139‑1 (MQ=255)                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  2:1139176/1‑139 (MQ=255)                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  2:513960/1‑139 (MQ=255)                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑cccccc                                                                                                                             >  1:84632/1‑63 (MQ=255)                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑cccccc                                                                                                                             <  2:84632/63‑1 (MQ=255)                                                   gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                              <  2:441054/139‑1 (MQ=255)                                                    aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                               <  2:646185/137‑1 (MQ=255)                                                    aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                               >  1:646185/1‑137 (MQ=255)                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                      >  2:276284/1‑88 (MQ=255)                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                      <  1:276284/88‑1 (MQ=255)                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc                                   >  2:1198838/1‑139 (MQ=255)                                                                  ttGCTTTTGCGTTTACTGCCATAAACGCCTCGCGAACCCC‑‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  2:1254814/139‑1 (MQ=255)                                                                  ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  2:1389810/139‑1 (MQ=255)                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                                                            <  2:575088/87‑1 (MQ=255)                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                                                            >  1:575088/1‑87 (MQ=255)                                                                            gTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:513960/139‑1 (MQ=255)                                                                               tACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAg                  >  2:1444557/1‑139 (MQ=255)                                                                                       tAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTGGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGGGAATATTTTTCGGATTATCCGGAAAACAGGAGACCGTCTCTTACAGCACGCCTCCGCACAGTCCACgcc          >  2:960462/1‑139 (MQ=255)                                                                                       tAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               >  1:84289/1‑118 (MQ=255)                                                                                       tAAACGCCTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  2:84289/118‑1 (MQ=255)                                                                                              cTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgcca         <  2:745101/133‑1 (MQ=255)                                                                                              cTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgcca         >  1:745101/1‑133 (MQ=255)                                                                                              cTCGCGAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCcagccagc   >  2:1386640/1‑139 (MQ=255)                                                                                                cgcgAA‑CCTCCCCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACCCCTGCCCACAGTCCACGCCAGCCAGCa  >  2:1401615/1‑139 (MQ=255)                                                                                                 gcgAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                                   <  1:879614/88‑1 (MQ=255)                                                                                                 gcgAA‑CCTC‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                                   >  2:879614/1‑88 (MQ=255)                                                                                                          |                                                                                                                                   ACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCC‑ACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734432‑2734663

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