Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,069	0	T	C	43.3%	80.5 / 32.8	68	G517G (GGA→GGG)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base T (31/7);  new base C (0/29);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.08e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGC  >  NC_000913/1575925‑1576214                                                                                                                                                 |                                                                                                                                                   tGCACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑at                                                                                                                                              <  1:1110319/150‑1 (MQ=255) tGCACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑at                                                                                                                                              >  1:755654/1‑150 (MQ=255)   cACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATc                                                                                                                                           >  2:201607/1‑151 (MQ=255)       cgcAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCa                                                                                                                                          >  2:1024663/1‑148 (MQ=255)         cAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAat                                                                                                                                     >  2:705100/1‑151 (MQ=255)           aTAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAatct                                                                                                                                   <  1:433412/151‑3 (MQ=255)            tAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATAtt                                                                                                                                  >  2:938358/1‑151 (MQ=255)               cGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTccc                                                                                                                               >  1:411224/1‑151 (MQ=255)                   cTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGcc                                                                                                                            >  1:174814/1‑150 (MQ=255)                    ttACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCAt                                                                                                                          >  2:2885/1‑151 (MQ=255)                      aCCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCAt                                                                                                                          >  2:18395/1‑149 (MQ=255)                       ccGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGt                                                                                                                       >  1:371365/1‑151 (MQ=255)                                      aaTCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACctg                                                                                                        >  2:896884/1‑151 (MQ=255)                                      aaTCCCATTCCAATTTCATTTTTAAAAGGCGTTTTAAAACCCTAAAACCGCTAAAGTTTTTTATTCCCTTCATCAGTTCCCGCCCCGTTTTAGCGAATATCACGTAACCCC‑AAAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACctg                                                                                                        <  1:912187/151‑1 (MQ=255)                                         cccATTCCAATTTCATCTTTATAAGGCTTTTTCAAATCTTGATAACGCGTCAGATTTTTATTCCCTGATCCCGTTTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACctgct                                                                                                      <  2:407503/150‑1 (MQ=255)                                         cccATTCCAATTTCATCGTTAAAGGGTTTTTTAAAACCTTAATACCCCGCCAGTTTTTTTTTCCCTAATCCCTATTCCGCCCAGTTATTGCGAATTCCACAAAACCCC‑ATACAAAAAATATCCCCGCCATAGTAACGGTTATAACctgct                                                                                                      <  2:255499/150‑1 (MQ=255)                                             ttGCCAGTTCATCTTTATAGGGCTTTTTAAAATCCTAAAAACGGGGAAGTTTTTTATTCCCTAATCCCGCTCCCGCCCACCAATTGCGAATATACCGTAACCCA‑ATAAAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTaa                                                                                                  <  1:1103237/150‑1 (MQ=255)                                              tGCCAGTTCATCTTTATAGGGTTTTTAAAACTCCTTAAAACGCGTCAGTTTTTTATTCCCTAACCCCGATCCCGCCCACCTATTGGGGAAATCAGATAACCCC‑ATAAAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                 <  2:976676/150‑1 (MQ=255)                                                ctgtttcttttttATAAGGTTTTTAAAATTCTTGAAACCGGGTGAGTTTTTTTTTCCCTGATCCCGATCCCCCCCAGCTTTTGCGTATATCAC‑TAACCCCAATAAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCa                                                                                               <  1:1155227/139‑1 (MQ=255)                                                ctatttCATCTTTATAGGGCTTTTTAAAATCCTGAAACCGCGTAAGATTTTTATTCCCTGATCCCGATCCCGTCCAGCTTTTGGGTATATCACTTAACCCC‑ATACAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGt                                                                                                    <  2:5720/141‑1 (MQ=255)                                                ccATTTCATCGTTATAAGGCTTTTTCAAATCCTGATAACGCGTCAGATTTTTTTTCCCTAATCCCGATCCCGCCCAGCTATTGCGGATATCACGTAACCCC‑ATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                 <  2:116211/148‑1 (MQ=255)                                                ccAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCAt                                                                                              >  2:84342/1‑151 (MQ=255)                                                    ttcttcTTTATAAGGTTTTTTAAAATCCTGATAACGCTTCAGATTTTTATTCCCTGATACCGATTCCGTCCACCTATTGGGTATATCACTTAACCCC‑ATAACAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                 <  2:740829/144‑1 (MQ=255)                                                    ttaaTCTTTAAAGGGCTTTTCAAAATCCGAAAACCGCGTAAGATTTTTATTCCCTGATCCCGATCCCGTCCAGTTTTTGGGGATATCCCGAAACCCA‑ATACAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATc                                                                                                <  2:912558/142‑1 (MQ=255)                                                        tttttATAAGGCGTTTTAAATTCTTGATACCGCGTGAGATTTTTATTCCCTGATCCCGTTCCCGCCCATCTTTTGCGAATATCACATACTCCA‑ATATCAAGAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAgg                                                                                         <  1:921842/145‑1 (MQ=255)                                                        tctttATAGGGCTTTTTAAAATCCTGATAACGCGTAAGTTTTTTATTCCCTGATCCCGATCCCGTCCAGTTTTTGGGAATATCACGTAACCCC‑ATATCAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGtt                                                                                       <  1:938554/147‑1 (MQ=255)                                                           ttATAAGGCTTTTTAAAATCCTAATAACGCGTAAGATTTTTATTCCCTAATCCCAATCCCGTCCAGTTATGGCGGATATCAC‑TAACCCCAATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGttt                                                                                      <  1:120078/148‑1 (MQ=255)                                                                gggCTTTTTAAACTCCTGATACGGCGTAAGTTTTTTATTCCCGGATCCCGTTCCCGCCCGGTTATGGGGGATATCAGGTACCCCA‑ATAAAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                               <  1:779885/149‑1 (MQ=255)                                                                  gCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaa                                                                             >  2:781924/1‑150 (MQ=255)                                                                     ttttAAAATCCTGATAACGCGTCAGATTTTTATTACCTGATACCGATTCCGTCCAGTTATGGCGGATATCACGTAACCCC‑ATATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaaa                                                                          <  2:411227/150‑1 (MQ=255)                                                                      tttCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAat                                                                         >  2:743452/1‑150 (MQ=255)                                                                        tCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATc                                                                      >  1:219655/1‑151 (MQ=255)                                                                          aaaTCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATccc                                                                    >  2:853566/1‑151 (MQ=255)                                                                            aTCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCAt                                                                  >  1:546583/1‑151 (MQ=255)                                                                                    accgctccAGATTTTTTTTACCTGATCCCGATTCGGTCCACTTATTGCGGATATCACTTAATCCC‑ATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCAtt                                                                 <  2:1194827/137‑1 (MQ=255)                                                                                    accgcgTAAGATTTTTTTTCCCTAATCCCAATCCCGCCCAGTTATGGCGAATATCACTTAACCCA‑AAATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCAt                                                          <  2:810446/149‑1 (MQ=255)                                                                                        cgTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATaa                                                        >  2:993567/1‑149 (MQ=255)                                                                                              agTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCggg                                                >  2:585285/1‑151 (MQ=255)                                                                                               tttttttTCCCTGATCCCGTTCCCGTCCAGTTATTGCGAATATCACGTAACCCC‑ATATCAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGcgcg                                                  <  2:877813/147‑1 (MQ=255)                                                                                               gTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgggg                                               >  2:291809/1‑151 (MQ=255)                                                                                                ttttttTCCCTAATCCCGATCCCGTCCAGTTATTGCGAATATCACTTACCCCC‑ATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgg                                                 <  2:689811/148‑1 (MQ=255)                                                                                                ttttATTCCCTGATCCCGATCCCGTCCATTTTTTGGGGATATCACGTAACCCC‑ATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAgcg                                                    <  2:374955/145‑1 (MQ=255)                                                                                                ttttATTCCCTAATCCCGATCCCGCCCAGTTTTTGGGGATATCACGAAACCCA‑ATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAgcg                                                    <  1:1095865/145‑1 (MQ=255)                                                                                                 tttATGTCCTGATCCCGTTTCGCTCCAGCTATGGCGGATATCGGTTACTCTA‑ATAGAAAAAAAACGCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATGCCGTCATAAAGCGCggg                                                <  2:642691/148‑1 (MQ=255)                                                                                                        ccTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATtgt                                      >  1:832208/1‑151 (MQ=255)                                                                                                        ccTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATtgt                                      >  1:584619/1‑151 (MQ=255)                                                                                                        ccTGATACCGATTCCGTCCAGCTATGGCGAATATCACGAAAACAC‑ATAACAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATAtt                                        <  1:53743/149‑1 (MQ=255)                                                                                                             tcccGATTCCGCCCAGTTTTTGGGAATATCCCGAATCCCC‑ATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATAtt                                        <  2:1126254/142‑1 (MQ=255)                                                                                                               ccttttCCGTCCAGTTATGGCGAATATCACGTAACCCA‑ATATCAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                <  2:667353/146‑1 (MQ=255)                                                                                                               cctttccctcccAGTTTTTGCGAATATCACAAAACCCC‑ATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGttt                                 <  1:1111139/139‑1 (MQ=255)                                                                                                               ccGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                >  1:172299/1‑150 (MQ=255)                                                                                                               ccGATCCCTCCCAGTTATTGCGAATATCACGTAACCCA‑ATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATtgtg                                     <  1:111472/145‑1 (MQ=255)                                                                                                               ccGATCCCGTCCAGTTATTGGGAATACCACGTACCCCC‑ATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTgg                                    <  2:1209588/146‑1 (MQ=255)                                                                                                               ccGATCCCGCCCAGTTATTGGGAATATCACGTAACCCC‑ATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGtt                                  <  1:201607/148‑1 (MQ=255)                                                                                                                   ttCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                              >  2:135819/1‑148 (MQ=255)                                                                                                                   ttCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                              >  1:42271/1‑148 (MQ=255)                                                                                                                         ccAGTTATGGCGGATATCACGTAACCCC‑ATAAAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAg                        <  1:513766/148‑1 (MQ=255)                                                                                                                         ccAGCTATTGCGGATATCACGTAATCCC‑ATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                       <  2:820877/149‑1 (MQ=255)                                                                                                                         ccAGCTATGGCGGATATCACATAACCCC‑ATAACAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGtt                      <  1:2885/150‑1 (MQ=255)                                                                                                                           aGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtc                   >  2:359949/1‑151 (MQ=255)                                                                                                                              tATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcat                 >  1:850931/1‑150 (MQ=255)                                                                                                                                tggcgaatcccccGTAACCCC‑ATAAAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatag               <  2:386568/138‑1 (MQ=255)                                                                                                                                tggCGAATATCACGAACCCCA‑ATACAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtc                   <  1:1167631/144‑1 (MQ=255)                                                                                                                                  gggaatacCACGTACCCCC‑AAAACAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcat                 <  1:208042/142‑1 (MQ=255)                                                                                                                                      atatCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATAcc        >  1:714042/1‑151 (MQ=255)                                                                                                                                         tCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGc      >  1:119079/1‑150 (MQ=255)                                                                                                                                           aCGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACg   >  1:152534/1‑151 (MQ=255)                                                                                                                                            cGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGc  >  1:1208862/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TGCACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGC  >  NC_000913/1575925‑1576214

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