Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,789,056	T→C	F208L (TTC→CTC)	aroH →	3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,789,056	0	T	C	100.0%	273.5 / NA	73	F208L (TTC→CTC)	aroH	3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/47);  total (26/47)

GCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGAT  >  NC_000913/1788910‑1789202                                                                                                                                                   |                                                                                                                                                   gcgcGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCct                                                                                                                                                <  1:1159560/151‑1 (MQ=255)    cgATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCtctc                                                                                                                                             <  1:351460/151‑1 (MQ=255)     gATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCg                                                                                                                                            >  2:592151/1‑151 (MQ=255)     gATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCg                                                                                                                                            <  2:745067/151‑1 (MQ=255)     gATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCg                                                                                                                                            <  1:457994/151‑1 (MQ=255)     gATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCg                                                                                                                                            <  1:1195899/151‑1 (MQ=255)      atgggCGCACGTACTACCGAAAGGCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCGGGTAGGTTTTAAAAAGGGTACCGATGGGAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGc                                                                                                                                           <  1:704772/148‑1 (MQ=255)           cgcACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGAc                                                                                                                                      <  2:887088/151‑1 (MQ=255)            gcACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGAc                                                                                                                                      >  1:1197544/1‑150 (MQ=255)               cGTACTACCGAAAGTCAGATCCACCGCGAAATGGTTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACaaaa                                                                                                                                  <  1:363416/151‑1 (MQ=255)                    tacCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGt                                                                                                                             >  2:363640/1‑151 (MQ=255)                         aaaGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGAt                                                                                                                        <  1:947849/151‑1 (MQ=255)                          aaGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATg                                                                                                                       <  2:962131/151‑1 (MQ=255)                             tCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCa                                                                                                                                         <  2:1010935/130‑1 (MQ=255)                             tCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCa                                                                                                                                         >  1:1010911/1‑130 (MQ=255)                                gATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                                                                                                 >  2:1046478/1‑151 (MQ=255)                                   ccACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTa                                                                                                               <  2:192778/150‑1 (MQ=255)                                      ccGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATc                                                                                                             >  1:1086657/1‑149 (MQ=255)                                       cgggAAAGGGCTTCGGCACTGTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCGGGGCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAAGGGTCAGATGACTAt                                                                                                                  <  1:244667/140‑1 (MQ=255)                                         cgAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAcc                                                                                                        <  1:193453/151‑1 (MQ=255)                                         cgAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAcc                                                                                                        >  1:1098316/1‑151 (MQ=255)                                          gAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAcc                                                                                                        >  1:115491/1‑150 (MQ=255)                                          gAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                                                                                       <  1:548474/151‑1 (MQ=255)                                           aaaTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                                                                                      <  2:1003324/151‑1 (MQ=255)                                           aaaTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                                                                                      <  1:564249/151‑1 (MQ=255)                                             aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                    >  1:425692/1‑151 (MQ=255)                                             aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                    <  2:325105/151‑1 (MQ=255)                                               ggCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGc                                                                                                  >  1:235725/1‑151 (MQ=255)                                                gCTTCGGCCCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCgg                                                                                                   <  1:600133/149‑1 (MQ=255)                                                gCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                    >  2:1159236/1‑148 (MQ=255)                                                gCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCa                                                                                                 <  2:630723/151‑1 (MQ=255)                                                   tCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                                                                                                <  2:27179/149‑1 (MQ=255)                                                   tCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                                                                                               <  1:768141/150‑1 (MQ=255)                                                    cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                                                                                             <  1:126953/151‑1 (MQ=255)                                                    cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAcc                                                                                              <  1:744928/150‑1 (MQ=255)                                                     ggCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                                                                                            >  2:822322/1‑151 (MQ=255)                                                     ggCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                                                                                            <  2:494440/151‑1 (MQ=255)                                                        aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                         <  2:242700/151‑1 (MQ=255)                                                        aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                         <  1:87002/151‑1 (MQ=255)                                                           ctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGc                                                                                      >  2:863201/1‑151 (MQ=255)                                                                gTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca                                                                                  <  2:1000599/150‑1 (MQ=255)                                                                gTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACat                                                                                 >  1:1130776/1‑151 (MQ=255)                                                                    ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt                                                                             >  2:585436/1‑151 (MQ=255)                                                                       aGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat                                                                           <  1:425757/150‑1 (MQ=255)                                                                            ttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTg                                                                      <  2:958650/150‑1 (MQ=255)                                                                              aaaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTggcg                                                                   <  1:38708/151‑1 (MQ=255)                                                                              aaaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTggcg                                                                   >  1:304768/1‑151 (MQ=255)                                                                                 aaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTggcggc                                                                 >  2:149348/1‑150 (MQ=255)                                                                                  aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCa                                                                >  2:219511/1‑150 (MQ=255)                                                                                  aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCa                                                                >  2:28905/1‑150 (MQ=255)                                                                                    ggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaa                                                              <  1:300779/150‑1 (MQ=255)                                                                                      tACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                            >  1:72672/1‑150 (MQ=255)                                                                                       aCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAAc                                                          >  2:1021916/1‑151 (MQ=255)                                                                                        ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAAcc                                                         <  2:1197575/151‑1 (MQ=255)                                                                                              gCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAAtt                                                    >  2:765392/1‑150 (MQ=255)                                                                                               cAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTAt                                                  <  2:647839/151‑1 (MQ=255)                                                                                               cAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTAt                                                  <  2:551628/151‑1 (MQ=255)                                                                                                 aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATc                                                 >  2:601637/1‑150 (MQ=255)                                                                                                   aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                              <  2:324541/151‑1 (MQ=255)                                                                                                   aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                              >  2:966765/1‑151 (MQ=255)                                                                                                   aCGCGGATTGCGGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                              <  1:941577/151‑1 (MQ=255)                                                                                                         aTTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgat                                         <  1:754220/150‑1 (MQ=255)                                                                                                           tGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatga                                       <  2:524223/150‑1 (MQ=255)                                                                                                               gtgGATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGatac                                    <  1:350641/149‑2 (MQ=255)                                                                                                                  gATGCTATCCGCGCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATcgccg                               >  1:299426/1‑151 (MQ=255)                                                                                                                           cgcgCAGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCg                      <  1:747522/151‑1 (MQ=255)                                                                                                                                aGCCCGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACg                 <  2:888503/151‑1 (MQ=255)                                                                                                                                  cccGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCt               <  1:161929/151‑1 (MQ=255)                                                                                                                                  cccGCGCCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCt               <  1:1021891/151‑1 (MQ=255)                                                                                                                                     gcgcCAGCCATATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCa            <  1:118709/151‑1 (MQ=255)                                                                                                                                           gccaTATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGtt      <  2:299429/151‑1 (MQ=255)                                                                                                                                           gccaTATGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGtt      <  1:290408/151‑1 (MQ=255)                                                                                                                                               tatGCTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGAt  >  2:72058/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGAT  >  NC_000913/1788910‑1789202

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