Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	51.6%	52.0 / 109.7	62	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (5/25);  new base T (32/0);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 9.67e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-12

ACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT  >  NC_000913/1413595‑1413883                                                                                                                                             |                                                                                                                                                       aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAaca                                                                                                                                             <  2:584415/151‑1 (MQ=255) aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTAAA‑TTaa                                                                                                                                               >  1:370165/1‑149 (MQ=255) aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTTACacttaa                                                                                                                                               >  2:792306/1‑145 (MQ=255)   tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TTAacaac                                                                                                                                           >  1:300173/1‑149 (MQ=255)     tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCA‑TGAACACCTc                                                                                                                                         >  2:805678/1‑151 (MQ=255)        gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtttaaaataaaaaccccc                                                                                                                                        >  1:258181/1‑136 (MQ=255)        gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCa                                                                                                                                       <  2:512789/150‑1 (MQ=255)           ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑AGAAAACCTCCcttat                                                                                                                                   >  2:69794/1‑151 (MQ=255)              tGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATA‑CTAACACCACCattttttt                                                                                                                                >  1:351620/1‑149 (MQ=255)               ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCattattatt                                                                                                                               <  1:911878/151‑1 (MQ=255)                 aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTc                                                                                                                             <  2:286325/151‑1 (MQ=255)                      ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCA‑TTAAAACCTCCCTtttttttttcctt                                                                                                                         >  2:611310/1‑145 (MQ=255)                                 tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                                                                                                              <  1:940010/150‑1 (MQ=255)                                  ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTA‑TGAACACCCCCCCTATTTTTTTCACTGGtttggttttt                                                                                                             >  2:957931/1‑142 (MQ=255)                                  ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAA‑TTAACACCTCCCTTATTATTTTCAGTGGGTCGtttttt                                                                                                             >  1:654637/1‑150 (MQ=255)                                  ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                                                                                                             <  2:222931/150‑1 (MQ=255)                                    aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGATTTTTTATA‑TGAAAACCTCCCCTATTTTTTTTCGTGGGTCGttttttt                                                                                                            >  2:802922/1‑148 (MQ=255)                                    aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACA‑AGAACACCCCCCTTATTATTTTCAGTGGGTCGtttttttc                                                                                                           >  1:425509/1‑148 (MQ=255)                                        gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAA‑TGAAAACCCCCCTTATTTTTTTCAGTGGGTCGGTTTTTTCATCtt                                                                                                      >  1:341364/1‑151 (MQ=255)                                           agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                                                                                                    <  1:462813/150‑1 (MQ=255)                                            gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCTCCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCgt                                                                                                  <  2:714959/151‑2 (MQ=255)                                             cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                                                                                                 <  1:827609/151‑1 (MQ=255)                                              aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAA‑TTAAAACCCCCCTTATTTTTTTCAGTGGGTTGttttttttattttttttgt                                                                                                >  2:538736/1‑138 (MQ=255)                                                tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAA‑TTAAAACCCCCCTTTTTTTTTTCTGTGGTTTGTTTTTTTCATCTTTTGAGTg                                                                                               >  2:259769/1‑150 (MQ=255)                                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCCATTTTTTTCA‑TGAACACCTCCATTATTTTTTTCCGTGGGTCGTTTTTTTCATCTTTTTAGTGcttc                                                                                           >  1:84109/1‑151 (MQ=255)                                                    tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttc                                                                                           <  1:467497/150‑1 (MQ=255)                                                       cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAAACCCCCATT‑TTATTTTCCATGGGTTGTTTTTTTCATCTTTTTAGTGGTtttttt                                                                                       >  1:607951/1‑151 (MQ=255)                                                          tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑TGAACAACTCCCTTTTTATTTTCAGTGGTTTGTTTTTTTTCTCTTTCTAGGGGTTCttttt                                                                                      >  2:31077/1‑149 (MQ=255)                                                           cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTA‑TTAAAACCCCCCCTTTTTTTTTCATTGGGTTGTTTTTTTCATTTTTTCTGTGGttttttttt                                                                                     >  2:599522/1‑148 (MQ=255)                                                             gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TGAACACCCCCCTTTTTTTTTTCTATGGGTCGTTTTTTTCATCTTTTCAGTGCTTCtttttt                                                                                     >  1:449273/1‑146 (MQ=255)                                                             gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc                                                                                   <  1:775182/149‑1 (MQ=255)                                                               tttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TTAACACCCCCCTTATTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGGGCttttttttttctt                                                                                 >  1:335742/1‑144 (MQ=255)                                                                 tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑TAAAAAACTCCCTTATTATTTTCATTGGGTTTTTTTTTTCATCTTTTTCGTGCTTCtttttttctttccc                                                                             >  1:753480/1‑142 (MQ=255)                                                                 tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTcca                                                                             <  2:503147/150‑1 (MQ=255)                                                                     tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAA‑TGAAAACCCCCCTTATTTTTTTTAGTGGGTCGGTTATTTCATTTTTTTAGTGCTTCtttttttctttc                                                                               >  1:746425/1‑138 (MQ=255)                                                                        tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑TGAAAAACCCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGCttttttttttctttcccc                                                                            >  1:266075/1‑135 (MQ=255)                                                                              ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCATCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAA‑TTAACACCTCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTTCTCTTTTTAGTGGTTTTTTTTTTCTTTCCCCCTTAACCgtt                                                                  >  2:909941/1‑147 (MQ=255)                                                                                    atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAA‑TTAAAAACCCCCTTTTTTTTTTCCGTGGGTCGGTTTTTTCCTCTTTTGAGGGGTTTTTTTTTTC‑TTTCCCCCTAACCGGGTTGGGGtt                                                          >  1:302790/1‑150 (MQ=255)                                                                                          tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCa                                                      <  2:613908/148‑1 (MQ=255)                                                                                          tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAA‑AGAACAACCCCCCTATTATTTTCAGTGGGTTGTTTTTTTCCTCTTTTGTGTGCTTTTTTTTTTC‑TTTCCCCCTAAACCGGTCgggggttgccttc                                                   >  1:607359/1‑141 (MQ=255)                                                                                           gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATcc                                                   <  1:431189/150‑1 (MQ=255)                                                                                              cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑AGAAAACCCCCCCTTTTTTTTTCCGTGGGTTGTTTTTTTTCTCTTTTCTGTGCTTCTTTTTTTC‑TTCTCCCATAAACCGTTCcggggttgcctccctt                                                >  1:515545/1‑137 (MQ=255)                                                                                                  tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TGAAAAACCCCCTTATTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTGAGTGCTTCTTTTTTTC‑TTTCCCCCTAAACGGGTTGGGGTTTGCATCCgcttca                                             >  1:740402/1‑144 (MQ=255)                                                                                                  tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTATAAA‑TGAAAACCCCCCCTTTTTTTTTTAGTGGTTTGTTTTTTTCCTTTTTTCTGTGCTTTTTTTTCAC‑TTTCCCCCAAAACGGTTCGGGATTCGCCTCCGtt                                                >  2:541931/1‑146 (MQ=255)                                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                         <  2:266852/151‑1 (MQ=255)                                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                         <  2:197024/151‑1 (MQ=255)                                                                                                      agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAA‑TGAAAACCCCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCTTCTTTCGAGTGCTTTTTTTTTTC‑TTCCACCATAACCCGGTTGGGGGTCGCATTCCTTaca                                             >  1:720411/1‑145 (MQ=255)                                                                                                      agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCg                                        <  2:73337/150‑1 (MQ=255)                                                                                                          gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg                                   <  1:469909/151‑1 (MQ=255)                                                                                                            gCCACCTTCCTTTTCAATAGTGGCGGCCATTTTTTTCA‑AGAACACCTCCCCTATTTTTTCTTGGGGGTCGGTTTTTTCATCTTTCTGGTGCTTTTTTTTTAT‑TTTCACCACAACCCGATTGGGGGTCCCATTCGTTaaac                                            >  1:363811/1‑140 (MQ=255)                                                                                                                ccttccttTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAata                              <  2:654566/150‑1 (MQ=255)                                                                                                                        ttCAATAGTGGCGGTAATTTTTTAAA‑TGAAAACCTCCCTTTTTTTTTTTAGTGGTTCGGTTTTTTCTTTTTTTTGGTGCTTTTTTTTTTC‑TTTCCCCATAAACGGTTTCGGGTTTGCATCCGTTTCACTTCCCCCGGTAatattccgcaa                      >  2:957727/1‑143 (MQ=255)                                                                                                                         tCAATAGTGGCGGTAATTTTTTAAA‑TGAAAACCCCCCCTTTTTTTTTCAGTGGGTCGTTTTTTTCATTTTTTTAGTGCTTTTTTTTTTC‑TTTCACCATAACCCGTTTGGGGTTTGCCTCCGTTTCACTTCTCTCGGTAATATTCTGaaa                      >  2:400767/1‑149 (MQ=255)                                                                                                                          cAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATg                    <  1:805760/150‑1 (MQ=255)                                                                                                                              aGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                <  2:667769/150‑1 (MQ=255)                                                                                                                                  gCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc           >  2:621973/1‑151 (MQ=255)                                                                                                                                  gCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc           >  2:410655/1‑151 (MQ=255)                                                                                                                                  gCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAca           <  2:495778/151‑2 (MQ=255)                                                                                                                                    ggTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt         >  2:504674/1‑151 (MQ=255)                                                                                                                                     gTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt        <  2:247666/151‑1 (MQ=255)                                                                                                                                     gTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt        >  1:445801/1‑151 (MQ=255)                                                                                                                                       aaTTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc      <  1:640773/151‑1 (MQ=255)                                                                                                                                           ttATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGa   <  2:503572/150‑1 (MQ=255)                                                                                                                                            tataCA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt  <  1:370494/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                       ACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT  >  NC_000913/1413595‑1413883

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