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%	134.5 / NA	42	*432* (TAG→TAA)	menF	isochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/21);  total (21/21)

AATGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTACTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCG  >  NC_000913/2379218‑2379468                                                                                                                                   |                                                                                                                         aaTGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                  <  1:388452/139‑1 (MQ=255)                gCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACg                                                                                                   >  1:625148/1‑139 (MQ=255)                     ggTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGcc                                                                                              <  2:625148/139‑1 (MQ=255)                     ggTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGcc                                                                                              <  2:43314/139‑1 (MQ=255)                                cTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTttg                                                                                   >  1:936324/1‑139 (MQ=255)                                              tttGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                      <  2:155828/122‑1 (MQ=255)                                              tttGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                      >  1:155828/1‑122 (MQ=255)                                              tttGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTg                                                                     <  2:936324/139‑1 (MQ=255)                                                tGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTg                                                                                                               <  2:189277/95‑1 (MQ=255)                                                tGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTg                                                                                                               >  1:189277/1‑95 (MQ=255)                                                                       aaGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCa                                            >  2:41496/1‑139 (MQ=255)                                                                               ggAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGAtt                                                                          <  2:1127984/101‑1 (MQ=255)                                                                               ggAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGAtt                                                                          >  1:1127984/1‑101 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                  <  1:1275007/56‑1 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                  >  2:1275007/1‑56 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGccc                                                                                             <  2:911986/77‑1 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGccc                                                                                             >  1:911986/1‑77 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGccc                                                                                             >  1:273002/1‑77 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGccc                                                                                             <  2:273002/77‑1 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCAct                                                                >  1:1252599/1‑106 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCAct                                                                <  2:1252599/106‑1 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTc                                                         >  1:379594/1‑113 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTc                                                         <  2:379594/113‑1 (MQ=255)                                                                                               attattGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCatat                    >  2:612008/1‑139 (MQ=255)                                                                                                 tattGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAg                         >  2:760306/1‑132 (MQ=255)                                                                                                 tattGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAg                         <  1:760306/132‑1 (MQ=255)                                                                                                    tGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCg               <  1:612008/139‑1 (MQ=255)                                                                                                    tGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCg               <  1:41496/139‑1 (MQ=255)                                                                                                           aaTCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGc              <  2:122621/133‑1 (MQ=255)                                                                                                           aaTCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGc              >  1:122621/1‑133 (MQ=255)                                                                                                             tCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACa        >  2:76884/1‑137 (MQ=255)                                                                                                             tCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACAtt      >  2:815477/1‑139 (MQ=255)                                                                                                             tCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACGCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACa        <  1:76884/137‑1 (MQ=255)                                                                                                                gTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGcc   >  2:900704/1‑139 (MQ=255)                                                                                                                       tgCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCg  >  2:476587/1‑133 (MQ=255)                                                                                                                       tgCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCg  <  1:476587/133‑1 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                      <  1:1401691/45‑1 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCt                                                                                      >  2:1401691/1‑45 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCa                          >  1:136626/1‑105 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCa                          <  2:136626/105‑1 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCg  >  2:1202918/1‑129 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCg  <  1:1202918/129‑1 (MQ=255)                                                                                                                                   |                                                                                                                         AATGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTACTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCG  >  NC_000913/2379218‑2379468

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