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

AATGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTACTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCT  >  NC_000913/2379218‑2379480                                                                                                                                   |                                                                                                                                     gaTGACCGACGCCCAGCGTCGGGTAAACGCGCTTACGGACATGAGCTTGTCCGGTATCAATATGGCGGCGAAGTAGAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                              <  1:423563/138‑1 (MQ=255) aaTGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                              <  2:1124777/139‑1 (MQ=255) aaTGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                              <  1:7221/139‑1 (MQ=255)       cgccgcCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGtaa                                                                                                                        >  1:1168130/1‑139 (MQ=255)         ccgccCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGtaata                                                                                                                      >  2:1478531/1‑139 (MQ=255)           gccCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATaaa                                                                                                                    <  2:1016333/139‑1 (MQ=255)            cccAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAg                                                                                                                   >  2:1159738/1‑139 (MQ=255)                   tCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAg                                                                                                            <  2:1168130/139‑1 (MQ=255)                                                      gTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCtt                                                                         >  2:689280/1‑139 (MQ=255)                                                               ttGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGt                                                                <  1:1478531/139‑1 (MQ=255)                                                                 gCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTttg                                                                                               <  1:1281707/106‑1 (MQ=255)                                                                 gCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTttg                                                                                               >  2:1281707/1‑106 (MQ=255)                                                                 gCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCg                                                              >  2:1147249/1‑139 (MQ=255)                                                                 gCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGATAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGTTCg                                                              >  1:1021680/1‑139 (MQ=255)                                                                   ggCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGa                                                            <  1:1147249/139‑1 (MQ=255)                                                                     cTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                              >  2:395241/1‑71 (MQ=255)                                                                     cTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCAt                                                                                                                              <  1:395241/71‑1 (MQ=255)                                                                     cTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAAc                                                          >  1:1527236/1‑139 (MQ=255)                                                                       aaGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCa                                                        >  1:1475887/1‑139 (MQ=255)                                                                               ggAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCAGCCg                                                                                                    <  2:189772/87‑1 (MQ=255)                                                                               ggAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCAGCCg                                                                                                    >  1:189772/1‑87 (MQ=255)                                                                                  gCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGcc                                             <  2:1181026/139‑1 (MQ=255)                                                                                   cTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGccc                                            >  2:1288909/1‑139 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           >  2:618340/1‑139 (MQ=255)                                                                                    tACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCcg                                           >  1:306619/1‑139 (MQ=255)                                                                                       aaaaaTAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        <  2:306619/139‑1 (MQ=255)                                                                                       aaaaaTAGATTATTGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        <  2:1527236/139‑1 (MQ=255)                                                                                                ttattGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCAtata                               >  1:1263472/1‑139 (MQ=255)                                                                                                 tattGATATGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTTGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATaa                              >  1:887990/1‑139 (MQ=255)                                                                                                         tGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        <  2:1361269/121‑1 (MQ=255)                                                                                                         tGAATCGGTAATGATGCGACTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        >  1:1361269/1‑121 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCg                                                              <  2:466349/81‑1 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCg                                                              >  1:466349/1‑81 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTg    >  1:912595/1‑139 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTg    >  2:657111/1‑139 (MQ=255)                                                                                                                           aCTCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTg    >  1:609784/1‑139 (MQ=255)                                                                                                                             tCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCt  <  2:1475887/139‑1 (MQ=255)                                                                                                                             tCATTATTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCt  <  2:1021680/139‑1 (MQ=255)                                                                                                                               attattattCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCgcgc                                        >  2:1034434/1‑99 (MQ=255)                                                                                                                               attattattCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCa                         >  1:1571186/1‑114 (MQ=255)                                                                                                                               attattattCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCa                         <  2:1571186/114‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     AATGACCGCCGCCCAGCGTCGGTTAAATGCGCTTACTGACATGAGTTTGTCCGGTATCAATATTGCGGCTAAGTATAAGGAGCTACAAAAATAGATTATTGATATGAATCGGTAATGATGCGACTCATTACTATTCCATTTGTAATAAAGTACGCAGCCCTGCCGCTTTGTTGTCGATTTCCTGCCACTCTTGCTCGGGGTCGGAACCACGGACAATGCCCGCGCCAGCATATAATCGCACGACATTGCCGCTAATTTTTGCT  >  NC_000913/2379218‑2379480

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