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	94.1%	107.0 / ‑0.2	34	coding (500/732 nt)	yfiH	UPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (18/14);  total (19/15)
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

AAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCA  >  NC_000913/2734403‑2734644                                                                                                                                      |                                                                                                              aaaaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑cccccc                                                                                                         <  1:712491/139‑1 (MQ=255)                   aCGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGt                                                                                        >  1:334596/1‑139 (MQ=255)                      ttCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGAc                                                                                    >  2:1364460/1‑139 (MQ=255)                              aCGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                             <  1:874113/139‑1 (MQ=255)                                 ataagagacaggCTGATAAATATACACCAGATACTGATCGCCATGCTGAATGACAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTACGAACGCGCGTGGACCAATTGCCgg                                                                          <  2:607588/128‑1 (MQ=255)                                  tGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGAcc                                                                        <  1:1364460/139‑1 (MQ=255)                                                      atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtt                                                     >  1:114624/1‑139 (MQ=255)                                                       tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCAACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTCCGAACGCGCGGGGACGGAATGCCGGCCCTAAACAGGccatccctcat                                                    >  1:607588/1‑129 (MQ=255)                                                          gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                              <  2:283226/94‑1 (MQ=255)                                                          gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                              >  1:283226/1‑94 (MQ=255)                                                             aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                              <  2:334596/139‑1 (MQ=255)                                                               aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGAtt                                            <  2:219006/139‑1 (MQ=255)                                                                   ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                                        >  1:1415647/1‑139 (MQ=255)                                                                     aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGc                                      >  2:1156625/1‑139 (MQ=255)                                                                        aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgcgcg                                                                                         >  1:1518329/1‑85 (MQ=255)                                                                        aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgcgcg                                                                                         <  2:1518329/85‑1 (MQ=255)                                                                        aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                                   >  1:1240229/1‑139 (MQ=255)                                                                         ccATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                             <  1:539635/96‑1 (MQ=255)                                                                         ccATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                             >  2:539635/1‑96 (MQ=255)                                                                            tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                               >  1:623992/1‑139 (MQ=255)                                                                                 aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                             >  2:64013/1‑88 (MQ=255)                                                                                 aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                             <  1:64013/88‑1 (MQ=255)                                                                                 aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc                          >  1:1167876/1‑139 (MQ=255)                                                                                       aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc                    >  2:175349/1‑139 (MQ=255)                                                                                           gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc                >  2:472313/1‑139 (MQ=255)                                                                                                      aGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGCGAATATTTTCTGGATTATCAGCAAACCAGGAAACCGTCTCTTCCAGCACGCCTgc     >  2:1354916/1‑139 (MQ=255)                                                                                                         gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca  <  1:472313/139‑1 (MQ=255)                                                                                                         gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca  <  2:489889/139‑1 (MQ=255)                                                                                                         gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca  <  1:175349/139‑1 (MQ=255)                                                                                                            tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                                   >  2:235560/1‑103 (MQ=255)                                                                                                            tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                                   <  1:235560/103‑1 (MQ=255)                                                                                                                    tAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCt       >  2:265292/1‑123 (MQ=255)                                                                                                                    tAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCt       <  1:265292/123‑1 (MQ=255)                                                                                                                           cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                >  2:779607/1‑91 (MQ=255)                                                                                                                           cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                <  1:779607/91‑1 (MQ=255)                                                                                                                                      |                                                                                                              AAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCA  >  NC_000913/2734403‑2734644

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