Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,504	0	A	C	43.2%	124.2 / 32.6	83	T291P (ACC→CCC)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (36/10);  new base C (0/35);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 3.96e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGT  >  NC_000913/436355‑436642                                                                                                                                                        |                                                                                                                                            tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCAGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTa                                                                                                                                             <  2:581650/150‑1 (MQ=255)  ttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc                                                                                                                                           >  1:734305/1‑151 (MQ=255)  ttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc                                                                                                                                           >  1:915816/1‑151 (MQ=255)   tCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCt                                                                                                                                          >  1:648613/1‑151 (MQ=255)      gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCt                                                                                                                                          <  1:503596/148‑1 (MQ=255)      gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCt                                                                                                                                          <  1:753143/148‑1 (MQ=255)         gcgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCg                                                                                                                                    <  1:821920/151‑1 (MQ=255)           gcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCggg                                                                                                                                  >  1:192002/1‑151 (MQ=255)             tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                 >  2:363025/1‑150 (MQ=255)              ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGGTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGGTTACCTGTATCGGGCa                                                                                                                                >  2:24780/1‑150 (MQ=255)               tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCaa                                                                                                                               <  1:186248/150‑1 (MQ=255)                  cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑Ga                                                                                                                          <  2:265117/151‑1 (MQ=255)                   gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GAc                                                                                                                         >  2:101227/1‑151 (MQ=255)                    ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GAcc                                                                                                                        >  2:965142/1‑151 (MQ=255)                     cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GAccg                                                                                                                       >  1:797000/1‑151 (MQ=255)                           gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGat                                                                                                                  >  2:86385/1‑150 (MQ=255)                              cgggAACCGGCCCTGCGCCGGGCGCTTTCGGGCGG‑GAAGG‑TTTGGGGTTTGGTTTTCACTGTGCCGGAACTGAACTGTGGCGGGGTGGATGTGGCCCTCGGACACGTGGGCGTACCTTTTCCCTGAATCGGGCAAAT‑GACCGCCGATATc                                                                                                               <  2:146869/148‑1 (MQ=255)                              ccGGAACGGGCGCTGCGCGGGGCGCTCTCTGGCGGTGAAAA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGGGGCGCGCGGGATGTGGCTCTCGGACACTGGGGCGTCCCTTTTCCCTGTATCGGGAAAAT‑GACCGCCGATATCg                                                                                                              <  2:596410/151‑1 (MQ=255)                               cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCg                                                                                                              >  1:817070/1‑150 (MQ=255)                                  aaCGGGCGCTGCGCTGGGCGCTCTCGGGGGGTGAAGA‑TTACGAG‑TTGTGTTTCCCTGTGCGGGAACTGACCCGTGGCGCGTTGGTTGTGGCTCTCGGACCCGTGGGCGTACCTTTTCCTTGTATGGGGAAAAT‑GACCGCCGATATCGAAg                                                                                                           <  1:393733/150‑1 (MQ=255)                                   acgggCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAT‑TTGTTTTTCCCTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCTTACCTTTTCCCTGTATCGGGCAAAT‑GACCGCCGATATCGAAggg                                                                                                         <  1:220521/148‑1 (MQ=255)                                          ctgcgctgGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGtt                                                                                                  >  2:108905/1‑151 (MQ=255)                                            gcgcggGGCGCTCTCGGGCGGTGAAGA‑TTACGTG‑TTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGCCACGTGGGCGTCCCTTTTCCCTGTACCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGtttt                                                                                                <  2:707901/151‑1 (MQ=255)                                             cgctgGGGGCTCTCTGGGGGT‑AAGATTTACGAT‑TTGTTTTTCCGTGGGCCGAAACTGAACCGTGGCGGGCGGGATGTGGCTCCCGGCCCCCGGGGCGTACTTTTTCCTTGTAGCGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTg                                                                                                    <  2:890652/146‑1 (MQ=255)                                                tgGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                            >  1:263309/1‑151 (MQ=255)                                                 gggCGCTCTCGGGGGGTGAAGA‑TTAGGAG‑TTGTTTTTCCCTGTGCCGGAACTGACCGGGGGCGCGTGGGATGTGGCTCTCGGACCCCGGGGCGTCCCGTTTACCTGTATCGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                              <  1:38504/148‑1 (MQ=255)                                                 gggCGCTCTCGGGGGGTGAAGA‑TTACGAG‑TTGTTTTTCCCTGTGCGGGAACTAACCGGTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCTTTTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGt                                                                                                   <  1:871843/143‑1 (MQ=255)                                                 gggCGCTCTCGGGCGGTGAAGA‑TTACGAG‑TTTTGTTTCATTGTGCGGGAACTGACCGGTGGCGCGTTGGTTGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTTTTTCGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                             <  1:456393/149‑1 (MQ=255)                                                 gggCGCTCTCGGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACGGAACCGTGGCGCGCTGGATGTGGCTCTCGAACACGGGGGCGTCCCGTTTCCCTGTATCGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                            <  1:624616/150‑1 (MQ=255)                                                 gggCGCTCTCGGGCGGTAAATT‑TTACGAG‑TTTTGTTTCATTGTGCCGGAACTGACCCGGGGCGCGCTGGGGGTGGCTCCCGCACACCTGGGCGTCCCTTTTCCCTGTAGCGGAAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                              <  2:960361/148‑1 (MQ=255)                                                  ggCGCTCTCGGGGGGTGAAGA‑TTACGTG‑TTGTGTTTCACTGCGCGGGAACTGACCCGTGGCGCGGTGGATGTGGCTCTCGGACCCCGGGGCGTCCTTTTTCCCTTTATCGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                          <  2:715633/151‑1 (MQ=255)                                                     gcTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                       >  1:517913/1‑151 (MQ=255)                                                     gcTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                       >  1:597328/1‑151 (MQ=255)                                                            gggggTGAAGA‑TTACGAT‑TTGTTTTTCACTGGGCGGAAACTGACCCGGGGCGCGTTGGATGTGCCTCCCGGACCCGGGGGCGTCCTTTTTCCTTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                <  2:995752/148‑1 (MQ=255)                                                            gggggTGAAGA‑TTACGAG‑TTTTGTTTCCCGGGGCCGGACCTAACCCGTGGCGCGCTGGATGTGGCTTTGGAACCCCGGGGCGTCCCTTTTCCCTGTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                <  2:86159/148‑1 (MQ=255)                                                            gggggTGAAGA‑TTACGAG‑TTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGCCCCCGGGGCGTCCCTTTTCCTTGTATGGGCCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                       <  2:830930/141‑1 (MQ=255)                                                            ggCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                >  2:933815/1‑151 (MQ=255)                                                             gCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                               >  1:612689/1‑151 (MQ=255)                                                             gCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                               >  2:56075/1‑151 (MQ=255)                                                                   aaGA‑TTGCGAG‑TTGTTTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGTATGTGGCTCGCGGACACTTGGGCGTACCGTTTACCTGTTTCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                         <  2:168711/151‑1 (MQ=255)                                                                    aGA‑TTACGAG‑TTGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGCTGGTTGTGGCTCTCGGACCCCTGGGCGTCCCGTTTCCCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                         <  1:355235/150‑1 (MQ=255)                                                                        ttACGAG‑TTGTTTTTCACTGTGCCGGAACTGACCGGGGGCGCGCTGGTTGTGGCTCTCGGACCCCTGGGCGTACCTTTTCCCTTTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                       <  1:600336/149‑1 (MQ=255)                                                                          tacggtttGTTTTTCACTGTGCCGGAACTGACCGGTGGCGCGTTGGATGGGGCTCTCGGCCCCGGGGGCGTCCCGTTTGCCTGTAGCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                     <  2:594023/144‑1 (MQ=255)                                                                            ggg‑ttgttTTTCCCTGTGCGGAAACTAACCGGGGGCGCGCGGGATGTGCCTCTGGGCCCCTGGGGCGTCCCTTTTCCTTTTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                      <  1:369548/144‑1 (MQ=255)                                                                                ttGTGTTTCAGTGTGCAGAAACTAACCCGTGGCGGGTGGGATGGGGCCCTCGGACACGGGGGCGCCCCGTTTCCTTGTATCGGAAAAGT‑GGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                    <  1:80206/145‑1 (MQ=255)                                                                                ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                               >  2:130509/1‑150 (MQ=255)                                                                                ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                               >  1:890065/1‑150 (MQ=255)                                                                                ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                >  2:705048/1‑149 (MQ=255)                                                                                 tttttttCACTGGGCCGAAACTGACCCGGGGCGCGCGGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTTTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                              <  2:1061915/146‑1 (MQ=255)                                                                                   tttttAACTGTGCAGGATCTGCCCGGGGGCGCGTTGGTTGTGGCTCCGGGCCCCCGGGGCCTACGTTTTCCCTATATGGGTAAAAA‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                             <  2:1051375/147‑1 (MQ=255)                                                                                   tgtTTCCCTGTGCGGAAACAAACCCGGGGCGCGCGGGAGGTGGCTCGCGAACCCCTGGGCGTTCCTTTTACCTGTATCGGAAAAAA‑AACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                            <  2:866103/150‑1 (MQ=255)                                                                                    gtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                          >  1:160900/1‑151 (MQ=255)                                                                                     tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                          >  1:461387/1‑150 (MQ=255)                                                                                      ttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                        >  2:1027226/1‑151 (MQ=255)                                                                                            gggCCGGAACTGACCCGGGGCGCGCTGGATGGGGCTCTCGGACACTTGGGCGTACCGTTTACCTGTATCGGGAAAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                  <  1:60346/149‑1 (MQ=255)                                                                                              gCCGGAACTGAACCGTGGCGCGGTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATCGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                 <  1:363097/150‑1 (MQ=255)                                                                                               cggaaacggaccggTGGGGCGTGGGATGTGTCTCGCGTAACCCGGGGCGTACCTTTTCCCTGTCTGGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                               <  2:131841/138‑1 (MQ=255)                                                                                               ccGGAACTGAACCGGGGCGCGTTGGATGGGGCTCTCGGACCCCTGGGCTTACCTTTTCCCTGTTTCGGACAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                       <  1:915138/143‑1 (MQ=255)                                                                                               ccGGAACTGAACCGGGGCGCGCTGGTTGGGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                <  1:57048/150‑1 (MQ=255)                                                                                                 ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                             >  1:389055/1‑151 (MQ=255)                                                                                                   atttaacccGGGGCGCGCGGGTTGTGGCTCCGGGACCCCGGGGCGTCCTTTTTCCCTGTATGGGGAAAAT‑ATCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                           <  2:871052/144‑1 (MQ=255)                                                                                                   aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGCCACGTGGGCGTACCGTTTCCCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                           <  1:129609/151‑1 (MQ=255)                                                                                                    aCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                          >  2:365737/1‑151 (MQ=255)                                                                                                      tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                        >  2:679226/1‑151 (MQ=255)                                                                                                      tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                        >  2:187386/1‑151 (MQ=255)                                                                                                       gAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                       >  2:150995/1‑151 (MQ=255)                                                                                                          ccGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                     >  1:805414/1‑150 (MQ=255)                                                                                                              ggCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                 >  1:488848/1‑150 (MQ=255)                                                                                                               gcgggttGGATGTGGCTCTGGGACACGTGGGCGTCCTGTTTCCCTGTACGGGAAAAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                 <  2:35340/143‑1 (MQ=255)                                                                                                                   gcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                           >  1:327522/1‑151 (MQ=255)                                                                                                                          gtgGCTCGCGGCCCCCGGGTCGTCCCGTTTCCCTGTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                          <  1:634134/145‑1 (MQ=255)                                                                                                                                   ggACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg           >  1:421849/1‑151 (MQ=255)                                                                                                                                      cccttGGGCGTCCCGTTTCCTTGTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa         <  1:548206/146‑1 (MQ=255)                                                                                                                                      ccccggggCGTCCCTTTTCCCTGTATGGGAAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa          <  1:700627/144‑1 (MQ=255)                                                                                                                                       cccggggCGTCCTGTTTCCCGTGATGGGGAAAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                <  2:394655/138‑1 (MQ=255)                                                                                                                                        cgggggCGTCCCTTTTCCTTGTTTCGGAAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa       <  1:101245/147‑1 (MQ=255)                                                                                                                                          tGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg     >  2:501863/1‑150 (MQ=255)                                                                                                                                          tGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg     >  2:95164/1‑150 (MQ=255)                                                                                                                                          tGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg     >  1:103362/1‑150 (MQ=255)                                                                                                                                           gggggtcctttttACCTTTAGGGGGCAAGT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg     <  2:734184/139‑1 (MQ=255)                                                                                                                                           gggCTTCCCGTTTCCTTTTCTCTGACAAAT‑GTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAg   <  1:151284/151‑1 (MQ=255)                                                                                                                                            ggCTTACCTTTTCCCTGTGCGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAg   <  1:824896/150‑1 (MQ=255)                                                                                                                                            ggCGTCCCTTTTCCTTGTATGGGGAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt  <  1:1047022/151‑1 (MQ=255)                                                                                                                                             ggggcccGTTTTCCTTGGTCGGGGAAATGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt      <  2:888838/142‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                            TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGA‑TTACGAG‑TTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGT  >  NC_000913/436355‑436642

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