Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,462,678	T→A	11.3%	intergenic (+32/‑334)	fadL → / → yfdF	long‑chain fatty acid outer membrane transporter/uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,462,678	0	T	A	11.3%	139.6 / 7.7	44	intergenic (+32/‑334)	fadL/yfdF	long‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base A (2/3);  total (21/23)
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.92e-01

GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAT  >  NC_000913/2462530‑2462821                                                                                                                                                     |                                                                                                                                                 gACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCAt                                                                                                                                                  >  4:331145/1‑149 (MQ=255)     ttGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                             <  2:133139/149‑1 (MQ=255)       ggTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAgg                                                                                                                                           <  3:331145/149‑1 (MQ=255)       ggTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAgg                                                                                                                                           <  1:759522/149‑1 (MQ=255)        gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGt                                                                                                                                          <  3:515912/149‑1 (MQ=255)        gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                          <  3:428656/149‑1 (MQ=255)        gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                          <  4:661312/149‑1 (MQ=255)        gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                          <  4:410363/149‑1 (MQ=255)          gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                        <  1:508393/149‑1 (MQ=255)            ttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGACt                                                                                                                                      <  3:964624/149‑1 (MQ=255)               ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTa                                                                                                                                   >  1:541768/1‑149 (MQ=255)                                  gTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                >  2:1110977/1‑149 (MQ=255)                                             cGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCt                                                                                                     >  1:190466/1‑149 (MQ=255)                                                                   tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                               >  1:872800/1‑149 (MQ=255)                                                                   tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                               <  1:1110977/149‑1 (MQ=255)                                                                       aaGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaa                                                                           >  4:629826/1‑149 (MQ=255)                                                                         ggTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAAc                                                                         >  1:23923/1‑149 (MQ=255)                                                                           tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAg                                                                       <  4:98673/149‑1 (MQ=255)                                                                            aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                      <  2:541768/149‑1 (MQ=255)                                                                                 cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATAGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                 >  4:936910/1‑149 (MQ=255)                                                                                                 aaCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTc                                                 <  2:23923/149‑1 (MQ=255)                                                                                                           aCGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATaa                                       <  2:190466/149‑1 (MQ=255)                                                                                                             gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                     >  1:43827/1‑149 (MQ=255)                                                                                                             gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                     <  2:872800/149‑1 (MQ=255)                                                                                                               gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCgg                                   >  3:1068049/1‑149 (MQ=255)                                                                                                                    gATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                      >  4:583819/1‑125 (MQ=255)                                                                                                                    gATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                      <  3:583819/125‑1 (MQ=255)                                                                                                                              tCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAAc                    <  2:43827/149‑1 (MQ=255)                                                                                                                                  cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGa                >  1:702815/1‑149 (MQ=255)                                                                                                                                    ggATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAAc                                                                         <  4:230872/90‑1 (MQ=255)                                                                                                                                    ggATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAAc                                                                         >  3:230872/1‑90 (MQ=255)                                                                                                                                    ggATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg              <  4:1068049/149‑1 (MQ=255)                                                                                                                                      ttAAAGTCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAg                                                                                   >  4:481223/2‑78 (MQ=37)                                                                                                                                      ttAAAGTCACCTGC‑TATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAg                                                                                   <  3:481223/77‑1 (MQ=37)                                                                                                                                       tAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                              >  2:147575/1‑114 (MQ=255)                                                                                                                                       tAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                              <  1:147575/114‑1 (MQ=255)                                                                                                                                        aaaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa          >  4:777649/1‑149 (MQ=255)                                                                                                                                            tCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCtt                                                         <  2:233619/98‑1 (MQ=255)                                                                                                                                            tCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCtt                                                         >  1:233619/1‑98 (MQ=255)                                                                                                                                            tCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGc      >  1:354903/1‑149 (MQ=255)                                                                                                                                              aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat    >  4:813380/1‑149 (MQ=255)                                                                                                                                                cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                               >  4:212630/1‑120 (MQ=255)                                                                                                                                                cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                               <  3:212630/120‑1 (MQ=255)                                                                                                                                                cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat  >  3:127193/1‑149 (MQ=255)                                                                                                                                                cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat  >  1:279585/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAT  >  NC_000913/2462530‑2462821

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