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

CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413588‑1413885                                                                                                                                                    |                                                                                                                                                           cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtata                                                                                                                                                          <  3:222196/150‑1 (MQ=255) cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                                                                                                         <  4:414256/151‑1 (MQ=255)  aGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGAATATTtat                                                                                                                                                           >  2:824610/1‑148 (MQ=255)      ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                                                                                                     <  1:89259/150‑1 (MQ=255)       tACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACa                                                                                                                                                        >  1:836339/1‑146 (MQ=255)         cTTTCTGCTCTCTGGACACGCCTGCTTCTTTTTTACCTGAGAGAATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGATTTTtatttatgcaga                                                                                                                                                  >  2:145933/1‑141 (MQ=255)               gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCa‑t                                                                                                                                          <  4:250035/151‑1 (MQ=255)               gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATTAAAACCTcc                                                                                                                                             >  2:478174/1‑149 (MQ=255)                ctctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAACACcccc                                                                                                                                             >  2:115900/1‑148 (MQ=255)                  ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAAATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCa‑tt                                                                                                                                         <  1:1025150/149‑1 (MQ=255)                   tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATGAAAACcccccctttt                                                                                                                                       >  1:68234/1‑145 (MQ=255)                         aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTcc                                                                                                                                <  1:126790/151‑1 (MQ=255)                         aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTcc                                                                                                                                <  1:1151020/151‑1 (MQ=255)                          aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCa                                                                                                                               <  4:332596/151‑1 (MQ=255)                                gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTATTAAAACCCCCC‑TTTTTTTTTTCAGT‑gg                                                                                                                          >  2:278760/1‑149 (MQ=255)                                gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTTATATGAAAACCCCCC‑CTTTTAtttttcgtgggt                                                                                                                         >  1:1200112/1‑142 (MQ=255)                                   tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTc                                                                                                                       <  2:1197066/149‑1 (MQ=255)                                   tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCg                                                                                                                      <  4:57118/150‑1 (MQ=255)                                     ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATTAAAACCCCCC‑TTTTTTTTTTCatggggttgttt                                                                                                                    >  2:983866/1‑140 (MQ=255)                                      tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGttt                                                                                                                   <  3:129608/150‑1 (MQ=255)                                      tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGttt                                                                                                                   <  3:81625/150‑1 (MQ=255)                                        tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTAt                                                                                                                 <  4:19677/150‑1 (MQ=255)                                            ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCAt                                                                                                            <  2:662518/151‑1 (MQ=255)                                               gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATc                                                                                                           <  1:151415/149‑1 (MQ=255)                                               gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCt                                                                                                          <  4:106032/150‑1 (MQ=255)                                                 gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATAAAAACCCCCC‑TTTTTTTTTTCAGT‑GGGTTGTTTTTTTCCTCttt                                                                                                        >  1:1262909/1‑150 (MQ=255)                                                  agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGCGGCGGTAAATTTTTACATGAAAAACTCCC‑CTTTTATTTCCACT‑GGTTGGTTTTTTTCATCTTTct                                                                                                      >  1:718227/1‑150 (MQ=255)                                                     aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAAAACCCCCC‑TTTTTTTTTTCAGT‑GGTTTGTTTTTTTCTTCttt                                                                                                        >  2:940028/1‑146 (MQ=255)                                                     aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAAACCCCC‑TTTTTTTTTTCCCT‑GGGTTGTTTTTTTCCTCtttttagt                                                                                                   >  1:764213/1‑146 (MQ=255)                                                     aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTAAAATGACACCCCCA‑ATTTTTTTTCTAAT‑GGGTCGttttttttctcttttta                                                                                                     >  2:948378/1‑138 (MQ=255)                                                     aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCCTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCCCCC‑CTTTTATTTTCAGT‑GGGTTGTTTTTTTTCTCTTTTGAg                                                                                                    >  4:19494/1‑150 (MQ=255)                                                            cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGcttct                                                                                             <  1:249695/150‑1 (MQ=255)                                                                ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAACACCCCCC‑CTTTTTTTTTTCGT‑GGGTCTTTTTTTTCCTCTTTTTAGTGCttttttttt                                                                                        >  2:346928/1‑150 (MQ=255)                                                                 tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc                                                                                        <  1:81116/150‑1 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑tt                                                                                   <  2:772844/151‑1 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑t                                                                                    <  4:281029/150‑1 (MQ=255)                                                                    gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑t                                                                                    <  1:468430/150‑1 (MQ=255)                                                                       ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAACACCCCCC‑TTTTTTTTTTTAGT‑GGGTTGTTTTTTTCCTTTTTTTTGTGCttttttttttctt                                                                                    >  2:759850/1‑143 (MQ=255)                                                                         cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTccac                                                                               <  3:299156/150‑1 (MQ=255)                                                                          gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCCCCC‑TTTTTTTTTTCCGT‑GGGTTGTTTTTTTCCTCTTTTGAGTGGttttttttttctttcccc                                                                               >  4:368071/1‑140 (MQ=255)                                                                          gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGTACACCCCCC‑ATCTTTTTTTCAAT‑GGTTCGTTTTTTTTCTCTTTTTTGTGCTTCtttttttctttcccc                                                                               >  2:486677/1‑140 (MQ=255)                                                                          gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTccacca                                                                             <  1:983321/151‑1 (MQ=255)                                                                           ttAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAAACCCCC‑TTTTTTTTTTCATT‑GGGTCGGTTTTTTCTTCTTTTTAGTGCTTTTTTTTTTC‑TTtcac                                                                               >  2:1216254/1‑148 (MQ=255)                                                                             aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCCCCC‑CTTTTTTTTTCAGT‑GGGTCGTTTTTTTCATCTTTTTAGTGGTTTTTTTTTTC‑TTccccc                                                                              >  1:1081783/1‑147 (MQ=255)                                                                             aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATaa                                                                          <  1:660511/151‑1 (MQ=255)                                                                             aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAAATTTTTACAATAAAACCCCCACTTTTTTTTTTCGGG‑GGTT‑GTTTTTTTCATTTTTTTTGTGCTTCttttttattttctcc                                                                               >  1:889727/1‑137 (MQ=255)                                                                                   gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGtt                                                                   <  1:765691/151‑1 (MQ=255)                                                                                    aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTc                                                                  <  4:284043/151‑1 (MQ=255)                                                                                         gAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGa                                                              <  1:867576/150‑1 (MQ=255)                                                                                         gAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGa                                                              <  1:617101/150‑1 (MQ=255)                                                                                          aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGAtt                                                            <  2:399894/151‑1 (MQ=255)                                                                                              tCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATTAAAACCCCCC‑CTATTTTTTTTATG‑GGGTTGTTTTTTTCATTTTTTTCGTGCTTTTTTTTTTC‑TTCCACCATAACCCGGTTgggggt                                                             >  2:246428/1‑145 (MQ=255)                                                                                               cTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCCCCC‑TTTTTTTTTTCAGT‑GGGTCGTTTTTTTCATCTTTTGTGTGCTTTTTTTTTTC‑TTTCCCCCTAA‑CACGGTCGGGGTTCGCa                                                        >  1:297715/1‑150 (MQ=255)                                                                                                 tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATCAAAACCTCCC‑TTATTTTTTTCATT‑TGGTCGTTTTTTTCCTTTTTTCTGTGCTTCTTTTTTTC‑TTTCACCCTAA‑CCGGTTCcggggttgcaccc                                                     >  2:194431/1‑141 (MQ=255)                                                                                                 tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCa                                                        <  2:331020/148‑1 (MQ=255)                                                                                                  gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAACACCCCCC‑CTATTTTTTTCCAT‑GGGTTGTTTTTTCCATCTTTTGTGTGGTTTTTTTTTTC‑TTTCCCCCTAA‑ACCGTTCgggggttgccttc                                                     >  1:142357/1‑140 (MQ=255)                                                                                                   tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAAGAAAACCCCCA‑ATTTTATTTTCACT‑GGGTTGTTTTTTTCCTCTTTCCCGTGCTTCTTTTTTTC‑TTCCACCATAA‑CCCGTTCGGGGTTCGCATcc                                                     >  1:820664/1‑149 (MQ=255)                                                                                                    ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGt                                                   <  2:479534/150‑1 (MQ=255)                                                                                                       gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCAATTTTTTAAAATGACACCGCCA‑GTTTTTTTTTCAGG‑GGGTTGTTTTTTTCCTCTTTTGGGTGCCTCTTTTTTTC‑TTTCACCATAACCCGGTTTGGG‑TTTTCATCCGtt                                                  >  3:224552/1‑148 (MQ=255)                                                                                                         tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCCCCC‑CTTTTATTTTCAGTGGGTTGGTTTTTTCC‑TCTTTTTGGTGGTTTTTTTTTTCTTTCCCCCAAAA‑CCGTTTGGGGttttcctcccgt                                                   >  1:1252325/1‑137 (MQ=255)                                                                                                          agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACt                                             <  2:646982/150‑1 (MQ=255)                                                                                                              ggAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCt                                        <  2:25860/151‑1 (MQ=255)                                                                                                               gAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTc                                       <  1:1087041/151‑1 (MQ=255)                                                                                                                 gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg                                     <  2:1205750/151‑1 (MQ=255)                                                                                                                      aCCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAata                                <  1:598448/151‑1 (MQ=255)                                                                                                                       ccttccttTTCAATAGTGGCGGTAATTTTTTTAAATAACACCCCCC‑TTTTTTTTTTCAGT‑GGGTCGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTTC‑TTTCACCATAA‑CCCGTTCGGGGTTTGCCTCCCTTACACTTCCCTCGGTaa                                  >  4:192071/1‑148 (MQ=255)                                                                                                                        cttccttTTCAATAGTGGCGGTAATTTTATAAATTAACACCTCCC‑CTTTTTTTTTCATT‑GGGTCGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTAC‑TTCCACCATAA‑ACGGGTTGGGGTTTGCATCCGTTACACTTCTCTCGGTaa                                  >  2:93539/1‑147 (MQ=255)                                                                                                                          tccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATc                              <  2:563703/149‑1 (MQ=255)                                                                                                                                tCAATAGTGGCGGTAATTTTTTAAATTAACAACCCCC‑TTTTTTTTTTCAAT‑GGGTTGTTTTTTTCATTTTTTGAGTGCTTTTTTTTTTC‑TTTCACCACAA‑ACCGTTCGGGGTTTCCCTTCCTTTCACTTCTCTCGGTGatattcctaaa                        >  2:326681/1‑142 (MQ=255)                                                                                                                                tCAATAGTGGCGGTAATTTTTTAAATTAAAACCCCCC‑CTTTTTTTTTCCCT‑GGGTTGTTTATTTCATCTTTTTAGTGGTTCTTTTTTTT‑TTTCCCCATAA‑CCGGGTCGGGGTTTGCCTCCCTTTCACTTTCCTCCGCAataatcctaa                         >  2:967221/1‑141 (MQ=255)                                                                                                                                 cAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATgg                     <  1:418781/151‑1 (MQ=255)                                                                                                                                      gTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTc                <  1:1060882/151‑1 (MQ=255)                                                                                                                                           ggTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt           >  1:817221/1‑151 (MQ=255)                                                                                                                                             tAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt          >  2:199326/1‑150 (MQ=255)                                                                                                                                              aaTTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc        >  2:623415/1‑151 (MQ=255)                                                                                                                                                   tataCATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt    <  2:462971/150‑1 (MQ=255)                                                                                                                                                   tataCATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   >  1:872554/1‑151 (MQ=255)                                                                                                                                                    ataCATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  >  2:387610/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                           CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCA‑TTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413588‑1413885

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