Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	44.1%	67.2 / 31.2	59	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (25/8);  new base A (0/26);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 5.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436368‑436666                                                                                                                                                        |                                                                                                                                                       ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCa                                                                                                                                                       <  1:21988/150‑1 (MQ=255)    tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATg                                                                                                                                                   >  2:916923/1‑151 (MQ=255)     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGa                                                                                                                                                  <  1:921853/151‑1 (MQ=255)        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGAccg                                                                                                                                               <  1:956697/151‑1 (MQ=255)        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGGGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGAc                                                                                                                                                 >  2:257148/1‑149 (MQ=255)               acccGGAACAGGCGCTGCGCTGGGCGCTCTCGGGCGGTGAAGTTTACGGTTTTTTTTTCCCTGGGCCGGAACTGAC‑CGGGGGCGGGGTGGATGTGGCTCTCGGCCCCCTGGGCTTCCC‑TTTTCCCTGTATGGG‑ACAAATGACCGCCGata                                                                                                                                         <  1:712519/148‑1 (MQ=255)                      aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAggg                                                                                                                                 >  2:875602/1‑151 (MQ=255)                        gggcgcggcgcGGGGCGCTCTCGGGGGGGGAAGTTTAGGGGTTGTTTTTCTCTGGGCGGAAACTGAC‑CGGGGGCGGGTGGGATGTGCCTCTCGCACCCCGGGGGGTCCC‑TTTTCCCTGTAGGGG‑AAAAATGACCGCCGATATCGAAGGGCt                                                                                                                               <  1:24252/150‑1 (MQ=255)                         ggCGCTGCGGTGGGCTTTTTCGGGCGGGGAAGATTACGTGTTGTTTTTCCCTGCCCGGGACCTAAC‑CCGTGGCGCGCGGGTTGTGCCTCTCGGACCCCGGGGTTTCCC‑TTTTCCCTGTATCGG‑AAAAAGGCCCGCCGATATCGAAGGGCtt                                                                                                                              <  1:241861/151‑1 (MQ=255)                         ggCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCtt                                                                                                                              >  1:936238/1‑151 (MQ=255)                            gctgcgcggGGCGCTCTCTGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGAA‑CCGGGGCGCGTTGGAGGTGGCTCTCGAACCCCTGGGCGTCCC‑GTTTCCCTGTATGGG‑GAAAATGACCGCCGATATCGAAGGGCttt                                                                                                                             <  2:233655/149‑1 (MQ=255)                               gcggggggCGCTCTGGGGCGGTGAAGTTTAGGAGTTGTTTTCCACTGTGCCGAAACTGAC‑CGGGGGCGCGTTGGATGTGGCTCTCGGACCCCGGGGCTTCCC‑GTTTCCTTTTATCGG‑CAAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                         <  2:193437/145‑1 (MQ=255)                               gcgctgGGCGCTCTGGGGGGGTGAAGATTACGGGTTTTTTTTCACTGTGCCGGAACTGAC‑CCGTGGCGCGTTGGATGTGGCTCGCGGACCCCGGGGCGTGCC‑TTTTCCCTGTAGGGG‑GAAAATGACCGCCGATATCGAAGGGc                                                                                                                                <  1:932334/143‑1 (MQ=255)                               gcgctgGCCCCTTTCGGGGGGTAAAGATTCCGATTTGGGTTTCGCGGCGCGGGATGTGAC‑CGGTGGCGCGTTGGATGGGCCCCTCGGACGCGGGGGCGGCCC‑TTTTACCTG‑ATCGGAGCAAAAGACCGCCGATATCGAAGGGCttt                                                                                                                             <  1:1050411/146‑1 (MQ=255)                                cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                        >  1:217721/1‑150 (MQ=255)                                   tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    >  2:356946/1‑151 (MQ=255)                                    gggCGTTTTCGGGGGGGGAAGATTAGGTGTTTTTTTTCCCGGGGCCGGACCGGAA‑CCGTGGCGCGTTGGATGTGGCTCCCGGACCCTGGGGGGTCCC‑TTTTCCTTGTATGGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                     <  2:657630/149‑1 (MQ=255)                                       cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                                 >  2:556835/1‑150 (MQ=255)                                         ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                              >  1:491399/1‑151 (MQ=255)                                          tcactGGCGGGGAAGATTACGAGTTGTGTTTCACCGTGCCGGAACTGAA‑CCGTGTCGCGCGGGATGTGGCTCTCGGGCACCTGGGCGTACC‑GTGTACCTGTATCGG‑GCAGATGACCGCCGATCTCGAATGGCGTTGTTTTATTCGTGa                                                                                                                >  1:165187/4‑148 (MQ=255)                                              tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                          >  2:785036/1‑150 (MQ=255)                                               gggggTGAAGTTTTGGTGTTGTGTTTCTCTGTGCGGAAATTAAG‑CCGGGGCGCGTGGGATGTGGCTCACGGTCCCGGGGGCGTCCT‑TTTTACCTTTACCGG‑TAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           <  2:57691/145‑1 (MQ=255)                                                 gggTGAAGATTGCGAGTTGTGTTTCACTGTGCCGGAACTGAC‑CGGGGGCGCGGTGGAGGTGGCTCTCGGACCCGTGGGCGTACC‑TTTTCCCTGTATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                      <  1:415409/150‑1 (MQ=255)                                                   gTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGt                                                                                                    >  1:85352/1‑151 (MQ=255)                                                           tcggtgTTGTGTTTCCTGGTGCC‑GAAATGAACCCGGGGCGGGTGGGATGGGCCTCGGGGACCCCGGGGCGTCCT‑TTTTCCCTGTAGGGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                             <  1:377491/145‑1 (MQ=255)                                                              gttttGTGTTTCTCTGGGCCGAAACTGAA‑CCGGGGCGCGTGGGATGTGGCTCTGGGCCCCTTGGGCGTCCC‑TTTTCCTTGTATGGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                         <  1:641036/148‑1 (MQ=255)                                                                 ttGTTTTTCCCTGGGCGGAAACTAAC‑CGGGGGCGGGTTGGTTGTGCCTCTCGGCCCCCGGGGCTTACT‑TTTTCCGTGTAGGGG‑GAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                           <  1:538430/146‑1 (MQ=255)                                                                 ttGTGTTCCCCTTTCCCGGACTGGAA‑CAGTGGCGCGCTGTTTTTCACTCCGGCCCACTTGGGCGTACG‑TTTTCCTTTTGTCGG‑CAAAATCCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                       <  2:940593/150‑1 (MQ=255)                                                                  tgttttccccTGGGCGGGAACTAAC‑CCGGGGCGCGTTGGAGGTGTCTCTCGGACCCGTGGGCATCCT‑TTTTCCGTTAATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  1:1004775/139‑1 (MQ=255)                                                                    tttttCACTGTGCCGAAACTGAA‑CCGTGGCGCGTTGGATGGGGCTCCCGAACCCCGGGGCGTCCC‑TTTTACCTGTATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                      <  2:401657/146‑1 (MQ=255)                                                                    ttttccccTGTGCCGAAACTGAC‑CCGTGGCGCGCGGGTTGTGGCTCTGGGCCCCCTGGGCTTCCC‑GTTTCCTTTTCTCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                                    <  1:279177/143‑1 (MQ=255)                                                                      tttCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                                   >  1:1071653/1‑149 (MQ=255)                                                                       ttCACTGTGCCGGAACTAAA‑CCGTGGCGCGCTGGTTGTGGCTCCCGGCCCCCGGGGCGTCCC‑GTTTCCCTGTATGGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa                                                                                 <  2:43804/150‑1 (MQ=255)                                                                         cACTGGGCCGAAACTGAC‑CCGGGGCGCGTTGGAGGTGGCTCCCGGACCCCGGGGCGTCCT‑TTTTACCTGTATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                               <  1:920778/150‑1 (MQ=255)                                                                             gggcgggAACTGAC‑CGGGGGCGCTCTGAGTGTGGCTCACGAACCCCGGGGGGTACC‑TTTTCCCTGTAGGGG‑GCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                                <  1:813636/140‑1 (MQ=255)                                                                                cgggAACTGAC‑CGGTGGCGCGTGGGTTGGGGCTCTCGGCCCCGGGGGCTTCCC‑GTTTCCTTGTATGGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                                <  2:134619/140‑1 (MQ=255)                                                                                ccGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                       >  1:1032240/1‑151 (MQ=255)                                                                                 cGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                      >  1:767201/1‑151 (MQ=255)                                                                                    accTGAA‑CGGTGGCGCGCTGGATGTGGCTCTCGGCCCCCTGGGCGTCCC‑TTTTACTTGTATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     <  1:616850/147‑1 (MQ=255)                                                                                    aaCTGAC‑CCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCT‑TTTTCCTTGTATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                      <  1:940621/148‑1 (MQ=255)                                                                                         aa‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGAGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                              >  2:822084/1‑151 (MQ=255)                                                                                            cgggggCGCGCTGGTTGTGGCTCTCGGACCCCGGGGCGTCCC‑TTTTCCTTGTGTGGG‑AAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                              <  2:612034/145‑1 (MQ=255)                                                                                            cgggggCGCGCTGGTTGGGGCTCTGGGACCCCTGGGGGTCCC‑TTTTCCCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                             <  2:254019/146‑1 (MQ=255)                                                                                              gTGGCGCGTTGGATGTGGCTCTGGGACACCTGGGCGTACC‑TTTTCCCTGTATCGG‑GAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                             <  1:125432/148‑1 (MQ=255)                                                                                                     gcTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt                                                     >  2:329213/1‑149 (MQ=255)                                                                                                      cTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  >  2:492387/1‑151 (MQ=255)                                                                                                        ggttgtgGCTCTGGGAC‑CCGGGGGGTTCCTTTTTCCCTGTTTGGG‑AAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  <  2:30612/146‑1 (MQ=255)                                                                                                                       acacCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  1:346168/1‑150 (MQ=255)                                                                                                                           cTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             >  2:693650/1‑151 (MQ=255)                                                                                                                           cTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             >  2:284462/1‑151 (MQ=255)                                                                                                                             ggggttccc‑GTTTCCCTGTTGGGG‑AAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAg                           <  2:131941/144‑1 (MQ=255)                                                                                                                              gggcgctctttttGCCTGTAGGGG‑GAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             <  2:688832/139‑1 (MQ=255)                                                                                                                                       gTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGTCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                  >  2:520029/1‑151 (MQ=255)                                                                                                                                             cTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt            >  1:777839/1‑151 (MQ=255)                                                                                                                                                   cGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc      >  2:517130/1‑151 (MQ=255)                                                                                                                                                    gg‑gAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt            <  2:171069/144‑1 (MQ=255)                                                                                                                                                    gg‑gAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt         <  2:760681/147‑1 (MQ=255)                                                                                                                                                     g‑gCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTg    >  2:695471/1‑151 (MQ=255)                                                                                                                                                       gCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGt   >  1:769087/1‑151 (MQ=255)                                                                                                                                                        cAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  2:745786/1‑151 (MQ=255)                                                                                                                                                        cAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  1:1067168/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACC‑GTTTACCTGTATCGG‑GCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436368‑436666

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