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%	296.3 / NA	62	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 . (22/40);  total (22/40)

AGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCA  >  NC_000913/4352352‑4352637                                                                                                                                              |                                                                                                                                                 aGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaaga                                                                                                                                          >  1:995606/1‑148 (MQ=255)    tGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAata                                                                                                                                      <  1:470849/150‑1 (MQ=255)     gATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAatat                                                                                                                                     <  1:476912/150‑1 (MQ=255)      aTGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATAtct                                                                                                                                   <  2:738414/151‑1 (MQ=255)             agTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTatat                                                                                                                            >  2:793769/1‑151 (MQ=255)                   gagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTCCTGAGC‑AAAAAAAGAATATCTCCTATATGAGAat                                                                                                                      <  2:857025/151‑1 (MQ=255)                      tCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAatcat                                                                                                                   >  2:543695/1‑151 (MQ=255)                         cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAGAAGAATATCTCCTATATGAGAAtcatca                                                                                                                 >  2:139292/1‑150 (MQ=255)                         cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAAtcatca                                                                                                                 >  1:851779/1‑150 (MQ=255)                          gATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAAt                                                                                                               >  2:725808/1‑151 (MQ=255)                           atatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATc                                                                                                              <  1:389053/151‑1 (MQ=255)                               ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCggg                                                                                                           <  1:1034224/150‑1 (MQ=255)                                tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCgggg                                                                                                          <  2:109603/150‑1 (MQ=255)                                         ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaataa                                                                                                   <  1:655314/148‑1 (MQ=255)                                             gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAAttt                                                                                               <  2:580513/148‑1 (MQ=255)                                              cgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg                                                                                           <  1:945526/151‑1 (MQ=255)                                               gATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg                                                                                           <  2:32935/150‑1 (MQ=255)                                               gATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGt                                                                                          >  1:145336/1‑151 (MQ=255)                                                  gtgATGCAGGAGAGCATCTGAAGGGTAGGTGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcc                                                                                        <  2:357273/150‑1 (MQ=255)                                                        cAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGc                                                                                 >  2:31845/1‑151 (MQ=255)                                                         aGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCg                                                                                <  2:653647/151‑1 (MQ=255)                                                         aGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCg                                                                                <  2:557891/151‑1 (MQ=255)                                                           gagaGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTGTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGtt                                                                              <  2:808341/151‑1 (MQ=255)                                                             gagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTa                                                                            <  1:462379/151‑1 (MQ=255)                                                                     gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTg                                                                      <  1:285190/149‑1 (MQ=255)                                                                       aGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGAt                                                                    >  2:577686/1‑149 (MQ=255)                                                                             gggggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTc                                                             >  2:485501/1‑150 (MQ=255)                                                                              ggggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCat                                                           <  2:941218/151‑1 (MQ=255)                                                                              ggggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCat                                                           <  1:281277/151‑1 (MQ=255)                                                                               gggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCata                                                          <  1:290450/151‑1 (MQ=255)                                                                                    gCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCAtattat                                                      <  1:543726/150‑1 (MQ=255)                                                                                    gCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATg                                                     >  1:591083/1‑151 (MQ=255)                                                                                         aaGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAgg                                                >  2:835173/1‑151 (MQ=255)                                                                                          aGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGt                                               <  2:233652/151‑1 (MQ=255)                                                                                           gactGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAGAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGt                                               >  2:637988/4‑150 (MQ=255)                                                                                            aaTGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTaa                                             <  1:200862/151‑1 (MQ=255)                                                                                            aaTGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTa                                              >  1:280720/1‑150 (MQ=255)                                                                                                 gCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACa                                        <  2:93451/151‑1 (MQ=255)                                                                                                     agagCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGtt                                    >  2:276889/1‑151 (MQ=255)                                                                                                     agagCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGtt                                    <  2:1021187/151‑1 (MQ=255)                                                                                                     agagCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGt                                     <  1:134960/150‑1 (MQ=255)                                                                                                       agCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTc                                   >  1:329970/1‑150 (MQ=255)                                                                                                         cgcgTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAAt                                <  1:906910/151‑1 (MQ=255)                                                                                                          gcgTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAAt                                >  2:1005876/1‑150 (MQ=255)                                                                                                          gcgTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAAt                                >  1:826362/1‑150 (MQ=255)                                                                                                            gTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACt                             >  2:48675/1‑151 (MQ=255)                                                                                                               tttttGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTa                           <  2:89459/150‑1 (MQ=255)                                                                                                               tttttGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTa                           <  2:409564/150‑1 (MQ=255)                                                                                                                   tGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgttg                      >  1:732560/1‑151 (MQ=255)                                                                                                                   tGTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgtt                       <  1:696257/150‑1 (MQ=255)                                                                                                                    gTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgttgt                     <  1:50638/151‑1 (MQ=255)                                                                                                                    gTCCCAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTAtgttgt                     <  1:554881/151‑1 (MQ=255)                                                                                                                      cccAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATgttgtt                    <  1:256349/150‑1 (MQ=255)                                                                                                                        cAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCt                 <  1:306835/151‑1 (MQ=255)                                                                                                                        cAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCt                 >  1:216100/1‑151 (MQ=255)                                                                                                                        cAAGTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCt                 <  2:1046440/151‑1 (MQ=255)                                                                                                                           gTCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCtttt              <  1:880341/151‑1 (MQ=255)                                                                                                                            tCATCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTg             >  1:840157/1‑151 (MQ=255)                                                                                                                              aTCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCTGAGCGTTTAATATTGATAGGAGTCATATTAGGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTGTTGAt           <  1:242627/151‑1 (MQ=255)                                                                                                                              aTCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGAt           >  2:161337/1‑151 (MQ=255)                                                                                                                               tCCCCTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCAGCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGa            <  2:766234/149‑1 (MQ=255)                                                                                                                                 cccTTTACTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTAt        <  1:683804/151‑1 (MQ=255)                                                                                                                                       aCTGAGC‑AAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCcatca  <  1:170068/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                 AGTTGATGGGAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCGTTTAATATTGATAGGAGTCATATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCA  >  NC_000913/4352352‑4352637

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