Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,150,120	T→C	F153F (TTT→TTC)	yghA →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,120	0	T	C	92.5%	161.8 / 1.1	53	F153F (TTT→TTC)	yghA	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (20/29);  total (22/31)
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

AGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCG  >  NC_000913/3150024‑3150216                                                                                                 |                                                                                                 aGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaa                                                                                               >  2:886530/1‑100 (MQ=255) aGCGCACAAGGCGTTAGGCGGGCTGGATATAATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaa                                                                                               >  1:430008/1‑100 (MQ=255)  gcgaACAATGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaa                                                                                              <  2:413170/100‑1 (MQ=255)     cacaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGa                                                                                           <  2:187860/100‑1 (MQ=255)     cacaAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGa                                                                                           >  1:1461200/1‑100 (MQ=255)            cGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGc                                                                                    >  1:1143935/1‑100 (MQ=255)             gTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGcc                                                                                   <  2:45974/100‑1 (MQ=255)                aGGCGGGCTGGATATTATGGCGCTGGTCGCCAGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCAtt                                                                                <  1:345348/100‑1 (MQ=255)                          gATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTcg                                                                      <  2:1575870/100‑1 (MQ=255)                            tattatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTcgcg                                                                    <  1:1025520/100‑1 (MQ=255)                            tattatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTcgcg                                                                    >  1:1187727/1‑100 (MQ=255)                             attatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCgcgc                                                                   >  2:753233/1‑100 (MQ=255)                                   gcgcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCt                                                             <  2:1329401/100‑1 (MQ=255)                                   gcgcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCt                                                             <  2:1006332/100‑1 (MQ=255)                                    cgcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTg                                                            <  1:1015630/100‑1 (MQ=255)                                    cgcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTg                                                            <  2:835708/100‑1 (MQ=255)                                      cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                          <  1:1234772/100‑1 (MQ=255)                                       tGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCt                                                         <  2:208048/100‑1 (MQ=255)                                        ggTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTa                                                        >  1:208048/1‑100 (MQ=255)                                         gTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTaa                                                       >  1:1334306/1‑100 (MQ=255)                                          tCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAc                                                      <  2:1355503/100‑1 (MQ=255)                                            gCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAccc                                                    >  2:61694/1‑100 (MQ=255)                                              cGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                   >  1:506307/1‑99 (MQ=255)                                               gggAAACGGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAgg                                                 <  2:468863/100‑1 (MQ=255)                                               gggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCGTTAACGTTTTCGCGCTGTTCTGGCTAACCCAgg                                                 <  2:1501276/100‑1 (MQ=255)                                               gggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAgg                                                 <  2:418651/100‑1 (MQ=255)                                               gggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAgg                                                 <  2:1501933/100‑1 (MQ=255)                                                  aaaCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAg                                              >  2:258892/1‑100 (MQ=255)                                                           gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATcccct                                     >  2:515325/1‑99 (MQ=255)                                                           gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccccc                                     <  1:1125249/100‑1 (MQ=255)                                                            ccATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCt                                    >  2:776429/1‑100 (MQ=255)                                                              aTTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                  <  2:1218951/100‑1 (MQ=255)                                                                tCCGGATATTGCAGACCTCTCCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                  <  2:1218432/98‑1 (MQ=255)                                                                 ccGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAc                               <  2:1143930/100‑1 (MQ=255)                                                                  cGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                              <  1:100559/100‑1 (MQ=255)                                                                  cGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                              >  1:1343088/1‑100 (MQ=255)                                                                   ggATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCg                             >  1:587362/1‑100 (MQ=255)                                                                         tGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGt                       >  2:1255787/1‑100 (MQ=255)                                                                          gCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                      >  1:675151/1‑100 (MQ=255)                                                                          gCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                      <  1:886535/100‑1 (MQ=255)                                                                            aGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCa                    <  2:14068/100‑1 (MQ=255)                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                  <  2:1432340/100‑1 (MQ=255)                                                                              aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                  <  1:61694/100‑1 (MQ=255)                                                                                    ccAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcat            >  2:1536735/1‑100 (MQ=255)                                                                                     cAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcatc           <  1:246478/100‑1 (MQ=255)                                                                                      aGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTAtcatca          >  1:670610/1‑100 (MQ=255)                                                                                       gCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcac         <  1:1272898/100‑1 (MQ=255)                                                                                        cgtaCAGTTTCGAAAGACCTTTGCCAGTAACGTTTTCGCGCTGTTCTGGCTATCCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc        <  2:1098955/97‑1 (MQ=255)                                                                                         gAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacca       >  1:199360/1‑100 (MQ=255)                                                                                          aaCAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcaccac      <  1:1271557/100‑1 (MQ=255)                                                                                             aGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTc   <  2:67807/100‑1 (MQ=255)                                                                                              gTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCg  >  2:128472/1‑100 (MQ=255)                                                                                              gTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCg  >  1:1298088/1‑100 (MQ=255)                                                                                                 |                                                                                                 AGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCG  >  NC_000913/3150024‑3150216

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