Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	44.6%	94.4 / 39.1	65	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (31/5);  new base A (0/29);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.72e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436368‑436663                                                                                                                                                     |                                                                                                                                                    ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCa                                                                                                                                                    <  1:498664/150‑1 (MQ=255) ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCa                                                                                                                                                    >  2:479291/1‑150 (MQ=255) ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCa                                                                                                                                                    <  1:244351/150‑1 (MQ=255)     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGa                                                                                                                                               <  2:594288/151‑1 (MQ=255)      gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAc                                                                                                                                              >  2:817788/1‑151 (MQ=255)       ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                                             >  2:1045044/1‑151 (MQ=255)        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                                             >  2:264192/1‑150 (MQ=255)             tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                                       >  2:980316/1‑151 (MQ=255)             tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGa                                                                                                                                        <  1:62613/150‑1 (MQ=255)              gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                                       >  1:1017460/1‑150 (MQ=255)                   ggAACGGGCGCTGCGGTGGGCGCTCTCGGGCGGGGAAGATTAGGAGTTGTTTTTCCCGGTGCCGGAACTGACCCGGGGCGCGCTGGGGGGGGCCCCCGGACCCCGGGGCTTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGaa                                                                                                                                 <  2:757977/151‑1 (MQ=255)                   ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGatat                                                                                                                                     >  2:578039/1‑147 (MQ=255)                        aGGGGCTGCGCTGGGCGCTCTCTGCCGGTCAAGCTTACGAGTCGTGTTTCCCTGTTCCGGAACTGAACTGTGGCGCTGTGGTTTTTGCTCTCGGCAACCCGGCCGTACCGCCTACCTCTATTGGGCCTATGACCGTCGCCATCgccgg                                                                                                                               >  2:1187558/1‑144 (MQ=255)                         ggcgcggcgcGGGGCGCTCTCGGGCGGTGAAGATTAGGATTTTTTTTTCGCGGTGCGGGAACTGAACCGGGGCGGGTTGGATGTGGCTCCCGGACCCCGGGGCGTCCCGTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                           <  1:463424/151‑1 (MQ=255)                             cggcgctgGGCGCTCTCGGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCttt                                                                                                                          <  1:76409/146‑1 (MQ=255)                               gcgctgGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTTTTTCACTGTGTCGAAACTGAACGGTGGCGCGGTGGGTGGGGCTCTCGGACCCTTGGGCGTCCCGTTTCCCTGTATGGGCAAAAGGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                       <  1:773332/149‑1 (MQ=255)                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                   >  1:39078/1‑151 (MQ=255)                                    gggCGCTCTCTGGCGGTGAAGATTACGTGTTGTGTTTCCCTGTGCGGGAACTGAACGGTGGCGCGGGGGATGGGGCTCTGGGCCACTGGGGCGTCCCGTTTCCTTGTAGGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTa                                                                                                                    <  2:813093/147‑1 (MQ=255)                                    gggCGCTCTCGGGGGGGGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGACCGGGGGCGCGTTGGATGTGGCCCCCGGACACTGGGGCGTTCTTTTTTCCTTTGTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                   <  2:85435/148‑1 (MQ=255)                                         ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                           >  1:473407/1‑151 (MQ=255)                                           ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                          >  2:275910/1‑150 (MQ=255)                                              ggggggggAAGATTAGGGTTGGTGTTCCACTGTGCCGGATTTGACCCGTGGGGCGTTGGATGTGGCGCCCGGCCTCCGGGGTGGCCCTTTTCCCTTTAGCGGGAAATGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                                            <  2:62807/138‑1 (MQ=255)                                               gggggTAAAGATTAGGAGTTGTTTTTCCCTGTGCCGAAACTGACCGGGGGCCCGCTGGATGTGGCTCTCGGACACCGGGGGGTCCCGTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                         <  1:369179/144‑1 (MQ=255)                                               ggCGGTGAAGATTACGTGTTGTGTTTCACTGGGCCGAAACTGACCCGGGGCGCGCGGGATGTGGCTCTCGGACCCGGGGGCGTCCCTTTTCCCTGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                        <  2:636951/148‑1 (MQ=255)                                                  gggaaaGTTTCCGAGTTGTTTTTCACGGTGCCGGAACTGAGCCGTGGCGGGGGGGATGTGGCTCTCGGACCCGTGGGGGTCCCTTTTCCCTGTAGCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                   <  1:143366/146‑1 (MQ=255)                                                  ggTGAAGATTACGGTTTGTTTTCCACGGTGCGGGAACTGACCCGTGGCGCGTTGGATGTGGCTCTCGGACCCCGGGGCGTACTTTTTCCTGGTAGCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                  <  2:878688/151‑1 (MQ=255)                                                  ggTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                  >  1:20744/1‑151 (MQ=255)                                                      aaGATTACGGGTTGTGTTTCATTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGCCTCTCGGACCCGTGGGGTTCCGGTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                <  1:809277/149‑1 (MQ=255)                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                              >  2:298445/1‑150 (MQ=255)                                                        gATTACGAGCTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                            >  2:345830/1‑151 (MQ=255)                                                          ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                           >  1:407290/1‑150 (MQ=255)                                                                gTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                    >  2:231772/1‑151 (MQ=255)                                                                 ttGTTTTTCTCTGTGCCGGAACTGACCGGGGGCGCGTTGGAGGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                     <  1:514783/149‑1 (MQ=255)                                                                 ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGTGGGTTGTGGCTCTCGGACCCGTGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                         <  1:980323/145‑1 (MQ=255)                                                                 ttGTGTTTCACTGCGCCGGAACGGACCGGTGGCGCGCTGGTTGTGGCTCTCGGACCCGTGGGCGTCCCTTTTACTTGTCTGGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                    <  2:509765/150‑1 (MQ=255)                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                               >  1:1099241/1‑150 (MQ=255)                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                               >  1:47750/1‑150 (MQ=255)                                                                       ttCTCTGTGCCGGAACTAACCCGGGGCGCGGTGGAGGTGGCTCGCGGACCCCGGGGCGTACCGTTTACCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa                                                                              <  1:491674/150‑1 (MQ=255)                                                                         ccttgtgCCGAAACTAAACCGGGGCGCGCTGGATGGGGCTCTCGGACCCGTGGGGGTCCTTTTTCCTTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                           <  1:314429/148‑1 (MQ=255)                                                                           cTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                          >  2:59687/1‑150 (MQ=255)                                                                             gtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                       >  1:721494/1‑151 (MQ=255)                                                                              tgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                       >  2:819606/1‑150 (MQ=255)                                                                              tgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                       >  1:447853/1‑150 (MQ=255)                                                                                ccGGAACTGACCGGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCGTTTACCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                      <  1:979791/149‑1 (MQ=255)                                                                                ccGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                        <  2:696769/147‑1 (MQ=255)                                                                                 cGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                   >  2:898775/1‑151 (MQ=255)                                                                                  ggAACTGACCGGGGGCGCGCGGGATGTGGCTCCCGGACCCCGGGGCGTCCTTTTTCCCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                    <  2:397569/149‑1 (MQ=255)                                                                                  ggAACTGAACGGTGGCGCGCTGGATGTGGCTCTCGGACCCGGGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                   <  1:400603/150‑1 (MQ=255)                                                                                  ggAACTGAACCGGGGCGCGCTGGTTGTGGCTCTCGGACCCGGGGGCGTCCCTTTTACCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                  <  1:583381/151‑1 (MQ=255)                                                                                    aaCTGACCCGGGGCGCGGTGGAGGTGGCTCCCGGCCACGGGGGCGTCCTTTTTACCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg                                                                 <  2:229286/150‑1 (MQ=255)                                                                                     aCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                               >  1:336199/1‑151 (MQ=255)                                                                                         accgggggCGCGTTGGGTGGGGCTCTCGGACACTGGGGCGTCCTTTTTCCTTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                            <  2:610852/144‑1 (MQ=255)                                                                                          aCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                          >  1:653656/1‑151 (MQ=255)                                                                                            cGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                        >  2:969097/1‑151 (MQ=255)                                                                                             ggggCGGGGTGGATGTGGCTCCCGGACCCCTGGGCGTCCTTTTTCCTTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                        <  2:427710/148‑1 (MQ=255)                                                                                                   cgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt                                                  >  1:558573/1‑150 (MQ=255)                                                                                                       ggATGTGGCTCCAGACCACGGGGGCGTCCCTTTTACCTGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt                                                  <  2:926208/146‑1 (MQ=255)                                                                                                                    ggACCCCGGGGCGTTCTGTTTCCTTGTATGGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt                                      <  2:1039791/145‑1 (MQ=255)                                                                                                                        atcggggCCTACTGTTTTCTTATATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAg                             <  2:1130293/146‑1 (MQ=255)                                                                                                                        acCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                               >  1:1132581/1‑148 (MQ=255)                                                                                                                         ccTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                               >  1:1051387/1‑147 (MQ=255)                                                                                                                                      tttccttgtcggggacAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  2:1077424/135‑1 (MQ=255)                                                                                                                                       ttcccTGTTGCGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg              <  2:1153451/147‑1 (MQ=255)                                                                                                                                           cTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt         >  1:886636/1‑151 (MQ=255)                                                                                                                                                 cGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc   >  2:295429/1‑151 (MQ=255)                                                                                                                                                  gggAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt  <  2:39078/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436368‑436663

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