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

CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTA  >  NC_000913/1413588‑1413878                                                                                                                                                    |                                                                                                                                                 cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtata                                                                                                                                                <  1:609428/150‑1 (MQ=255)   gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑At                                                                                                                                            <  1:471528/151‑1 (MQ=255)     tttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtttaaa                                                                                                                                              >  2:454096/1‑146 (MQ=255)       tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAT‑ATTaa                                                                                                                                         >  1:126514/1‑150 (MQ=255)        aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAac                                                                                                                                        <  2:530290/150‑1 (MQ=255)           ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATTAacaac                                                                                                                                     >  1:700873/1‑148 (MQ=255)              tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATTAACAcccc                                                                                                                                   >  2:753786/1‑147 (MQ=255)                    ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCTACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCattat                                                                                                                             <  1:633467/149‑1 (MQ=255)                    ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTA‑ATTAATACCTCCatttt                                                                                                                             >  2:883047/1‑147 (MQ=255)                     tGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAT‑ATTAACACCTCCcttattat                                                                                                                          >  2:199306/1‑151 (MQ=255)                        aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑ATGAAAACcccccctttttttt                                                                                                                        >  2:393164/1‑140 (MQ=255)                         aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATTAAAACCCCCCTtttttttt                                                                                                                        >  1:288778/1‑149 (MQ=255)                                gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAT‑ATGAACAACTCCCTTTTTATTTTCAGt                                                                                                                   >  1:183795/1‑147 (MQ=255)                                   tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAT‑ATGAAAACCTCCCTTTTTTTTTTCCGTgggttgt                                                                                                            >  2:629132/1‑146 (MQ=255)                                     ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAA‑ATGAACACCTCCCTTATTATTTTCCGTGGGTCGttt                                                                                                          >  2:208402/1‑151 (MQ=255)                                       ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTA‑ATGAAAACCCCCCTTATTTTTTTCCGTGgtttgttttt                                                                                                        >  1:939982/1‑149 (MQ=255)                                        tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑ATTAACAACCCCCTTTTTTTTTTCCGTGGTTTGtttttt                                                                                                       >  2:7314/1‑151 (MQ=255)                                            ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATGAACACCTCCCTTATTTTTTTTAGTGGGTCGtttttt                                                                                                       >  1:221705/1‑147 (MQ=255)                                                agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTTGTAGGAGTGCCTCCTTCCTTTTCAATAGTGGCGGTAATTTTATTC‑TTGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAGCttt                                                                                               <  2:1007652/151‑1 (MQ=255)                                                agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                               <  1:48589/151‑1 (MQ=255)                                                 gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑ATTAACACCCCCATTATTTTTTTCAGTGGTTCGTTTTTTTCATCttt                                                                                               >  2:395147/1‑150 (MQ=255)                                                      ttttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATTAAAACCCCCCTTTTTTTTTTCAGTGGTTTGTTTTTTTCATCTTTTTAGTg                                                                                         >  1:448345/1‑151 (MQ=255)                                                       tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑AATAAAACCCCCCTTTTTATTTTCCGTGGGTTGTTTATTTCCTCTTTcttgtg                                                                                         >  2:94718/1‑145 (MQ=255)                                                             gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑ATGAACACCTCCCTTTTTATTTTCAGTGGGTTGTTTTTTTCATCTTTTTCGTGCtttttt                                                                                  >  1:700975/1‑151 (MQ=255)                                                               aTTCTGCTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAAAACATCCCTTATTAATTTCAGGGGTTCCGTTTTTCTCACTTTCGAGGGGCTCtttt                                                                                 >  2:730956/1‑150 (MQ=255)                                                               aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTA‑ATGAACACCCCCCTTTTTATTTTCCGTGGGTTGTTTTTTTTATCTTTTGAGTGCtttttttt                                                                                >  1:1074136/1‑151 (MQ=255)                                                                 tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAA‑ATGAAAACCTCCATTTTTATTTTTCGTGGGTCGTTTTTTTTATCTTTTTAGTGCttttttttt                                                                               >  1:888512/1‑149 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa                                                                              <  1:146690/148‑1 (MQ=255)                                                                     aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTA‑ATTAACAACCCCCTTATTTTTTTCAGTGGGTCGTTTTTTTCCTCTTTTTCGTGCTTCttttt                                                                                >  1:13102/1‑145 (MQ=255)                                                                     aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTc                                                                          <  1:867779/151‑1 (MQ=255)                                                                         cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccac                                                                       <  2:317472/150‑1 (MQ=255)                                                                         cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTccac                                                                       <  2:407579/150‑1 (MQ=255)                                                                          gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑ATGAAAACCCCCCTTTTTTTTTTCCGGGGGTTGTTTATTCCATCTTTTGAGTGCCTCTTTTTTACtttcccc                                                                      >  2:876128/1‑145 (MQ=255)                                                                           ttAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCCTTATTATTTTCAGTGGTTTGTTTTTTTCATCTTTTTCGTGCTTTTTTTTTTATTccacc                                                                      >  2:152825/1‑149 (MQ=255)                                                                            tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAt                                                                    <  1:760839/150‑1 (MQ=255)                                                                             aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑AATAACAGCTCCATTTTTATTTTCCGAGGGTCGTTTTTTCCTTCTTTTCCGTGGTTTTTTTTTACTTccccc                                                                      >  2:456888/1‑147 (MQ=255)                                                                                ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAcc                                                                <  2:802635/150‑1 (MQ=255)                                                                                 tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCg                                                               <  2:920962/150‑1 (MQ=255)                                                                                  aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCgg                                                              <  2:940005/150‑1 (MQ=255)                                                                                   gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGtt                                                            <  1:961939/151‑1 (MQ=255)                                                                                      tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCg                                                          <  2:826921/150‑1 (MQ=255)                                                                                          aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGAtt                                                     <  2:772178/151‑1 (MQ=255)                                                                                            tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTA‑AAGAAAACCCCCATTTTTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTTGTGCTTTTTTTTTACTTCCCCCCTAAACCGGTCggg                                                        >  1:1139310/1‑146 (MQ=255)                                                                                            tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAT‑ATGAAAACCCCCCTTTTTTTTTTCCGTGGGTCGTTTTTTTCTTCTTTTGTGTGCTTCTTTTTTTCTTCCACCATAACACGTTCGGGGtt                                                     >  1:1143231/1‑149 (MQ=255)                                                                                             atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAA‑ATGAACACCTCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCCTCTTTTTTGTGCTTTTTTTTTTCTTCCCCCCCAACCCGGTCggggt                                                      >  1:1021498/1‑145 (MQ=255)                                                                                             atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCg                                                   <  2:526195/150‑1 (MQ=255)                                                                                               cTTGGCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTAAACAATGAAAACC‑CTCTTAATAATTTTAGGGGGTTGGGTATTTCATCTCTCCAGTGGTTTTTTTTTTATTTTACCCTAACCGGTGTGGGagt                                                     >  2:414277/1‑144 (MQ=255)                                                                                                    ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATTAACAACCCCATTATTATTTTCAGTGGGTCGTTTTTTCCATCTTTTGGGTGGTTCTTTTTTACTTTCACCCTAACCGGTTCGGGGTTTCCATCCGtt                                           >  2:918137/1‑151 (MQ=255)                                                                                                       gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATGAAAACCCCCCTTATTATTTTCAGTGGGTTGTTTATTTCATCTTTCGAGGGCTTTTTTTTTTCTTTCCCCCTAAACGGGTCGGGGTTTGCATCCGtt                                           >  2:818521/1‑148 (MQ=255)                                                                                                         tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTacac                                       <  1:189314/150‑1 (MQ=255)                                                                                                             agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAC‑ATTAAAAACTCCCTTTTTTTTTTCCATGGTTTGTTTTTTTCCTCTTTTGTGTGCTTCTTTTTTTCTTTCACCCTAAACCGTTTGGGGTTCGCCTCCCTTtcac                                       >  1:1117524/1‑146 (MQ=255)                                                                                                             agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCg                                   <  2:666109/150‑1 (MQ=255)                                                                                                                   gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTa                            <  1:838979/151‑1 (MQ=255)                                                                                                                       ccttccttTTCAATAGTGGCGGTAATTTTTTTA‑ATTAACACCCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTGCTTTTTTTTTACTTCCCCCCTAACCGGGTCGGGGTTTCCATCCGTTACACTTCCCTCGGGaaaa                         >  2:885026/1‑148 (MQ=255)                                                                                                                       ccttccttTTCAATAGTGGCGGTAATTTTTTAC‑ATGAACACCCCCATTATTTTTTTCAGTGGGTCGGTTTTTTCCTCTTTTGGGGGCTTCTTTTTTACTTCCCCCATAAACGGTTCGGGGTTCGCATCCCTTTCACTTCCCTCGGGaa                           >  2:825716/1‑148 (MQ=255)                                                                                                                       ccttccttTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAata                         <  1:968623/150‑1 (MQ=255)                                                                                                                         ttccttTTCAATAGTGGCGGTAAATTTTTTA‑ATGAAAACCTCCCTTTTTATTTTCAGTGGTTTGTTTTTTTCATCTTTCTCGTGCTTTTTTTTTTCTTTCACCCTAAACCGTTTGGGGTTTGCCTCCCTTTCTCTTCCCTCGGTaaaa                         >  1:219451/1‑146 (MQ=255)                                                                                                                           ccttTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTg                    <  2:563273/151‑1 (MQ=255)                                                                                                                                       tGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc         >  1:979696/1‑150 (MQ=255)                                                                                                                                          cGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCg     >  1:684851/1‑151 (MQ=255)                                                                                                                                           ggTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt    >  2:489261/1‑151 (MQ=255)                                                                                                                                           ggTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt    <  2:1145309/151‑1 (MQ=255)                                                                                                                                            gTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt   <  1:987517/151‑1 (MQ=255)                                                                                                                                              aaTTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTa  <  2:91556/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAC‑ATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTA  >  NC_000913/1413588‑1413878

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