Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,352,500	(A)8→7	intergenic (+143/‑84)	yjdJ → / → ghoS	GNAT family putative N‑acetyltransferase/antitoxin of GhoTS toxin‑antitoxin pair; endonuclease for ghoT mRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,352,493	0	A	.	100.0%	281.7 / NA	59	intergenic (+136/‑91)	yjdJ/ghoS	GNAT family putative N‑acetyltransferase/antitoxin of GhoTS toxin‑antitoxin pair; endonuclease for ghoT mRNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (26/33);  total (26/33)

ATATTCGCAGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCAT  >  NC_000913/4352344‑4352638                                                                                                                                                      |                                                                                                                                                  atatTCGCAGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCa                                                                                                                                                   >  2:416339/1‑150 (MQ=255)   atTCGCAGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaa                                                                                                                                                <  1:94398/151‑1 (MQ=255)   atTCGCAGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaa                                                                                                                                                 >  2:426032/1‑150 (MQ=255)         aGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaaga                                                                                                                                           >  1:87700/1‑148 (MQ=255)          gTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaagaat                                                                                                                                         <  1:9207/151‑5 (MQ=255)           ttGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaagaat                                                                                                                                         >  2:1077667/1‑146 (MQ=255)            tGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAatat                                                                                                                                      >  1:898285/1‑151 (MQ=255)             gATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAatat                                                                                                                                      <  1:647748/150‑1 (MQ=255)              aTGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATAtct                                                                                                                                    <  1:715305/151‑1 (MQ=255)              aTGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATAtct                                                                                                                                    <  1:739313/151‑1 (MQ=255)                 ggAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCt                                                                                                                                 <  2:458273/151‑1 (MQ=255)                 ggAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCt                                                                                                                                 <  1:678487/151‑1 (MQ=255)                  gagaGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCta                                                                                                                                <  2:647588/151‑1 (MQ=255)                  gagaGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCta                                                                                                                                <  2:627176/151‑1 (MQ=255)                  gagaGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCt                                                                                                                                 >  2:728304/1‑150 (MQ=255)                  gagaGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCt                                                                                                                                 >  1:1044714/1‑150 (MQ=255)                      gTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATg                                                                                                                            <  2:561437/151‑1 (MQ=255)                         cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATg                                                                                                                            >  2:988359/1‑148 (MQ=255)                            agTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGAATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATTAGAatc                                                                                                                      <  1:1192611/151‑1 (MQ=255)                                aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGGAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAAtcatca                                                                                                                  <  2:995720/151‑1 (MQ=255)                                   atatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATc                                                                                                               >  1:197069/1‑151 (MQ=255)                                          tCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                                                                                                         >  2:232304/1‑150 (MQ=255)                                           cATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaa                                                                                                       <  2:275967/151‑1 (MQ=255)                                            aTTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaat                                                                                                      >  2:109030/1‑151 (MQ=255)                                                tctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGt                                                                                                  >  1:930710/1‑151 (MQ=255)                                                 ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGtt                                                                                                 >  1:381368/1‑151 (MQ=255)                                                  tcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGttt                                                                                                >  2:1018959/1‑151 (MQ=255)                                                    cgcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg                                                                                              <  1:75743/151‑1 (MQ=255)                                                     gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGc                                                                                             <  1:728294/151‑1 (MQ=255)                                                     gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGc                                                                                             <  2:498964/151‑1 (MQ=255)                                                     gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGc                                                                                             <  1:633521/151‑1 (MQ=255)                                                       gATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGt                                                                                           <  2:32587/151‑1 (MQ=255)                                                          gtgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcc                                                                                         <  2:1045628/150‑1 (MQ=255)                                                             aTGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcccc                                                                                       >  1:834146/1‑149 (MQ=255)                                                               gCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCagag                                                                                   <  1:187204/151‑1 (MQ=255)                                                                    agagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGttt                                                                              <  2:430726/151‑1 (MQ=255)                                                                             gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGAt                                                                     >  2:373285/1‑151 (MQ=255)                                                                              aaGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATAtt                                                                        <  1:88748/147‑1 (MQ=255)                                                                              aaGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGAt                                                                     <  2:630699/150‑1 (MQ=255)                                                                              aaGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGAt                                                                     >  1:920678/1‑150 (MQ=255)                                                                                 ggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGa                                                                 <  1:691574/151‑1 (MQ=255)                                                                                  gTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAgg                                                                  >  2:785536/1‑149 (MQ=255)                                                                                    aGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTc                                                              >  2:149409/1‑151 (MQ=255)                                                                                      ggggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCat                                                            <  2:319892/151‑1 (MQ=255)                                                                                      ggggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCat                                                            >  2:893515/1‑151 (MQ=255)                                                                                         gATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCAtatt                                                         <  1:975921/151‑1 (MQ=255)                                                                                            gCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATg                                                      >  2:604399/1‑151 (MQ=255)                                                                                                  aGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGt                                                <  1:135282/151‑1 (MQ=255)                                                                                                               agCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGATAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCa                                   <  1:439266/151‑1 (MQ=255)                                                                                                               agCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTc                                    >  2:1122886/1‑150 (MQ=255)                                                                                                               agCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTc                                    >  2:1008801/1‑150 (MQ=255)                                                                                                                 cgcgTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAAt                                 <  1:293826/151‑1 (MQ=255)                                                                                                                         tttGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgt                         >  1:481591/1‑151 (MQ=255)                                                                                                                           tGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCa                                                                                      <  1:211927/88‑1 (MQ=255)                                                                                                                           tGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCa                                                                                      >  2:211930/1‑88 (MQ=255)                                                                                                                            gTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgttgt                      <  1:232302/151‑1 (MQ=255)                                                                                                                             tCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgtt                        <  2:867509/148‑1 (MQ=255)                                                                                                                                cAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGT‑AAAACAAGTTCAATTCTTATGTTGTTTCtt                 >  2:1144093/1‑151 (MQ=255)                                                                                                                                     cATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGa             <  2:930728/151‑1 (MQ=255)                                                                                                                                      aTCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGAt            >  1:655543/1‑151 (MQ=255)                                                                                                                                         cccTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTAt         >  2:849619/1‑151 (MQ=255)                                                                                                                                            tttACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATcc       >  2:1044486/1‑150 (MQ=255)                                                                                                                                               aCTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCcatca   <  2:821296/151‑1 (MQ=255)                                                                                                                                               aCTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCcatca   <  2:928376/151‑1 (MQ=255)                                                                                                                                                cTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCcatcat  >  2:1164867/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ATATTCGCAGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCAT  >  NC_000913/4352344‑4352638

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