Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	42.2%	95.6 / 32.4	65	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (33/4);  new base A (0/27);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 3.72e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436370‑436663                                                                                                                                                    |                                                                                                                                                    ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAAt                                                                                                                                                 >  2:1132279/1‑151 (MQ=255)   cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAAt                                                                                                                                                 <  2:489484/149‑1 (MQ=255)    gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAc                                                                                                                                              <  2:363498/151‑1 (MQ=255)       aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccg                                                                                                                                            >  2:302300/1‑150 (MQ=255)         gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgccg                                                                                                                                         >  1:272490/1‑151 (MQ=255)           tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGat                                                                                                                                       >  1:257926/1‑151 (MQ=255)             aaCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGatat                                                                                                                                     >  2:598190/1‑151 (MQ=255)                 ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                  >  2:782791/1‑150 (MQ=255)                   aaCAGGCGCTGCGCGGGGCGCTCTCTGGCGGTGAAGATTACGGGTTGTTTTTCACTGTGCCGGAACTAAACCGGGGCGCGTTGGATGTGGCTCGCGGACCCGGGGGCGTACCGTTTCCC‑TTTATCGGGAAAATGACCGCCGATATCGAAgg                                                                                                                               <  2:316222/151‑1 (MQ=255)                    aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                              >  1:105599/1‑151 (MQ=255)                      aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                            >  1:851559/1‑151 (MQ=255)                       ggcgcggcgcGGGGCGCTCTCGGGCGGTGAAAATTACGGGTTGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGTTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCt                                                                                                                            <  1:243408/150‑1 (MQ=255)                          gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                        >  2:960863/1‑151 (MQ=255)                               gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                   >  2:1062657/1‑151 (MQ=255)                                 tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                  >  2:1092455/1‑150 (MQ=255)                                   ggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                               >  1:1039258/1‑151 (MQ=255)                                       ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGCTATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                           >  2:268177/1‑151 (MQ=255)                                          tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                        >  1:1164/1‑151 (MQ=255)                                             gggggTGAAGTTTACGATTTGTTTTTCCCTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGCCCCTCGGACCCCTGGGCGTCCCTTTTCCT‑TTTATCGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                           <  2:636454/142‑1 (MQ=255)                                             ggCGGTGAAGTTTACGGGTTGTGTTTCCCGGGGCGGAAACTGACCCGGGGCGCGTGGGAGGGGGCTCCCGGCCCCCTGGG‑GTTCCTTTTTCCTTTTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa                                                                                                             <  2:27108/143‑1 (MQ=255)                                             ggCGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCCCGGCCCCCTGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                          <  1:509365/146‑1 (MQ=255)                                             ggCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACGGTGGCGCGCTGGTTGGGGCTCTCGGACCCCGGGGCGTCCCGTTTACT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                              <  1:948843/142‑1 (MQ=255)                                              gCGGTGAAGATTACGGGTTGTTTTTCTCTGGGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGGGTCCCGTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                                    <  1:30601/151‑1 (MQ=255)                                                     aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                              >  2:112167/1‑150 (MQ=255)                                                      gtttACGGTTTTTGTTTCCCTGTGCCGAAACTGACCGGGGGCGCGCTGGATGTGCCTCGCGAACCCTTGGGCGTCCCTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACa                                                                                             <  1:822011/148‑1 (MQ=255)                                                        ttACGAGTTGTGTTTCACTGTGCGGGAACTGACCGGTGGCGCGTTGGATGTGGCTCTCGGACACCTGGGCGTCCCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACa                                                                                             <  1:156142/148‑1 (MQ=255)                                                          acgttttttttttCCCTGTGCCGAAACTGACCCGTGGCGCGCTGGAGGTGCCTCTCGGACCCCGGGGCGTACTTTTTCCT‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                         <  2:20468/140‑1 (MQ=255)                                                          acggtttGTTTTTCCCTGTGCCGAAACTAACCGGTGGCGCGCTGGATGTGGCTCTCGGCCCCCTGGGCGTACCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                        <  1:1107229/146‑1 (MQ=255)                                                           cGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                                       >  2:935108/1‑151 (MQ=255)                                                            gatttGTGTTTCCCTTTGCCGGAGCTGATCCGGCGCGCGTTGGGTGTCGCTCGAGGACACCGGGGCTTACCGTTTGCT‑TTTCTCGGGTAATTTACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      <  1:163536/148‑1 (MQ=255)                                                             aGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                     >  1:367406/1‑151 (MQ=255)                                                              gTTGTGTTTCACTGGCCCGGACCTGACCGGGGGCGCGCTGGTGGTGCCTCTCGGACCCTGGGGCGTCCCGTTTCCC‑TGTCTCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                         <  2:785439/146‑1 (MQ=255)                                                               ttttttttCCCTGTGCCGGAACTAAACGGTGGCGCGCTGGTTGTGGCTCTCGGACCCGTGGGCGTACCGTTTCCT‑TTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                     <  2:256953/144‑1 (MQ=255)                                                                           gtgCGGGAACTGACCCGGGGCGCGTTGGATGGGGCTCTCGGACCCCTGGGCGTCCCTTTTACC‑TGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                             <  2:1127462/145‑1 (MQ=255)                                                                              ccGGAACTGACCCGTGGCGCGTTGGTTGTGCCTCCCGGACCCCGGGCCGTCCCTTTTACT‑TTTTTGGGAAAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                    <  2:275479/151‑1 (MQ=255)                                                                                      gAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                            >  2:224701/1‑151 (MQ=255)                                                                                       acccGGGGCGCGCTGGTGGTGCCTCTCGGACCCCGGGGCTTCCCGTTTCCC‑TGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                            <  2:437923/148‑1 (MQ=255)                                                                                          cGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                        >  1:762370/1‑151 (MQ=255)                                                                                          cGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                        >  2:885275/1‑151 (MQ=255)                                                                                           gTGGCGCGGGGGTGGGGGCTCTCGGACCCCTGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                        <  1:508976/150‑1 (MQ=255)                                                                                             ggCGCGTTGGTTGTGGCTCCCGACCCCCGGGGCGTCCCTTTTCCC‑TTTTGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                        <  2:1164985/148‑1 (MQ=255)                                                                                                     ggATGGGGCTCTCGGACCCCGGGGGGCCCCTTTTCCG‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                               <  1:515949/149‑1 (MQ=255)                                                                                                             ctctCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt                                      >  1:814252/1‑150 (MQ=255)                                                                                                               ctcggccccgggggCGTACTTTTTCCT‑TTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                       <  2:790604/136‑1 (MQ=255)                                                                                                               cccgaacacCGGGGCGTCCCTTTTTCC‑TGTATCGCGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                       <  1:800980/142‑1 (MQ=255)                                                                                                               cccGGACCCCGGGGCGTCCCGTTTCCT‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                     <  1:152337/147‑1 (MQ=255)                                                                                                                tcGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTTGCGAAAAGTCGcc                                  >  1:676132/1‑151 (MQ=255)                                                                                                                tcGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc                                  >  1:384112/1‑151 (MQ=255)                                                                                                                     cccttgggcttcctttttCCC‑TTTCTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa                              <  2:1064267/135‑1 (MQ=255)                                                                                                                      cccgggacgacctgTTTCCC‑TGTTGGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAg                             <  2:233231/137‑1 (MQ=255)                                                                                                                      acCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                               >  1:595920/1‑148 (MQ=255)                                                                                                                          gggggtcctttttCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa                         <  1:971643/140‑1 (MQ=255)                                                                                                                             cGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa                      >  2:751944/1‑150 (MQ=255)                                                                                                                              gTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                    >  1:37794/1‑151 (MQ=255)                                                                                                                                tcccttttccttttATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTaa                     <  1:933268/136‑1 (MQ=255)                                                                                                                                 ccGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                 >  2:72250/1‑151 (MQ=255)                                                                                                                                  cgattACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                 >  2:636070/4‑150 (MQ=255)                                                                                                                                     tttccttGTATGGGGAAAAGACCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                 <  1:860218/142‑1 (MQ=255)                                                                                                                                     ttACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc             >  1:979521/1‑151 (MQ=255)                                                                                                                                      tACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc             >  2:596138/1‑150 (MQ=255)                                                                                                                                         c‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt         >  1:508884/1‑151 (MQ=255)                                                                                                                                            gTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTAc       >  2:283808/1‑151 (MQ=255)                                                                                                                                               tCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTg    >  1:964498/1‑151 (MQ=255)                                                                                                                                                 gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt     <  2:105590/148‑1 (MQ=255)                                                                                                                                                 gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt  <  1:85308/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436370‑436663

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