Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,379,348	C→T	*432* (TAG→TAA)	menF ←	isochorismate synthase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,379,348	0	C	T	100.0%	123.2 / NA	39	*432* (TAG→TAA)	menF	isochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

AATGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTACTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCT  >  NC_000913/2379218‑2379480                                                                                                                                   |                                                                                                                                     aaTGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAACGATGCGACTCATTATTATTCCAt                                                                                                                              <  2:940454/139‑1 (MQ=255)                 cGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCattattatt                                                                                                                                  <  2:1262342/119‑1 (MQ=255)                 cGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCattattatt                                                                                                                                  >  1:1262342/1‑119 (MQ=255)                                    cTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTtgttgt                                                                                           >  2:1145061/1‑139 (MQ=255)                                         aTGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                                  >  2:665037/1‑127 (MQ=255)                                         aTGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                                  <  1:665037/127‑1 (MQ=255)                                                tGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGcc                                                                               >  1:555522/1‑139 (MQ=255)                                                  tCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCAc                                                                             >  1:561026/1‑139 (MQ=255)                                                      gTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGc                                                                                <  2:967833/132‑1 (MQ=255)                                                      gTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGc                                                                                >  1:967833/1‑132 (MQ=255)                                                                  cGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGta                                                                                                                         >  2:922915/1‑79 (MQ=255)                                                                  cGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGta                                                                                                                         <  1:922915/79‑1 (MQ=255)                                                                    gCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTTCGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTc                                                                     <  2:951321/129‑1 (MQ=255)                                                                    gCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTc                                                                     >  1:951321/1‑129 (MQ=255)                                                                       aaGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCa                                                        <  2:555522/139‑1 (MQ=255)                                                                               ggAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAAt                                                >  2:740033/1‑139 (MQ=255)                                                                                   cTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGccc                                            >  1:1128926/1‑139 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCa                                                        <  1:1227929/126‑1 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCa                                                        >  2:1227929/1‑126 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           >  1:347371/1‑139 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           >  1:1093105/1‑139 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           >  2:58095/1‑139 (MQ=255)                                                                                       aaaaaTAGCTTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        <  2:817656/139‑1 (MQ=255)                                                                                                   ttGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGAtt                                                                                      <  1:487651/81‑1 (MQ=255)                                                                                                   ttGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGAtt                                                                                      >  2:487651/1‑81 (MQ=255)                                                                                                      atatGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCa                         <  2:1128926/139‑1 (MQ=255)                                                                                                           aaTCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACa                    <  1:740033/139‑1 (MQ=255)                                                                                                             tCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACAtt                  <  2:241765/139‑1 (MQ=255)                                                                                                             tCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACAtt                  >  2:914113/1‑139 (MQ=255)                                                                                                                gTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                                  >  2:1226726/1‑56 (MQ=255)                                                                                                                gTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                                  <  1:1226726/56‑1 (MQ=255)                                                                                                                gTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        >  2:32322/1‑114 (MQ=255)                                                                                                                gTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        <  1:32322/114‑1 (MQ=255)                                                                                                                    tgatgCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        <  2:250053/110‑1 (MQ=255)                                                                                                                    tgatgCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        >  1:250053/1‑110 (MQ=255)                                                                                                                      atgCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAAt         >  2:301180/1‑139 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTc                                                                     >  1:1147185/1‑74 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTc                                                                     <  2:1147185/74‑1 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           <  2:300414/100‑1 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           >  1:300414/1‑100 (MQ=255)                                                                                                                             tCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTCCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCt  <  1:1145061/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     AATGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTACTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCT  >  NC_000913/2379218‑2379480

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