Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	49.4%	64.5 / 72.4	88	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (13/31);  new base T (43/0);  total (57/31)
Fisher's exact test for biased strand distribution p-value = 1.87e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413586‑1413885                                                                                                                                                       |                                                                                                                                                           aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTtat                                                                                                                                                           <  1:955087/151‑1 (MQ=255) aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTta                                                                                                                                                            <  2:17399/150‑1 (MQ=255) aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTta                                                                                                                                                            <  2:299175/150‑1 (MQ=255)    aGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAc                                                                                                                                                         >  2:925653/1‑150 (MQ=255)     gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑At                                                                                                                                                      <  2:802500/151‑1 (MQ=255)          aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGaaa                                                                                                                                                  >  2:139754/1‑149 (MQ=255)          aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAac                                                                                                                                                  <  1:523620/150‑1 (MQ=255)           cTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATTAaca‑c                                                                                                                                               >  1:1072310/1‑151 (MQ=255)           cTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAaca                                                                                                                                                 >  1:124795/1‑150 (MQ=255)            tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATTAACA‑cc                                                                                                                                              >  2:1446456/1‑151 (MQ=255)             ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTC‑CTGAACA‑ccc                                                                                                                                             >  2:412676/1‑150 (MQ=255)             ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAaca‑c                                                                                                                                               <  1:1144549/149‑1 (MQ=255)              tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTNTAN‑NNNNACA‑CCt                                                                                                                                             <  1:1137135/150‑1 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑cc                                                                                                                                          <  1:321208/149‑1 (MQ=255)                 gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCat                                                                                                                                        <  2:1106679/151‑1 (MQ=255)                     tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TCATTTTTTTG‑ATGAAAA‑Ccc‑cctatctt                                                                                                                                    >  1:120265/1‑145 (MQ=255)                        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCatt‑attat                                                                                                                                 <  1:623545/150‑1 (MQ=255)                        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCatt‑atta                                                                                                                                  <  2:506406/149‑1 (MQ=255)                         gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTATTTTTTAA‑ATTAAAA‑CCT‑CCCTT‑Atttttt                                                                                                                               >  1:208747/1‑151 (MQ=255)                            aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCa                                                                                                                            <  1:1293552/151‑1 (MQ=255)                               ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAT‑ATTAAAA‑CCT‑CCCTT‑TTTATTTTCAATg                                                                                                                         >  2:990423/1‑151 (MQ=255)                                    ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTATTTTTTAA‑ATAAACA‑CAT‑CCCTT‑TTTATTTTCCGTGGTTCg                                                                                                                    >  2:1380333/1‑151 (MQ=255)                                     tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTc                                                                                                                     <  2:1153064/149‑1 (MQ=255)                                      cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAACA‑CCT‑CCCTT‑ATTATTTTCAGTGGTTCGtt                                                                                                                  >  2:407511/1‑151 (MQ=255)                                        tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGttt                                                                                                                 <  1:719945/150‑1 (MQ=255)                                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCC‑CCATT‑ATTATTTTCAGTGGGTCGTTTAt                                                                                                               >  2:1232987/1‑151 (MQ=255)                                         ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTa                                                                                                                <  2:1266808/150‑1 (MQ=255)                                          tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTAt                                                                                                               <  2:117135/150‑1 (MQ=255)                                            tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCT‑CCATT‑ATTTTTTTCATTgggttttttttttt                                                                                                            >  1:1286723/1‑139 (MQ=255)                                            tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTc                                                                                                             <  1:741955/150‑1 (MQ=255)                                             aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTcc                                                                                                            <  2:1396639/150‑1 (MQ=255)                                              ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATTAACA‑CCC‑CCATT‑TTTTTTTTCCGTGGGTCGttttttttat                                                                                                          >  1:656525/1‑147 (MQ=255)                                                  agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCT‑CCCTT‑ATTTTTTTCCGTGGTTTGGTTTTTTCATCttt                                                                                                      >  2:174525/1‑151 (MQ=255)                                                  agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCttt                                                                                                      <  1:284792/151‑1 (MQ=255)                                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                                                                                                    <  1:929533/151‑1 (MQ=255)                                                     gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                                                                                                   <  1:260012/151‑1 (MQ=255)                                                         tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATTAAAA‑CCC‑CCCCT‑TTTATTTCCCGTGGGTCGTTTTTTcccttctttttttgg                                                                                               >  2:1037464/1‑138 (MQ=255)                                                          ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑TTTTTTTTCAGTGGTTCGTTTTTTTCATCTTTTTAGTGct                                                                                              >  1:1167906/1‑151 (MQ=255)                                                           tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAAAA‑ACT‑CCATT‑TTTTTTTTCAGTGGTTCGTTTTTTTCATCTTTTTAGTGctt                                                                                             >  1:630877/1‑151 (MQ=255)                                                            ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATAAAAA‑ACC‑CCCTT‑ATTATTTTCCGTGGTTCGTTTATTTCCTCTTTCGAGTGcttc                                                                                            >  2:274845/1‑151 (MQ=255)                                                            ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAT‑AATAAAA‑CCT‑CCCTT‑ATTATTTTTTGTGGGTCGGTTTTTCCTTCTTTCggattcttc                                                                                            >  2:494305/1‑143 (MQ=255)                                                              cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑ACC‑CCCTT‑TTTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGcttttt                                                                                          >  1:578789/1‑151 (MQ=255)                                                              cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑AACAAAA‑CCTCCCATA‑TTTTTTTCCAAGGGTT‑GTTTTTTTCAACTTTctgggggttttt                                                                                          >  2:128922/1‑140 (MQ=255)                                                               gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttt                                                                                         <  2:173677/151‑1 (MQ=255)                                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATTAAAA‑CCT‑CCCTT‑ATTATTTTCCGTGGGTCGTTTTTTTCCTCTTTTTAGTGCtttttttt                                                                                       >  2:182064/1‑151 (MQ=255)                                                                 aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                       <  1:159647/151‑1 (MQ=255)                                                                   tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TCATTTTTTTA‑ATGAAAA‑CCC‑CCCCT‑TTTTTTTTTTGTGGTTCGGTTATTTCCTCTTTTGAGTGCCTCttttttt                                                                                     >  1:1120081/1‑149 (MQ=255)                                                                    cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTAATTTTTAA‑ATGAACA‑CCT‑TCATT‑ATTATTTTCAGGGGTTTTTTTTTTTCACCTTTTCCGGGTTTCTTTTTCAc                                                                                    >  2:660427/1‑151 (MQ=255)                                                                     tGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑AAGAAAA‑ACC‑CCCTT‑TTTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTGAGTGGCTCttttttaat                                                                                   >  1:1226529/1‑147 (MQ=255)                                                                      gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TTATTTTTTTC‑ATGAACA‑CCT‑CCCTT‑TTTTTTTTCACAGGTGTGTTTTTTTCATCTTTCCAGTTCTTCTTTTTTACtt                                                                                  >  1:123877/1‑151 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATTAACA‑CCT‑CCCTT‑ATTATTTTCAGTGGTTCGTTTTTTTCCTCTTTCTAGGGCTTTTTTTTTACTTccacc                                                                             >  2:1369874/1‑151 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATTT‑AATAACA‑ACC‑CCCTT‑ATTATTTCCAGTGGGTCGTTTTTTTCTTCTTTTCGGGGGTTCTTTTTTACTTccacc                                                                             >  2:1211761/1‑151 (MQ=255)                                                                           cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCC‑CCCTT‑ATTATTTTCCGTGGGTTGTTTTTTTCATCTTTTGAGTGCTTTTTTTTTACTTccccc                                                                             >  2:751608/1‑151 (MQ=255)                                                                            gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATTAACA‑ACT‑CCCTT‑TTTTTTTTCAGTGGTTTGTTTTTTTCCTTTTTTGGGTGCCTCTTTTTTTCTTtcacc                                                                             >  1:1359334/1‑150 (MQ=255)                                                                               aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGATAATTTT‑TAA‑AATAACA‑CCT‑CCAAT‑TTTTTTTTCAGGGGTTGGTGTTTTATATTCTTTGTAAGCTTCTTTTTCCCTTcctccaca                                                                          >  2:317675/1‑148 (MQ=255)                                                                               aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑TTTATTTTCAGTGGGTTGTTTATTTCCTCTTTTGAGTGGTTTTTTTTTTCTTCCACCATa                                                                          >  2:1167004/1‑150 (MQ=255)                                                                                 tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑GAAATTTTTAC‑ATGAAAA‑CCC‑CCCTT‑ATTTTTTTCACTGGTTCGTTTTTTTCCTCTTTTGAGTGGTTTTTTTTTACTTTCACCCTaaac                                                                       >  2:773314/1‑149 (MQ=255)                                                                                  ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATTAACA‑CCC‑CCCCT‑ATTTTTTCCAGCGGGTTGTTTATTCCACCTTTTCAGTGCTTCTTTTTTACTTCCCCCAtcaacc                                                                      >  2:182430/1‑146 (MQ=255)                                                                                  ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAAcc                                                                       <  1:848891/150‑1 (MQ=255)                                                                                    aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCgg                                                                     <  1:1098802/150‑1 (MQ=255)                                                                                    aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGt                                                                    <  1:1118772/151‑1 (MQ=255)                                                                                     gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATTAAAA‑CCC‑CCCAT‑ATTATTTTCAATGGGTTGTTTTTTTCATTTTTTCAGGGCTTTTTTTTTTCTTCCCCCCCAAcccggt                                                                   >  1:1203874/1‑147 (MQ=255)                                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAACA‑CCT‑CCCTT‑ATTTTTTTCCGTGGGTTGTTTTTTTCATCTTTTGAGTGCTTCTTTTTTTCTTTCCCCCTAAACGGGTCgggg                                                              >  1:454672/1‑150 (MQ=255)                                                                                          tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAACA‑CCC‑CCATT‑ATTATTTTCAATGGTTTGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTACTTTCCCCATAAACGGGTTgggg                                                              >  2:814336/1‑150 (MQ=255)                                                                                            aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑AATAAAA‑CCC‑CCCTT‑ATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCACCCTAACCGGGTCGGGGtt                                                            >  2:1001154/1‑151 (MQ=255)                                                                                            aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAATTTATTA‑ATGAACA‑CCC‑CCCTT‑TTTATTTTCAGTGGTTTGTTTTTTTCCTCTTTTGAGGGCTTCTTTTTTACTTTCACCATAACCGGTTTgggggt                                                            >  1:332997/1‑148 (MQ=255)                                                                                               atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑AAGAACACCCC‑CCTTT‑TTTTTTTCCCGGGGTCCG‑TTATTCCATTTTTTGAGGGGTTCTTTTTTACTTCCACCATAAACGGCTCGGGGTTTGc                                                         >  2:1419609/1‑151 (MQ=255)                                                                                                 cTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCC‑CCATT‑TTTATTTTCAGTGGTTTGTTTTTTTCATCTTTTGAGTGCTTTTTTTTTTCTTCCACCCTAAACCGGTCGGGGTTCGCAt                                                       >  2:218038/1‑151 (MQ=255)                                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAAAA‑ACT‑CCCTTAATTTTTTTCAG‑GGTGCGTTTTTTTCTTCTTTTTAGTGCTTCTTTTTTACTTTCACCATCAACGCGTTGGGGTTCGCCTcc                                                     >  1:18429/1‑151 (MQ=255)                                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑TTTTTTTCCCGTGGGTTGTTTTTTTCCTCTTTTTAGTGGTTCTTTTTTTCTTTCCCCATAACCCGGTTGGGGTTTGCAttc                                                     >  1:72475/1‑149 (MQ=255)                                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGc                                                         <  2:1340501/147‑1 (MQ=255)                                                                                                    gTCCAGTNAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTANNTNNNGTGGTTCGTTTATTNNNNNNNTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCg                                                    <  2:105962/151‑1 (MQ=255)                                                                                                      ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAACA‑CCT‑CCCTT‑ATTTTTTTCCGTGGGTTGTTTTTTTCATTTTTTTAGTGCTTCTTTTTTTCTTTCACCATAACCGGTTTGGGGTTCGCCTCCGtt                                                  >  1:751654/1‑151 (MQ=255)                                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAC‑ATGAAAA‑CCC‑CCCTT‑ATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTGAGTGCTTCTTTTTTACTTCCCCCCTAAACGGGTCGGGGTTTCCATCCCTTac                                                >  2:1411708/1‑150 (MQ=255)                                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAT‑AAGGAAA‑CCC‑CCCTT‑TTTTTTTTCCCTGGTTTGTTTTTTTCCTCTTTTTTGGGCTTCTTTTTTTCTTCCCCCCTAAACCGTTCGGGGTTTGCCTCCCtttcc                                               >  1:736404/1‑148 (MQ=255)                                                                                                         gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTaca                                               <  1:1410352/151‑1 (MQ=255)                                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                           <  1:209045/151‑1 (MQ=255)                                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                           <  1:925649/151‑1 (MQ=255)                                                                                                                  aGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCg                                      <  1:407519/151‑1 (MQ=255)                                                                                                                   gTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAATTTTTAA‑ATGAACA‑CCT‑CCCTT‑ATTATTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCCCCCTAAACGGGTTGGGGTTCGCATCCCGTACACCTCGCTCgg                                     >  2:51893/1‑151 (MQ=255)                                                                                                                     gCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTTA‑ATGAACA‑CCC‑CCCCT‑TTTTTTTTCAGTGGTTTGTTTTTTTCATCTTTCGCGGGCTTCTTTTTTTCTTCCACCATAACCGGTTTCGGGTTCGCCTCCGTTACACCTCGCTCggg                                    >  2:700145/1‑149 (MQ=255)                                                                                                                      ccACCTTCCTTTTCAATAGTGGCGG‑GAATTTTTTACAATTAACC‑CCC‑CCCTT‑TTTTTTTCCA‑TGGTTCGTTTATTTCTTCTTTCGTGTGCTTCTTTTTTTCTTTCTCCATTACCGGGTTGGGGTTGTCTTCCGTTACACTTCGCTCGTTa                                   >  1:182619/1‑150 (MQ=255)                                                                                                                             ccttTTCAATAGTGGCGG‑TAATTTTTTAA‑ATGAAAA‑CCC‑CCCTT‑ATTTTTTTCAGTGGGTCGTTTATTTCATCTTTTGAGTGCTTCTTTTTTACTTCCCCCATAAACGGTTTGGGGTTTGCATCCGTTACACTTCGCTCGGGAATATCCTg                           >  2:303610/1‑151 (MQ=255)                                                                                                                                   cAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATgg                     <  1:1186173/151‑1 (MQ=255)                                                                                                                                       aGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                  <  2:124788/150‑1 (MQ=255)                                                                                                                                        gTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                <  2:1008895/151‑1 (MQ=255)                                                                                                                                        gTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGGTCGTTTATTTCATCTTTTGAGTGCTTCTTTTTTACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtttt                >  1:134264/1‑150 (MQ=255)                                                                                                                                           gCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc             >  2:1381507/1‑151 (MQ=255)                                                                                                                                              g‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt          >  1:1043222/1‑151 (MQ=255)                                                                                                                                                 aaTTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc        >  1:476681/1‑151 (MQ=255)                                                                                                                                                     ttATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt    >  1:1373104/1‑151 (MQ=255)                                                                                                                                                     ttATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt    >  1:386179/1‑151 (MQ=255)                                                                                                                                                      tataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt    <  1:139761/150‑1 (MQ=255)                                                                                                                                                      tataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   >  1:323660/1‑151 (MQ=255)                                                                                                                                                       ataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  1:1437566/151‑1 (MQ=255)                                                                                                                                                       ataC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  2:1202449/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                           ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTATAC‑ATGAACA‑CCT‑CCATT‑ATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413586‑1413885

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: ‑