Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,566,266	C→A	8.8%	G207C (GGC→TGC)	eutA ←	reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,566,266	0	C	A	8.8%	149.2 / 5.0	46	G207C (GGC→TGC)	eutA	reactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base A (2/2);  total (22/24)
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.33e-01

CCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGT  >  NC_000913/2566127‑2566406                                                                                                                                            |                                                                                                                                             ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGa                                                                                                                                      <  4:293138/148‑1 (MQ=255)  cTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAg                                                                                                                                    >  1:241866/1‑149 (MQ=255)  cTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGacg                                                                                                                                    <  1:402693/149‑3 (MQ=255)               aCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCGGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCcaca                                                                                                                        >  4:248249/1‑148 (MQ=255)                     ttGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatcatc                                                                                                                  >  4:45616/1‑148 (MQ=255)                         aTCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGccc                                                                                                             >  2:220889/1‑149 (MQ=255)                              tGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCa                                                                                                                                          <  3:124170/115‑1 (MQ=255)                              tGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCa                                                                                                                                          >  4:124170/1‑115 (MQ=255)                                 cTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTa                                                                                                      >  4:220266/1‑148 (MQ=255)                                 cTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTAt                                                                                                     <  1:190820/149‑1 (MQ=255)                                  tGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATg                                                                                                    <  2:19202/149‑1 (MQ=255)                                  tGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATg                                                                                                    <  2:282972/149‑1 (MQ=255)                                      cgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGc                                                                                                 <  4:292318/148‑1 (MQ=255)                                              gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAc                                                                                        >  1:356533/1‑149 (MQ=255)                                               agagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcg                                                                                       >  2:250666/1‑149 (MQ=255)                                                gagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcg                                                                                       <  3:59133/148‑1 (MQ=255)                                                    gTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGccc                                                                                  <  2:241866/149‑1 (MQ=255)                                                      tCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCg                                                                                 >  4:272752/1‑148 (MQ=255)                                                         gTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                             <  2:402278/149‑1 (MQ=255)                                                         gTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                             >  1:352222/1‑149 (MQ=255)                                                                 cTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGtt                                                                     >  1:233672/1‑149 (MQ=255)                                                                   tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                    >  3:68947/1‑148 (MQ=255)                                                                   tCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGttt                                                                    <  3:220266/148‑1 (MQ=255)                                                                    cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCGCCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCTTGCACCGAAGCATTCATCCAAAAGCATCTGCCCCGGTTTATGAGAGTAAACCACGCGCCCGTGGCTGTAGGttgct                                                                  <  2:397346/149‑4 (MQ=255)                                                                    cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTcc                                                                  <  2:52279/149‑1 (MQ=255)                                                                     gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCca                                                                 <  2:356533/149‑1 (MQ=255)                                                                           cAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCa                                                                                                                              <  3:286300/82‑1 (MQ=255)                                                                           cAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCa                                                                                                                              >  4:286300/1‑82 (MQ=255)                                                                           cAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAgg                                                            >  3:396280/1‑148 (MQ=255)                                                                                             ggTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGagg                                          >  3:186661/1‑148 (MQ=255)                                                                                               tCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggca                                        <  4:186661/148‑1 (MQ=255)                                                                                                  ccTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAggcaggca                                    >  2:281784/1‑149 (MQ=255)                                                                                                   cTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGgcaggcag                                   <  1:281784/149‑1 (MQ=255)                                                                                                      cACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTg                                 >  3:329456/1‑148 (MQ=255)                                                                                                           gCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCt                            <  4:329456/148‑1 (MQ=255)                                                                                                             tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGcc                               >  1:358329/1‑143 (MQ=255)                                                                                                             tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGcc                               <  2:358329/143‑1 (MQ=255)                                                                                                             tGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGa                          <  4:68947/148‑1 (MQ=255)                                                                                                                   ccGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTc                    <  3:272752/148‑1 (MQ=255)                                                                                                                          gCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCGAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCg            >  1:70811/1‑149 (MQ=255)                                                                                                                                gggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg      >  2:406469/1‑149 (MQ=255)                                                                                                                                gggCGTCGGTGCCTGCACCGAAACACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACATTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAg       <  3:248249/148‑1 (MQ=255)                                                                                                                                     tCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGt  <  4:328669/148‑1 (MQ=255)                                                                                                                                     tCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTg                           <  3:365582/123‑1 (MQ=255)                                                                                                                                     tCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTg                           >  4:365582/1‑123 (MQ=255)                                                                                                                                     tCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGt  >  3:256935/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                             CCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGT  >  NC_000913/2566127‑2566406

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: ‑