Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	243,207	A→G	W33R (TGG→CGG)	fadE ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	243,207	0	A	G	97.6%	141.3 / NA	41	W33R (TGG→CGG)	fadE	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/22);  minor base C (0/1);  total (18/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 = 1.00e+00

ACCGCGAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCTGACCACTTGT  >  NC_000913/243079‑243343                                                                                                                                 |                                                                                                                                         aCCGCGAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAATGCCCAGGGCGGCTCTCCGGGCGAGCaaaa                                                                                                                               >  1:258415/1‑140 (MQ=255)     cgAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAg                                                                                                                           >  1:551088/1‑140 (MQ=255)     cgAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAg                                                                                                                           >  2:160968/1‑140 (MQ=255)      gAAATACCGGCGCGGAAATCATTGACTTACGTATATGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGa                                                                                                                          <  2:288002/140‑1 (MQ=255)          tACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGActgc                                                                                                                      <  1:87351/140‑1 (MQ=255)            ccGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGActgctg                                                                                                                    <  2:479360/140‑1 (MQ=255)                 gcgGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGataaa                                                                                                               <  1:105041/140‑1 (MQ=255)                  cgGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGataaat                                                                                                              >  1:510239/1‑140 (MQ=255)                       atcatcGACTTACGCTTAGGCGCAATGTTTAATGGCATGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCt                                                                                                         <  1:236385/140‑1 (MQ=255)                       atcatcGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGTAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCt                                                                                                         <  1:88883/140‑1 (MQ=255)                       atcatcGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCt                                                                                                         <  1:78507/140‑1 (MQ=255)                           tcGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCAcg                                                                                                     <  1:578175/140‑1 (MQ=255)                           tcGACTTACGCATAGACGCAAACTTCAATTGCACGAGCATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCAcg                                                                                                     <  1:170343/140‑1 (MQ=255)                            cGACTTACGCATAGGCGCAAAGCTAAATGGCGCGACGATTATGGCCAGTGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCAcgc                                                                                                    <  2:176132/140‑1 (MQ=255)                                        aGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGaaca                                                                                        >  2:657626/1‑140 (MQ=255)                                              aaaGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGAATGCGGATACATAAGCTCACGTACTGATAGAACAACgcgc                                                                                  >  2:313349/1‑140 (MQ=255)                                                  ttAAATGGCACGAGGATCATGGTCAGAGGCACCAGCACCTACGCCGACTACAGACCAGCAACGCCGAGGGCTGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAg                                                                              <  2:276292/140‑1 (MQ=255)                                                  ttAAATGGCACGAGGATAATGGGCATAGGTACCAGTATCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAg                                                                              <  1:285908/140‑1 (MQ=255)                                                          cACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACaa                                                                      >  1:496452/1‑140 (MQ=255)                                                           aCGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAAc                                                                     >  2:104875/1‑140 (MQ=255)                                                                     tGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCgaga                                                           >  1:492437/1‑140 (MQ=255)                                                                       gCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAAt                                                         >  2:98144/1‑140 (MQ=255)                                                                                          aCGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACg                                      >  1:61069/1‑140 (MQ=255)                                                                                            gCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGaa                                    <  1:104875/140‑1 (MQ=255)                                                                                            gCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGaa                                    >  1:126035/1‑140 (MQ=255)                                                                                                   acaGACCAGCAACGCTGAGGGCTGCAGTCGGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGTGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGcccc                             <  1:313349/140‑1 (MQ=255)                                                                                                    caGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCt                            <  2:126035/140‑1 (MQ=255)                                                                                                    caGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCt                            <  2:664247/140‑1 (MQ=255)                                                                                                    caGACCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCt                            <  1:657626/140‑1 (MQ=255)                                                                                                       aCCAGCAACGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTAc                         >  2:416946/1‑140 (MQ=255)                                                                                                             aaCGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTag                   <  1:160968/140‑1 (MQ=255)                                                                                                              aCGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTagg                  <  1:574518/140‑1 (MQ=255)                                                                                                               cGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTagga                 >  1:540375/1‑140 (MQ=255)                                                                                                               cGCCGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTagga                 >  1:166029/1‑140 (MQ=255)                                                                                                                 ccGAGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACACAAAATCGTCATAACGGAAAGACCCTTACTTGTaggagg               >  2:169282/1‑140 (MQ=255)                                                                                                                  cGAGGGCGGCTGTCCCGGCGAGCAATATCTGACTGCTCATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGt              <  1:169282/140‑1 (MQ=255)                                                                                                                    aGGGCGGCTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAGAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCt            >  2:375753/1‑140 (MQ=255)                                                                                                                      ggcggcTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCTGa          >  1:62583/1‑140 (MQ=255)                                                                                                                          gcTGTCCGGGCGAGCTAAATCACACTGCTAATAAATATGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCTGACCAc      <  1:168120/140‑1 (MQ=255)                                                                                                                           cTGTCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCTGACCACt     <  2:505718/140‑1 (MQ=255)                                                                                                                              tCCGGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAGAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCTGACCACTtgt  <  1:375753/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                         ACCGCGAAATACCGGCGCGGAAATCATCGACTTACGCATAGGCGCAAAGTTAAATGGCACGAGGATAATGGCCAGAGGCACCAGTACCCACGCCGACCACAGACCAGCAACGCCGAGGGCGGCTGTCCAGGCGAGCAAAATCAGACTGCTGATAAATAAGCTCACGCGGTGATAGAACAACGCGCCGAGCAGGACAACCGTAGCGAGAATACTCAAAATCATCATAACGAAAAGCCCCTTACTTGTAGGAGGTCTGACCACTTGT  >  NC_000913/243079‑243343

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