Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,599,757	A→G	F2L (TTC→CTC)	dapA ←	4‑hydroxy‑tetrahydrodipicolinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,599,757	0	A	G	97.9%	163.6 / ‑4.3	47	F2L (TTC→CTC)	dapA	4‑hydroxy‑tetrahydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (19/27);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACC  >  NC_000913/2599633‑2599886                                                                                                                             |                                                                                                                                  ccAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGcc                                                                                                                           <  2:571739/133‑1 (MQ=255) ccAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGcc                                                                                                                           >  2:432291/1‑133 (MQ=255)      aGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCct                                                                                                                      <  2:639115/133‑1 (MQ=255)                gCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAAc                                                                                                            >  2:55558/1‑133 (MQ=255)                      gTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAgtgt                                                                                                      <  2:19161/133‑1 (MQ=255)                             tGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATg                                                                                               <  1:423503/133‑1 (MQ=255)                              ggCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATgg                                                                                              >  2:69659/1‑133 (MQ=255)                              ggCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATgg                                                                                              <  2:486029/133‑1 (MQ=255)                                  aCATGATACTCAATCAGTCTTTCCAAGCTAACCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACg                                                                                          <  1:180190/133‑1 (MQ=255)                                        tAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGt                                                                                    <  1:69659/133‑1 (MQ=255)                                             aatcaGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTATGTTTGGTATGGc                                                                               >  1:96414/1‑133 (MQ=255)                                                caGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCAtt                                                                            >  2:126895/1‑133 (MQ=255)                                                caGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCAtt                                                                            >  1:447764/1‑133 (MQ=255)                                                        tCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaa                                                                    <  1:55558/133‑1 (MQ=255)                                                        tCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaa                                                                    <  1:121556/133‑1 (MQ=255)                                                         cAAGCTAGCCCGACAAACATTACCCTTTGCATTCATCGGAGGAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaag                                                                   <  2:260055/133‑1 (MQ=255)                                                             cTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGcaga                                                               >  2:121237/1‑133 (MQ=255)                                                             cTAGCCCGACAGACATTACCTTTTACATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGcaga                                                               >  2:591314/1‑133 (MQ=255)                                                                 cccGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAGAGCAAGcagacaga                                                           >  2:480599/1‑133 (MQ=255)                                                                 cccGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGcagacaga                                                           <  1:126895/133‑1 (MQ=255)                                                                 cccGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGcagacaga                                                           >  2:527204/1‑133 (MQ=255)                                                                 cccGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGGTTGGTATGGCATTAAAAGCAATcagacaga                                                           >  1:260055/1‑133 (MQ=255)                                                                  ccGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGaa                                                          <  1:121237/133‑1 (MQ=255)                                                                    gacagacaTTGCCTTTTTCATCCATCGGTGTATCAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAAcc                                                        <  2:281679/133‑1 (MQ=255)                                                                            ttACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                <  1:591314/133‑1 (MQ=255)                                                                            ttACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                <  2:518430/133‑1 (MQ=255)                                                                              aCCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGa                                                 >  1:142354/1‑130 (MQ=255)                                                                              aCCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGa                                                 <  2:142354/130‑1 (MQ=255)                                                                                   tttCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGTACGGTTCTGATTGTTGTa                                         >  1:204471/1‑133 (MQ=255)                                                                                       atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCa                                     <  1:66503/133‑1 (MQ=255)                                                                                         ccatcGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATg                                   >  2:17611/1‑133 (MQ=255)                                                                                          catcGGAGGAACAATCGCCACATTACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGt                                  <  1:239101/133‑1 (MQ=255)                                                                                           atcGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                 >  2:372558/1‑133 (MQ=255)                                                                                             cGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACTGAACCGTTCTGCTTGTTGTATGCATGtttt                               >  1:306061/1‑133 (MQ=255)                                                                                             cGGAGTAACAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                               >  2:10450/1‑133 (MQ=255)                                                                                                    aCAATCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                        <  1:432291/133‑1 (MQ=255)                                                                                                      aaTCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCt                      <  2:365063/133‑1 (MQ=255)                                                                                                        tCGCGACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCttt                    <  1:534043/133‑1 (MQ=255)                                                                                                         cgcgACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTc                   >  1:639483/1‑133 (MQ=255)                                                                                                          gcgACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTTATTGTCGCATGCATGTATTTTATATGCTTTcc                  >  2:292302/1‑133 (MQ=255)                                                                                                            gACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAGAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCtt                <  1:480599/133‑1 (MQ=255)                                                                                                            gACAATACCTCCCGTGAGTATGTGCCATCCCCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCtt                <  1:292302/133‑1 (MQ=255)                                                                                                               aaTACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAg             <  1:527204/133‑1 (MQ=255)                                                                                                                    ttCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGTATGTTTTGTTTATGCTTTCCTTAAGaacaa        >  1:258154/1‑133 (MQ=255)                                                                                                                       ccGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTc     <  1:372558/133‑1 (MQ=255)                                                                                                                        cttGAGCATGGGCCATCCTCTGTCCATAAAAGTGTCTCAATTGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa    <  2:306061/131‑1 (MQ=255)                                                                                                                          tGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAcc  <  1:10450/133‑1 (MQ=255)                                                                                                                             |                                                                                                                                  CCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACC  >  NC_000913/2599633‑2599886

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