Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	48.5%	71.7 / 37.3	66	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (29/5);  new base A (0/32);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.33e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATC  >  NC_000913/436366‑436655                                                                                                                                                         |                                                                                                                                            cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGc                                                                                                                                             >  2:304099/1‑151 (MQ=255)  tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGc                                                                                                                                             >  2:1071431/1‑150 (MQ=255)   ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCa                                                                                                                                            >  2:810119/1‑150 (MQ=255)       cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATg                                                                                                                                        <  1:352449/150‑1 (MQ=255)       cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGa                                                                                                                                       <  1:963750/151‑1 (MQ=255)                gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGat                                                                                                                               >  1:1001943/1‑150 (MQ=255)                   ccGGAACAGGCGCTGCGCGGGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGAAGCTGAACCGTGGCGCGTGGGATGTGCCTCTCGGACACCTGGGCGTCCGG‑TTTACC‑TGTATCGGGAAAATGACCGCCGATATc                                                                                                                            <  2:551310/150‑1 (MQ=255)                       aaCAGGCGCTGCGCTGGGCGCTCTCGGGGGGTGAAGATTACGAGTTGTGTTTCACTGGGCGGAAACTGACCCGGGGCGCGTGGGATGTGGCTCTCGGACCCCTGGGGGTCCCT‑TTTACC‑TGTATCGGGAAAATGACCGCCGATATCGAAg                                                                                                                        <  2:159933/150‑1 (MQ=255)                              gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                >  2:610619/1‑151 (MQ=255)                                   gttgGGCGCTCTCGGGCGGTGAAGATTACGGTTTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGTTGGATGTGGCTCTCGGACCCCGGGGGGCCCCT‑TTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                               <  2:505073/145‑1 (MQ=255)                                   gctgGGCGCTCTCTGGCGGTGAAGATTCCGAGTTGTGTTTAACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCTGACCG‑TTTACC‑TGAATCGGGCAACGGACCGACGATATCGAAGGACTTTGctt                                                                                                              >  2:325057/1‑148 (MQ=255)                                     tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                         >  1:503618/1‑151 (MQ=255)                                              tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                >  1:163475/1‑151 (MQ=255)                                              tcggggggTAAAGTTTGCGTGTTGTGTTTCCCTGTGCGGGAACTGACCGGGGGCGGGGTGGGTGTGGCTCTC‑GACCCCTGGGGGTTCCC‑TTTTCCTTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                   <  2:571696/142‑1 (MQ=255)                                                 gggggTGAAGATTACGAGTTGTGTTTCCCTGTGCCGAAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCT‑TTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                   <  1:1115406/142‑1 (MQ=255)                                                 ggCGGTGAAGATTACGTTTTGTTTTTCCCTGTGCCGGACCTGACCTGTGGCGGGTGGGATGGGGTTCTCGGCCACCGGGGCGTCCCG‑TTTCCT‑TGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                 <  1:346326/147‑1 (MQ=255)                                                   cGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                           >  2:471542/1‑151 (MQ=255)                                                    ggggAAGATTCCGTGTTGTGTTCCCCTGGGCCGAAACTGAACGGTGGCGCGCTGGTTGTGGCTCGCGGACGCGGGGGCGTACGT‑TTTACT‑TTTCAGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                           <  2:152454/147‑1 (MQ=255)                                                    ggTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                          >  1:585782/1‑151 (MQ=255)                                                        aaGATTACGGTTTGTTTTTCCTGGTGCCGAATCTAACCCGTGGCGCGTGGGATGTGTCTCTCGGACACCGGGGCGTCCCT‑TTTCCT‑TGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                      <  2:369918/151‑1 (MQ=255)                                                         aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                      >  2:1139775/1‑150 (MQ=255)                                                          gATTACGGTTTTTGTTTCTCTGTGCCGAAATTGAACCGGGGCGCGTGGGATGTGGCTCTGGGACCCCGGGGCTTACTT‑TTTCCT‑TTTTTGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACa                                                                                     <  2:1143637/150‑1 (MQ=255)                                                              aCGAGTTGTTTTTCACTGGGCCGGAACGGAACCGTGGCGCGCTGGGTGTGGCTCTCGGACACCTGGGGTTCCCG‑TTTCCC‑TGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                <  1:808606/151‑1 (MQ=255)                                                                   ttGTTTTCCACTGGGCGGGAACTGACCGGTGGCGCGCTGGATGTGGCTCTCGGACCCGTGGGCGTACCG‑TTTCCC‑TTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                               <  1:848530/147‑1 (MQ=255)                                                                    tgtgTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                            >  1:405287/1‑149 (MQ=255)                                                                      tttttCACTGTGCCGGAACTGAACGGTGGCGGGCTGGATGTGGCTCTCGGACACCTGGGCGTCCCT‑TTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                         <  1:773828/148‑1 (MQ=255)                                                                      tgttccgcggtgCGGAAACTGACGCGGGGCGCGGTGAGTGTGTCTCTGGAACCCCTGGGGGTGCCG‑TTTCCC‑TATAGCGGTAAAGTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                         <  2:605307/141‑1 (MQ=255)                                                                      tgtTTCACTGTGCCGGAACTGAACGGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTCCT‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                         <  1:831272/150‑1 (MQ=255)                                                                      tgtTTCACGGTGCCGAAACTGACCCGTGGCGCGTTGGATGTGGCTCGGGGACCCGTGGGCGTCCCG‑TTTCCT‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                         <  2:1116229/150‑1 (MQ=255)                                                                         ttCACGGGGCCGAAACTGACCCGTGGCGCGTTGGTGGGGGCTCGCGGACCCCGGGGCGTACTT‑TTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa                                                                      <  2:443535/150‑1 (MQ=255)                                                                           cACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                   >  2:835650/1‑151 (MQ=255)                                                                            aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                  >  2:444281/1‑151 (MQ=255)                                                                               gtgCCGGAACTGAACCGGGGCGCGTTGGATGTGGCTCTGGGACCCGTGGGCGTACCG‑TTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                               <  1:252607/151‑1 (MQ=255)                                                                                tgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                               >  2:1012366/1‑150 (MQ=255)                                                                                  ccGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                            >  1:137816/1‑151 (MQ=255)                                                                                   cGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                           >  1:550217/1‑151 (MQ=255)                                                                                    ggAACTGAACCGTGGCGCGGTGGATGTGCCTCTCGGCCACCTGGGCGTACCG‑TTTCCC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                           <  1:1039106/150‑1 (MQ=255)                                                                                         tGACCCGGGGCGCGGGGGATGGGGTTCGGGGACCCGGGGGCGTCCCG‑TTTCCC‑TTTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                       <  2:896600/149‑1 (MQ=255)                                                                                         tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                     >  1:1119501/1‑151 (MQ=255)                                                                                           agcggTGGCGCGCTGGATGTGGCTCTCGAACGCTGGGGGGTACCG‑TTTCCT‑GGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                   <  1:283959/147‑1 (MQ=255)                                                                                           acccGGGGCGCTGTGGATGTGGCTCTCGGCCACCTGGGCGTCCCT‑TTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                     <  2:365372/147‑1 (MQ=255)                                                                                           acccGGGGCGCGTTGGATGTGGCTCCCGGACCCTGGGGCGTCCCT‑TTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                    <  1:778820/148‑1 (MQ=255)                                                                                            aCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                  >  1:866869/1‑151 (MQ=255)                                                                                             cggcggCGCGTTGTATGTGCCCGTCGGCCACTTGGGTGTACCG‑TTTTCT‑TGTCCGGGGAAAACAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                   <  2:802724/145‑1 (MQ=255)                                                                                             ccGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                  >  1:372469/1‑150 (MQ=255)                                                                                                 ggCGCGGGGGATGGGGCTCTGGGCCCCTGGGGCGTCCCT‑TTTCCT‑TTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                               <  1:620577/149‑1 (MQ=255)                                                                                                    gcgcTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                              >  1:824412/1‑147 (MQ=255)                                                                                                      gcTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                        >  1:335200/1‑151 (MQ=255)                                                                                                         ggttgtgGCTCTGGGAC‑CCTGGGGCTTCCGTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgc                                     <  1:1123942/148‑1 (MQ=255)                                                                                                         ggttgtgGCTCTCGGACCCGTGGGGGTACCT‑TTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                       <  1:625893/146‑1 (MQ=255)                                                                                                          gATGTGGCTCTGGGCCCCTTGGCGGTACTT‑TTTTCT‑TTTTGCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                              <  2:24760/141‑1 (MQ=255)                                                                                                          gATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGCAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                                    >  2:817967/1‑151 (MQ=255)                                                                                                                  tctcGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                            >  1:781893/1‑151 (MQ=255)                                                                                                                     cGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                         >  1:382053/1‑151 (MQ=255)                                                                                                                       gACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTTGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                        >  2:106574/1‑150 (MQ=255)                                                                                                                       gACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                        >  1:285825/1‑150 (MQ=255)                                                                                                                         cccgggggCGT‑CCG‑TTTTCCTTTTATCGGAAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa                      <  1:80285/145‑1 (MQ=255)                                                                                                                           cttGGGCGTCCCG‑TTTCCC‑TGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa                    <  1:403926/148‑1 (MQ=255)                                                                                                                              gggCGTCCGT‑TTTACT‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt                   <  2:503603/148‑1 (MQ=255)                                                                                                                              gggCGTACTT‑TTTCCC‑TTTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                        <  2:767218/143‑1 (MQ=255)                                                                                                                              gggCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                  >  1:609851/1‑149 (MQ=255)                                                                                                                                 cGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa              >  2:354447/1‑150 (MQ=255)                                                                                                                                  gTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt            >  2:591127/1‑151 (MQ=255)                                                                                                                                         tttACC‑TTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg       <  1:589161/150‑1 (MQ=255)                                                                                                                                             cc‑TGTGTCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg      <  1:731497/147‑1 (MQ=255)                                                                                                                                             cc‑TGTATGGGGAAATTAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt   <  1:887741/150‑1 (MQ=255)                                                                                                                                             cc‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc  >  2:150007/1‑151 (MQ=255)                                                                                                                                                tGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCa    <  1:1088055/147‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                            CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATC  >  NC_000913/436366‑436655

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