Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,859,685	0	A	C	42.9%	78.8 / 17.9	57	V182G (GTA→GGA)	ydjK	putative MFS sugar transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (23/9);  new base C (0/24);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 9.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTATGACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGA  >  NC_000913/1859536‑1859823                                                                                                                                                        |                                                                                                                                           ggCGTATGACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTa                                                                                                                                            <  2:96059/150‑1 (MQ=255) ggCGTATGACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTa                                                                                                                                            <  2:391551/150‑1 (MQ=255)   cGTATGACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTAccc                                                                                                                                         <  1:314438/151‑1 (MQ=255)      aTGACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGcc                                                                                                                                      >  1:766280/1‑151 (MQ=255)        gACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGcc                                                                                                                                      >  1:576040/1‑149 (MQ=255)                 ttctATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCActct                                                                                                                           >  2:897669/1‑151 (MQ=255)                    tataCTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGc                                                                                                                         >  1:1318/1‑150 (MQ=255)                        cTCCGCATCTCTTTCTCGCCTTCCTGATCCCGCCCGCGCGATCCCCCCCGGCGGGGGGTTTAAGA‑AAAGTAGCGCCGGCCCACCGCCGTAGCAAAAAGCGACAATATTGCGGGA‑ATAAGCATTTGTCCCCGCCAGTTCCACTCTGCAc                                                                                                                       <  1:123913/148‑1 (MQ=255)                          ccGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATc                                                                                                                  >  2:944580/1‑151 (MQ=255)                            gaaTGACTTTCTCGCCTTCGTGTTCCCGGCCGCGCTTCCCCCCCGAGGGCGGGAATTAAGA‑AAAGTGGCGCCAGCCCGGCGCTGGAGTGATAAGCAAAAAAGTGGAAGGG‑ATAGGTATTGTTACCCGCCAGTTCCACTCTGCACTAATCAg                                                                                                                <  2:4596/149‑1 (MQ=255)                               tGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAgcgg                                                                                                             >  2:636420/1‑151 (MQ=255)                                gtcTTTCTCGCCTCCCTGACCCCGCCCGCGCGTTTCCAGCCAGCGGGGGGTTTAAGG‑AAGGTGGCGCAGGGCCGCCGCCGTAGCGAAAAGCAAAAATGTGGCGGGA‑AAAAGCAGTGTTCCCCGCCAGTTCCACTCTGCACTAATCAgcggc                                                                                                            <  1:867335/149‑1 (MQ=255)                                gtcTTTCCCGGCTCCCTGATCCCGCCCGCGCGATTCCAGCCAGCGGGGGGATTCAGG‑AAAGTAGCGCCAGGCCGGCGCCGTAGCGATAAGCGAAAATATGGCGGGA‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAgc                                                                                                               <  1:993180/146‑1 (MQ=255)                                   tttCCCGGCTCCCTGACCCCGCCCGCGCGATCCCGCCCGGCGGGGGGATTAAGA‑AAAGAGGCCCCAGGCCGGCGCGGTAGCGAAAAGCGACAATATGGCAGGA‑ATAAGCAGTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGt                                                                                                          <  1:130399/150‑1 (MQ=255)                                          gCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTccc                                                                                                   >  1:30927/1‑150 (MQ=255)                                                 gttcccccccGCGCGTTCCCACCCAGCGGGGGGTTTAAGA‑AAAGGAGCGCCGGCCCGGCGCCGTGGCAAAAAGCAAAATTTTGGGGGGG‑ATAAGCTGTTTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTAt                                                                                           <  1:658713/143‑1 (MQ=255)                                                    ccccccggcgcaATCCCAGCCGGCGGGGTGATCAAGA‑AAAATGGCCCCAGCCCAGCGCGGGAGCAAAAAGAAACAATATGGCGGGG‑ATAACAAGTTTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCa                                                                                         <  1:490887/143‑1 (MQ=255)                                                    aCCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCa                                                                                         >  2:878490/1‑150 (MQ=255)                                                     ccgcccgcGCGCTCCCAGCCGGCGGGGGTATTGAGA‑AAAGTGGCGCCGGCCCGGCGCCGGAGCAAAAAGCGAAAATATGGCGGGG‑AAAGGAAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCa                                                                                         <  2:546366/149‑1 (MQ=255)                                                         cggcgcgTTCCCAGCCGGGGGGGGGATTAAGA‑AAAGTAGCGCAGGGCCGGCGCGTTGGCAAAAGGCGAAAATATTGAGGGA‑AAAAGCATTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAAc                                                                                   <  1:221994/149‑1 (MQ=255)                                                         ccgggggttcccagccggGGGGGGGATTAAGA‑AAAGTAGCGCCGGGCCGGCGCGGTAGCAATAAGCAAAAATATTGCAGGA‑ATAAGCATTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAAc                                                                                   <  1:715456/141‑1 (MQ=255)                                                         ccGCGCGATCCCAGCCGGGGGGGGGATGAAGA‑AAAGGAGCGCCAGGCCGGCGCCGTAGCAATAAGCGAAAATATGGCGGGG‑AAAGGGGGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAAc                                                                                   <  2:621632/151‑1 (MQ=255)                                                          cgcgcgATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAAc                                                                                   >  1:749010/1‑150 (MQ=255)                                                          cgcgcgATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAAc                                                                                   >  1:262012/1‑150 (MQ=255)                                                               gATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCg                                                                             >  1:697601/1‑151 (MQ=255)                                                                    cagccggcGGGGTAATTAAGA‑AAGGTAGCGCCGGGCGAGCGCCGGAGCAATAAGCAAAAATATGGGGGGG‑ATAAGCATTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACagag                                                                               <  2:1317/144‑1 (MQ=255)                                                                      gccagcGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATg                                                                       >  2:721008/1‑150 (MQ=255)                                                                       ccgggggggggATAAAGG‑AAAGGGGCGCCGGCCCGGCGCGGTGGAGAAAAGCAAAAATATGGAGGGG‑AAAAGCGGTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGa                                                                           <  1:564500/138‑1 (MQ=255)                                                                          gcggGGGGATAAAGG‑AAAGTCGCCCCGGCCCAGCCCGGTGGGAAAAAACAAAAATATGGGAGGG‑AAAAGCATTTGGACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAg                                                                  <  2:933766/151‑1 (MQ=255)                                                                             gggggAATCAAGAAAAGGAGCGCCGGGCCGGCCCCGTAGCAAAAAGCAAAAATTTTGAAGGA‑AAAGGCATTTTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCa                                                                   <  1:769720/147‑1 (MQ=255)                                                                             cgGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTg                                                               >  1:459276/1‑151 (MQ=255)                                                                              ggggATTCAGA‑AAAGTCCCCCCAGCCCAGCGCGGTAGCAATAAGCAAAAATTTGGCGGGA‑AAAAGCGTTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTg                                                               <  1:392999/150‑1 (MQ=255)                                                                              ggggATAAAGA‑AAGGTAGCCCCGGCCCGGCGCGGTGGCAATAAGCAAAAATTTTGGGGGA‑AAAAGCAGTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTg                                                               <  1:351387/150‑1 (MQ=255)                                                                                    ttaagaaaaggggcgcCAGGCCGGCGCCGTAGCGAAAAGCAAAAATATGGCGGGA‑AAAAGCAGTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGTGGATATGACCAGTTGCCAAt                                                          <  1:579111/138‑1 (MQ=255)                                                                                      aGG‑AAAGTAGCGCCGGCCCGGCGCAGTAGCAAAAAGCGAAAATTTGGCAGGG‑TTAAGCGTTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATa                                                         <  2:190660/148‑1 (MQ=255)                                                                                        gaaaaGTGGCCCCGGCCCGGCGCCGGAGCAAAAAGCGA‑AAAACTGGAGGGAAAAAGAGGTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAg                                                      <  2:965381/148‑1 (MQ=255)                                                                                          aaaggggcgcCGGCCCAGCGCTGTGGCAATAAGCAAAAATATGGGAGGG‑AAAAGCAGTTTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATaaa                                                       <  2:417551/141‑1 (MQ=255)                                                                                          aaaGTGGCGCCGGGCCGGGGCTGAGGGGAAAAGCAAAAATGTGGAAGGG‑AGAAGTGGTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCaa                                                           <  2:347398/143‑1 (MQ=255)                                                                                           aaGTGGCGCCAGGCCAGCGCGGTAGAGATAAGCAAAAAAATGGGAGGG‑ATAGGCAGTGTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGaaa                                                  <  2:50723/151‑1 (MQ=255)                                                                                            aggaGCGCCGGCCCGGCGCGGTGGAGAAAAGCAAAATTATTGCGGGA‑ATAGGCGGTTGCGCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGTGGATATGACCAGTTGCCAATAAAGGAAAc                                                 <  1:999507/148‑1 (MQ=255)                                                                                                 cgcCGGGCCGGCGCCGGAGCGATAAGCAAAAAT‑TTGGAGGGAAAAAGCGTTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGa                                             <  2:94624/150‑1 (MQ=255)                                                                                                     aggccagCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTg                                         >  2:397085/1‑150 (MQ=255)                                                                                                               cGTAGCGATAGCCGAAAATATGGCGGGG‑ATAACCTGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCa                              <  1:595235/151‑1 (MQ=255)                                                                                                                gtggCGAAAAGCAACAATATGGAGGGG‑AAAGGCAGTGTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCa                              <  2:476368/147‑1 (MQ=255)                                                                                                                  aGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACg                            >  2:650361/1‑150 (MQ=255)                                                                                                                   gaaaaaaGCAAAAATATTGAGGGA‑AAAAGCATTTGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACg                            <  1:659199/144‑1 (MQ=255)                                                                                                                    cGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGAtt                         >  1:913902/1‑151 (MQ=255)                                                                                                                       tAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCt                      >  1:188889/1‑151 (MQ=255)                                                                                                                        aaGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCt                      >  2:400863/1‑150 (MQ=255)                                                                                                                             aaaaTATGGCAGGG‑AAAGGCATTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCgg                 <  1:42439/148‑1 (MQ=255)                                                                                                                               aaTATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGc               >  1:409007/1‑150 (MQ=255)                                                                                                                               aaTATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGc               >  1:183271/1‑150 (MQ=255)                                                                                                                                 tatTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCat             >  2:484939/1‑150 (MQ=255)                                                                                                                                 tatTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCat             >  1:697839/1‑150 (MQ=255)                                                                                                                                     ggcgggaaaaaGCTGTTGCCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCata            <  2:154188/139‑1 (MQ=255)                                                                                                                                        ggg‑AAAAGCTTTTTTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATAtt        <  1:115918/148‑1 (MQ=255)                                                                                                                                          gaataggcttttGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCgg     <  2:193417/140‑1 (MQ=255)                                                                                                                                             tAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGa  >  1:489086/1‑150 (MQ=255)                                                                                                                                                gcgtttGTCCCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGt    <  1:465994/141‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                           GGCGTATGACGCCTTCTTCTATACTCCGCATGACTTTCTCGGCTTCCTGATACCGTCCGCGCGATTCCAGCCAGCGCGGGGATTCAGG‑AAAGTAGCGCCAGGCCAGCGCCGTAGCGATAAGCGACAATATTGCAGGG‑ATAAGCAGTTGTACCCGCCAGTTCCACTCTGCACTAATCAGCGGCGTGAGTCCCATCGCTATCAATGAACAGAGCGGATATGACCAGTTGCCAATAAAGGAAACCCGACTTGACCACGTTCCACGATTTCTACCGGGCATATATTCGGTGA  >  NC_000913/1859536‑1859823

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: ‑